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ABSTRACT 
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The  levels  m  Ar  below  3.30  MeV  excitation 

37  37 

have  been  studied  using  the  Cl(p,n)  Ar  and  the 
37  37 

Cl(p,ny)  Ar  reactions.  These  levels  were  identified 
from  neutron  time-of-f light  spectra  and  their  energies 
from  deexcitation  y-ray  spectra.  A  level  at  3185  keV 
excitation  has  been  positively  confirmed.  The  mean  life¬ 
times  of  these  levels  were  determined  by  the  Doppler- 
shift  attenuation  method  with  the  y  rays  measured  in 
coincidence  with  neutrons  and  in  singles.  Gamma- ray 
angular  distributions  were  obtained  at  various  incident 
energies  and  compared  with  the  statistical  compound 
nuclear  (SCN)  model.  Spins , parities  and  mixing  ratios 
have  been  deduced  in  conjunction  with  transition  strengths 
based  on  Weisskopf  estimates.  Neutron  differential  cross 
sections  were  also  measured  using  the  time-of-f light 
method  and  compared  with  the  statistical  (Hauser-Feshbach) 
calculations . 


The  spin  and  parity  assignments  and  mean  lifetimes 

3  7 

have  been  extracted  for  the  following  levels  of  Ar: 

1410  keV  (J77  =  l/2+ ,  x  =  943*iq0  fs)  ,  1611  keV  (J77  =  7/2", 
T>3000  fs),  2217  keV  (J77  =  7/2+,  x  =  510*^  fs)  '  2488  keV 
(J77  =  3/2",  x  =  660^°  fs),  2795  keV  (J77  =  5/2+,  x  = 

fs),  3171  keV  (J77  =  5/2+,  x  =  85^4°  fs)  ,  3185  keV 


-12 


(J77  =  9/2  ,  x  =  285*™  fs)  and  3274  keV  (J^  =  5/2^  x  = 

+  21 

50_^g  fs) .  The  results  of  the  present  investigation  are 


IV 


compared  with  previous  experimental  and  theoretical 
works . 

29  29 

A  study  of  the  Si(p,p'y)  Si  reaction  was 
undertaken  to  investigate  statistical  limitations  of 
the  SCN  model.  Using  a  thin  target  angular  distribu¬ 
tions  of  y-rays  deexciting  from  the  1st  and  2nd  states 
29 

of  Si  were  measured.  Proton  energies  were  varied 
from  3.00  to  3.60  MeV  in  10  keV  intervals.  Angular 
distributions  for  progressively  thicker  targets  were 
constructed  from  the  thin  target  yields  and  compared 
with  predictions  of  the  SCN  model.  An  estimate  of  the 
number  of  levels  required  for  successful  application 
of  the  model  was  thus  obtained. 
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CHAPTER  1 


INTRODUCTION 

Nuclei  near  doubly-closed  shells  are  of  consider¬ 
able  interest  due  to  the  relatively  simple  shell  model 
configuration  expected  and  because  experimental  data  can 

be  readily  compared  in  detail  with  theoretical  predic- 

3  7 

tions.  The  Ar  nucleus  which  lacks  a  neutron  and  two 

40 

protons  from  the  'magic'  Ca  nucleus  is  of  particular 

interest  because  its  low-lying  excited  states  extend  into 

the  next  shell;  hence  states  of  both  parities  are 

observed.  Its  odd  mass  can  also  permit  studies  involving 

models  of  coupling  a  single  particle  to  an  even-even  core. 

In  addition  to  its  interest  from  a  nuclear  struc- 

37 

ture  point  of  view,  the  Ar  nucleus  also  has  astrophy- 

sical  implications.  It  has  been  reported  that  the 

3  7  -  3  7 

detection  of  solar  neutrinos  by  the  Cl(ve,e  )  Ar 

3  7 

reaction  involves  some  excited  states  in  Ar  (Ba  64, 

Da  68  and  Ba  68) .  Thus  knowledge  of  the  low-lying  levels 
could  lead  to  better  understanding  of  the  solar  energy 
generation  process. 

At  the  time  when  this  experiment  was  initiated 

37 

most  of  the  experimental  information  on  the  Ar  nucleus 
had  been  obtained  by  charged  particle  spectroscopy  (En  67 
and  Na  68) .  As  a  result,  excitation  energies  contained 


1 


2 


large  uncertainties  (±15  keV)  and  electromagnetic 
properties  were  incomplete;  only  the  lifetime  of  the 
1611  keV  level  had  been  reported  (Go  67) . 

The  present  investigation  was  undertaken  to 

provide  a  more  complete  experimental  description  of 

3  7 

the  low-lying  states  of  Ar.  One  of  the  intentions 
was  to  determine  more  precisely  the  energies  of  the 
excited  states  in  order  to  separate  perhaps  previously 
unresolved  states.  The  study  of  y-ray  deexcitation  of 
these  states  could  then  lead  to  information  of  their 
decay  modes,  branching  ratios  and  spins  and  parities 
(J  ) .  Also  measurement  of  Doppler-shifted  y  rays  could 
yield  mean  lifetimes  for  the  excited  states  as  well  as 
infer  limits  on  their  J  values.  Several  different 
experimental  observations,  however,  were  necessary  to 
obtain  the  desired  information  and  these  will  be  dis¬ 
cussed  separately. 

The  37Cl(P,n)37Ar  reaction  (Q=  -1.596  MeV)  was 
chosen  for  this  investigation  for  several  reasons.  The 
(p,n)  reaction  would  populate  all  energetically  allowed 
states  and  its  negative  Q-value  would  also  permit  one  to 
excite  only  up  to  the  levels  of  interest  unlike  the  (d,p) 
reaction  for  example.  Furthermore,  experimental  infor¬ 
mation  on  this  reaction  had  been  lacking;  neutron  emission 

threshold  energies  of  the  lowest  four  excited  states  in 
3  7 

Ar  and  the  angular  distribution  of  the  ground  state 
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neutron  group  had  been  the  only  observations  reported 
(Fe  59  and  Iy  66).  Finally  at  low  incident  proton 
energies  (<6  MeV)  the  (p,n)  reaction  would  proceed 
mainly  through  the  compound  nuclear  mechanism  and 
experimental  conditions  could  be  suitably  arranged  to 
justify  the  statistical  analysis.  Such  results  are 
necessary  for  a  comparison  with  those  obtained  using 
the  direct  interactions  or  model  independent  approaches 
to  provide  a  better  understanding  of  the  SCN  model. 

The  lifetimes  of  all  the  levels  up  to  3274  keV 
3  7 

m  Ar  have  been  studied  using  the  Doppler  shift  atten- 

3  7 

uation  method  (DSAM) .  Enriched  Cl  targets  were 
bombarded  by  5.70  MeV  protons  from  the  van  de  Graaff 
accelerator  at  the  Nuclear  Research  Centre  of  the 
University  of  Alberta.  The  Doppler-shifted  y  rays  from 
the  excited  states  were  observed  using  a  Ge(Li)  detector 
at  0°,  90°  and  120°  with  respect  to  the  beam  direction. 
The  spectra  were  collected  in  singles  and  in  coincidence 
with  neutrons  detected  at  180  degrees  by  an  annular 
detector.  From  the  energy  shifts  of  the  detected  y  rays, 
mean  lifetimes  could  be  determined  according  to  the 
method  of  Warburton  (Wa  67).  The  description  of  the 
experiment  and  results  are  contained  in  Chapter  2. 

Chapter  3  reports  the  y-ray  yields  and  angular 
distributions  measured  with  different  target  thickness 
and  incident  energies.  The  results  were  compared  with 
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the  statistical  compound  nuclear  (SCN)  model  to  predict 
the  spins  and  parities  of  the  level  and  to  find  multipole 
mixing  ratios  of  their  deexcitation  y  rays;  lifetime 
results  were  utilized  to  restrict  some  of  these  assign¬ 
ments.  A  comparison  of  the  SCN  model  and  model-indepen¬ 
dent  analysis  was  also  made  for  the  1611  keV  level. 
Branching  ratios  were  extracted  from  the  zero-order  term 
of  the  Legendre  polynomial  expansion  of  the  angular  dis¬ 
tributions  . 

Since  the  reaction  was  assumed  to  proceed  mainly 

through  the  compound  nuclear  (CN)  process,  it  was  of 

interest  to  measure  the  neutron  angular  distributions 

in  search  of  symmetry  around  90°  in  the  center-of-mass 

system  which  is  characteristic  of  the  CN  interaction 

(En  66).  Neutron  angular  distributions  also  provided 

3  7 

knowledge  of  the  average  Ar  recoil  velocity  which  was 
helpful  in  our  DSAM  singles  analysis.  Of  course,  neutron 
spectra  would  also  show  which  excited  states  are  populated, 
thereby  aiding  the  proper  identification  of  the  y  rays. 

This  resulted  in  the  confirmation  (from  the  neutron  spec¬ 
trum)  of  the  existence  of  a  new  level  at  3185  keV  which 
did  not  decay  to  the  ground  state.  The  time— of —flight 
technique  was  used  to  record  the  neutron  spectra  and  is 
presented  in  Chapter  4. 

3  7 

The  summary  of  the  Ar  nucleus  and  comparison  with 
Shell  model  calculations  are  contained  in  Chapter  5. 
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In  the  analysis,  the  SCN  model  and  model- 

independent  calculations  have  been  compared.  The 

agreement  is  good;  however,  it  was  felt  that  a  more 

detailed  investigation  of  the  applicability  of  the  SCN 

model  was  desirable.  An  extensive  angular  distribution 

29  29 

study  using  the  Si(p,p'y)  Si  reaction  was  therefore 

undertaken.  The  selection  of  this  reaction  and  the 

result  of  the  investigation  are  described  in  Chapter  6. 

Basically,  the  experiment  consists  of  measuring  y-ray 

angular  distributions  with  a  thin  target  over  a  range  of 

bombarding  energies,  in  small  energy  increments.  Hence 

variations  of  the  angular  distribution  as  a  function  of 

yield  and  target  thickness  may  be  studied  to  afford  an 

insight  into  the  statistical  process. 

During  the  course  of  this  project  several  inves- 
3  7 

tigations  on  the  Ar  nucleus  have  been  published  includ¬ 
ing  part  of  this  dissertation  (Wo  72  and  Appendix  E) . 
Inspite  of  the  numerous  lifetime  determinations  of  the 
low-lying  states,  agreement  has  not  been  satisfactory. 
The  level  energies,  decay  modes  and  J77  assignments  have, 
however,  agreed  more  consistently.  Detail  comparison  of 
these  and  present  results  will  be  made  in  the  discussion 
of  the  respective  subjects. 

Several  calculations  within  this  study  require 
modification  of  standard  formula  and  computer  codes  to 
expedite  calculations.  A  short  description  of  these 
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together  with  published  papers  are  presented  in  the 
Appendices . 

Appendix  A  describes  modification  of  the  Doppler- 

shift  attenuation  factor  to  include  both  the  direct 

(particle)  and  cascade  (y  ray)  feeding  of  the  state. 

A  program,  ’MAGPOP',  written  for  the  model  independent 

3  7 

study  of  the  1611  keV  level  in  the  Ar  nucleus  (Section 
3.4.2),  is  presented  in  Appendix  B.  In  Chapter  6,  the 
procedure  to  simulate  yields  of  thicker  targets  by  summing 
and  averaging  yields  from  a  thin  target  is  described  by 
the  code,  'WAVY',  in  Appendix  C. 

40 

A  paper  on  the  positive  parity  states  m  K  which 
introduced  the  author  to  (p,ny)  studies  and  the  SCN  model 
is  contained  in  Appendix  D.  Finally,  part  of  this  dis¬ 
sertation  which  has  been  published  is  included  in 
Appendix  E  and  the  listing  of  the  code  ’ 0XDS’  (Ro  69)  is 
given  in  Appendix  F. 


CHAPTER  2 


MEAN  LIFETIMES  OF  3?Ar 


2 . 1  Doppler  Shift  Attenuation  Method 

Electromagnetic  transition  strength  is  usually 
expressed  in  terms  of  mean  lifetime,  x,  for  the  emission 
of  a  y  ray  of  a  particular  multipolarity  between  a  pair 
of  nuclear  energy  levels  or  in  terms  of  corresponding 
radiative  width,  T,  according  to  the  uncertainty  princi¬ 
ple,  Tx  =  fu  The  measured  strength  of  the  transition 
can  then  give  important  clues  as  to  the  properties  of 
the  nuclear  energy  levels  involved. 

The  Doppler  shift  attenuation  method  (DSAM)  (Wa  67a) 

has  been  adopted  in  the  present  investigation  to  determine 

3  7 

lifetimes  of  the  Ar  levels.  This  method  relies  on  the 
rapid  slowing  down  of  swiftly  moving  atoms  in  a  solid. 

In  short,  y  rays  emitted  by  an  excited  nucleus  in  motion 
will  exhibit  the  familiar  Doppler  shift  in  energy.  The 
average  (centroid)  energy  of  the  y  rays  emitted  by  recoil¬ 
ing  nuclei  of  varying  velocities  is  given  by  (Wa  63  and 


01  66) 


v 

E  (0)  =  Ey{)[l  +  F(T)  cos  6] 


(2.1) 


where  E^Q  specifies  the  y-ray  energy  emitted  by  nuclei 
at  rest  and  0  is  the  angle  between  the  gamma  emission  and 
the  bombarding  direction.  vq  is  the  average  component  of 
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the  initial  recoil  velocity  along  this  axis  and  c  is 
the  velocity  of  light.  F(x)  is  the  Doppler-shif t 
attenuation  factor  whose  deviation  from  unity 
(0  <  F (t)  <  1)  describes  the  effect  on  E^(0)  due  to 
the  slowing  down  process  of  the  recoiling  nuclei. 

This  attenuation  factor  may  be  obtained  experi¬ 
mentally  by  observing  the  average  y-ray  energies  at  two 
detection  angles,  0^  and  0£.  From  the  energy  difference 
between  these  measurements  one  gets. 


F(t) 


E  (0,  )  -  E  (0O) 
y  1  y  2 


ey0  (cos  ei  “  cos 


(2.2) 


The  denominator  may  be  recognized  as  the  maximum  energy 
difference  (F(i)  =  1)  and  it  can  be  calculated  from  re¬ 
action  kinematics.  Thus  one  rewrites  Equation  2.2 


F(t) 


[ey<V  '  Ve2)] 

lVV  ■  Vb2)J 


observed 

maximum 


(2.3) 


3  7  3  7 

In  the  present  experiment,  Cl(p,n)  Ar,  the 

3  7 

initial  recoil  velocity  of  the  Ar  nuclei  was  restricted 
to  the  forward  direction  (in  the  same  direction  as  the 
incident  protons)  by  requiring  the  y  rays  to  be  in  coin¬ 
cidence  with  the  neutrons  detected  in  the  backward 
direction.  This  requirement  ensured  maximum  initial 
recoil  velocity  and  hence  maximum  Doppler  shift  according 
to  Equation  2.1.  However,  for  the  higher  excited  states 
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(which  had  low  yields) ,  y-ray  shifts  were  measured  in 
singles  only.  The  mean  initial  recoil  velocity  in  this 
case  was  derived  from  the  independently  measured  neutron 
angular  distributions  (Chapter  4)  as  well  as  from  the 
method  of  Warburton  (Wa  6  7)  .  For  the  lower  energy 
levels,  the  lifetimes  have  been  obtained  by  both  coinci¬ 
dence  and  singles  methods. 

The  theoretical  derivation  of  the  attenuation 
factor,  F (t)  has  been  described  extensively  (Wa  63,  B1  66 
Ro  69  and  Gu  70)  and  hence  only  a  brief  outline  will  be 
given  here.  The  F(t)  that  was  given  by  Warburton  et  al 
(Wa  63)  was  modified  by  Blaugrund  (B1  66)  to  include  the 
cos  (p  (t)  term  and  can  be  written  as 


F  ( i ) 


v  ( t ) 


cos  ( t ) 


(2.4) 


where  t  is  the  lifetime,  vq  and  v(t)  are  the  initial  and 

instantaneous  velocities  respectively,  ( t )  is  the  angle 

-t/r 

between  v  and  v(t)  and  e  '  is  the  probability  of  decay 
o 

The  term,  cos~~^TtT,  was  introduced  to  account  for  the 
large  angle  scattering  of  the  recoiling  nucleus  and  can 
be  calculated  from  multiple  scattering  consideration 


(Le  50)  . 


cos  cp  =  exp(- 


x 


0 


K  .  dx) 

<P 


(2.5) 


in  which  x  is  the  distance  travelled  in  the  stopping 


medium  and 
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K  ,  =  N 

4> 


(1  -  cos  (J))  da((j)) 


(2.6) 


<J> 


where  N  represents  the  number  of  scatterers  per  cm  and 
da(<j))  the  differential  scattering  cross  section  of  the 
recoil . 

It  has  been  shown  (Po  61)  that  the  velocity 

function  v(t)  in  Equation  2.4  varies  linearly  with 

dE 

stopping  power,  ^  ,  of  the  stopping  material  for  the 
recoiling  nuclei  and  can  be  expressed  in  terms  of  it. 

The  stopping  power,  however,  may  be  divided  into  elec¬ 
tronic  and  nuclear  stopping  components.  The  electronic 
process  involves  energy  loss  of  the  recoiling  ion  due 
to  collision  with  electrons  in  the  stopping  material  and 
the  nuclear  process  involves  collision  with  atoms  in  the 
stopping  medium.  According  to  Lindhard  et  al  (Li  61  and 


Li  63) 


(dE) 


electronic 


k  (— ) 
K  'M 


1/2 


keV. cm  /mg 


(2.7) 


where  E  is  the  energy  of  the  moving  ion,  M  is  its  mass 


and 


k  = 


73.0  Z^/6  Z1  Z2 


M. 


(z3/3  +  z2/3)3^2 


.  .  e"  „2/3  1 

for  v  <  -3r  Z ,  =  — 

n  1 


_  cz2/3 

137  CZ1 


(2.8) 


where  Z  and  M  are  the  charge  and  mass  and  the  subscripts, 
1  and  2,  denote  the  moving  ion  and  the  stopping  material 

respectively. 
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In  similar  notation,  the  nuclear  stopping  is 


given  by 


.  =  K  Bn  ln 

nuclear 


2  E  M  'h 


2 


(Mi+M2 )  Z±  Z 2  e^  {Z^/3+Z* 2 32/3)1/2 


(2.9) 


where 


B 

n 


2 


M 


2 


E 


4 

e 


The  total  stopping  power  is  simply  the  sum  of 
Equations  2. 7  and  2.9.  Calculations  of  the  theoretical 
F(x)  from  the  total  stopping  power  were  performed  using  a 
computer  code  ’OXDS'  (Ro  69  and  Appendix  F) .  The  program 
has  provision  to  vary  the  electronic  and  nuclear  contri¬ 
butions  to  the  stopping  power.  Thus  mean  lifetimes  may 
be  obtained  from  comparison  of  the  observed  attenuation 
factor  and  the  theoretical  F(t)  vs  t  curves. 


2.2  Experimental  Details 

3  7 

Lifetimes  of  the  low-lying  levels  in  '  Ar  below 

3.274  Mev  excitation  have  been  determined  by  the  Doppler- 

shift  attenuation  method.  These  levels  were  populated 
3  7  3  7 

via  the  Ci(p,ny)  Ar  reaction  using  a  5.70  MeV  proton 

2 

beam.  The  target  consisted  of  0.5  mg/cm  of  PbCl^ 

3  7 

(enriched  to  87%  in  Cl)  evaporated  onto  a  0 . 5  mm  thick 
Au  backing.  Gamma  rays  were  detected  both  in  coincidence 

3 

with  neutrons  and  in  singles  using  a  45  cm  Ge(Li)  y-ray 


■ 


detector  placed  10  cm  from  the  target  at  0°,  90°  and 
120°  with  respect  to  the  incident  beam  direction.  The 
detector  provided  a  resolution  of  4  keV  at  1.33  MeV 
under  experimental  conditions.  Neutrons  were  detected 
in  a  5.1  cm  thick  annular  (17.8  cm  O.D.,  25  cm  I.D.) 

NE  218  liquid  scintillation  counter  which  subtended 
angles  between  157°  and  177°  to  the  beam  direction. 

It  was  placed  15.8  cm  from  the  target  (Figure  2.1). 

The  y  flux  in  the  neutron  detector  was  suppressed  by  a 
1.95  cm  thick  lead  shield. 

In  order  to  detect  only  the  y  rays  emitted  by  the 
3  7 

excited  states  of  Ar,  stringent  neutron  coincidence 
restrictions  were  imposed  on  the  Ge(Li)  signals.  These 
requirements  consisted  of  distinguishing  neutrons  and  y 
rays  arriving  at  the  neutron  detector  by  their  time-of- 
flight  separation  and  also  by  their  rise  time  differences 
in  the  detector.  A  detailed  description  of  the  n-y 
discrimination  and  all  electronic  circuitry  will  follow 
after  a  brief  outline  of  the  general  experimental 
arrangement . 

Gamma-gamma  coincidence  between  the  Ge(Li)  detec¬ 
tor  and  a  7.5  cm x  7.5  cm  Nal  (Tl)  y-ray  detector  was  also 
arranged.  This  allowed  y  rays  from  known  sources  placed 
nearby  to  be  recorded  simultaneously  during  the  experi¬ 
ment  to  serve  as  internal  calibration  lines  as  well  as  a 
continuous  monitor  of  the  system  stability. 
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Figure  2.1  Experimental  arrangement  of  target  and  detectors 
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The  n-y  and  y-y  coincidence  signals  were  accumu¬ 
lated  in  parallel  and  added  so  that  a  single  analog-to- 
digital  converter  (ADC)  was  used.  The  advantage  of 
using  one  ADC  for  all  the  y  rays  detected  by  the  Ge(Li) 
whether  they  satisfied  the  n-y  or  the  y-y  coincidence 
requirement,  was  that  all  the  y  rays  would  then  be 
subjected  to  the  same  conversion  gain.  Thus  the  known 
y-ray  peaks  from  the  sources  may  be  used  directly  as 
internal  calibration  lines  without  gain  correction  aris¬ 
ing  from  the  use  of  another  ADC.  The  n-y  and  y-y 
coincidences  were,  however,  separated  into  different 
regions  of  the  computer  memory  by  using  the  y-y  coinci¬ 
dence  signals  as  routing  pulses. 

A  block  diagram  of  the  electronics  is  shown  in 
Figure  2.2  and  its  description  follows. 

Neutron  and  y-ray  timi ng  coincidence 

The  Ge(Li)  detector  has  two  output  signals;  an 
energy  and  a  timing  signal.  A  combination  of  timing 
filter  amplifier  (TFA)  and  constant  fraction  pulse  height 
trigger  (CFPHT  1)  were  used  to  obtain  the  starting  time 
signal  for  the  time-to-amplitude  converter  (TAC  1) .  The 
stop  pulses  were  obtained  from  the  fast  output  of  the 
neutron  detector  with  the  use  of  another  (CFPHT  2) .  Two 
gate  and  delay  generators  (G&D  1  and  2)  were  employed  to 
sychronize  the  start  and  stop  signals  so  that  they  fall 
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Figure  2.2  Block  diagram  of  electrons 
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within  the  time  range  (50  ns)  of  the  (TAC  1) .  The 
different  velocities  between  the  neutrons  and  the  y  rays 
caused  their  arrival  at  the  liquid  scintillator  to  be 
separated  by  their  flight  time  difference  and  this 
yielded  essentially  two  groups  of  pulses  at  the  output 
of  the  (TAC  1) .  A  single  channel  analyzer  (SCA  1)  window 
placed  around  the  group  of  pulses  corresponding  to  the 
neutron  time  peak  of  the  n-y  time  spectrum  provided  one 
of  the  coincidence  signals. 

Neutron  and  y-ray  shape  discrimination 

The  different  rise  times  of  the  neutron  and  y-ray 
signals  from  the  liquid  scintillator  provided  another 
method  of  discriminating  between  the  neutrons  and  the  y 
rays.  From  the  slow  output  of  the  liquid  scintillator, 
the  signals  were  sent  through  a  pre-amplifier  and  ampli¬ 
fier  and  fast  zero  cross-over  discriminator  (FZCD)  to 
obtain  starting  signals  for  the  second  (TAC  2) .  The 
same  stopping  signals  derived  for  the  (TAC  1)  were  used 
to  stop  (TAC  2).  The  output  of  (TAC  2)  consisted  of  two 
pulse  height  groups;  one  corresponding  to  the  neutrons 
and  the  other  to  the  y  rays.  The  second  coincidence 
signal  for  (COIN  1)  was  obtained  by  placing  a  (SCA  2) 
window  around  the  peak  corresponding  to  the  neutrons. 
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Gamma-gamma  coincidence 

To  record  simultaneously  known  y-ray  lines  from 

the  sources  and  to  use  them  as  internal  calibration 

lines,  a  y-y  coincidence  arrangement  was  employed  using 

a  Nal(Tl)  crystal  as  a  second  y-ray  detector.  The  two 

88  60 

y-ray  sources  Y  and  Co,  placed  between  the  two  y-ray 
detectors  provided  calibration  lines  at  898,  1836,  1173, 
and  1332  keV.  The  Nal(Tl)  detector  was  mounted  with  the 
Ge(Li)  detector,  8  cm  directly  above  it.  A  third  (TAC  3) 
was  used  to  generate  a  coincidence  time  peak  between  the 
signals  of  these  two  detectors  and  it  was  triggered  by 
the  same  starting  pulses  that  triggered  (TAC  1) .  The 
timing  signals  extracted  from  the  Nal  crystal  through  a 
pre-amplifier,  amplifier  and  timing  single  channel  analy¬ 
zer  (TSCA )  system  provided  the  stopping  pulses.  Once 
again  the  coincidence  signals  were  derived  by  placing 
a  (SCA  3)  window  around  the  time  peak. 

Coincidence  requirements 

It  was  essential  that  all  coincidence  signals 
arrived  at  the  universal  coincidence  (COIN  1  and  3)  at 
the  same  time.  This  was  accomplished  by  using  the  true 
start  signals  from  (TAC  1)  and  (G&D  3)  to  strobe  all  the 
(TAC'S).  The  linear  signals  from  the  Ge(Li)  detector 
were  self-gated  (TSCA  4)  in  order  to  eliminate  low  energy 
part  of  the  spectrum  and  they  provided  additional  coinci¬ 
dence  signals  for  (COIN  1  and  3) .  Thus  a  triple 
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coincidence  between  the  linear,  n-y  time  and  n-y  shape 
signals  was  required  for  each  n-y  event  accepted  by 
the  (COIN  1) .  The  linear  and  y-y  time  signals  consti¬ 
tuted  the  coincidence  requirement  for  the  y-y  events 
accepted  by  the  (COIN  3) .  The  outputs  of  both  (COIN  1 
and  3)  were  sent  to  another  (COIN  2)  whose  output  was 
used  to  gate  the  analog-to-digital  converter  (ADCA) . 

This  caused  the  (ADCA)  to  accept  the  n-y  and  y-y  events. 
The  output  of  (COIN  3)  was  also  used  for  routing  the  y-y 
events  into  a  separate  region  of  the  computer  memory. 

The  ungated  n-y  time  spectra  were  obtained  from 
the  output  of  the  (TAC  1)  and  accumulated  during  all 
runs  in  a  second  analog-to-digital  converter  (ADCB) . 

Both  (ADC'S)  we re  interfaced  to  a  Honeywell  516  on-line 
computer. 


2 . 3  A-nalysis 


The  centroids  of  all  y-ray  peaks  of  interest  were 

carefully  calculated  and  checked.  Each  spectrum  was  then 

Calibrated  using  the  centroid  positions  of  the  known  y 

rays  from  the  60Co  (1173  and  1332  keV)  and  8dY  (898  and 

1836  keV)  sources.  Since  our  known  y  rays  were  available, 

the  coefficients,  a  ,  of  calibration  were  determined  by 

v 

polynomial  fits  to  the  equation. 


E 

Y 


2 


I 

v  =  0 


v 


a  x 
v 


(2.10) 
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where  E^  is  the  energy  of  the  y  ray  and  x  is  the  channel 
number.  There  were  no  significant  deviations  from  the 
linear  relationship,  E^  =  aQ  +  a-^  x  for  all  spectra. 
Nevertheless,  to  ensure  that  the  calibrations  were  indeed 
linear  even  at  the  high  energy  region  (>1836  keV) ,  the 
channel-separations  of  first,  second  and  full  energy 
peaks  of  all  identified  high  energy  y  rays  were  checked 
and  found  to  be  consistent.  Hence  the  calibrations  were 
extended  to  the  regions  outside  the  calibration  lines. 

Having  determined  the  average  Doppler  shifted 
energies  of  all  the  y  rays  concerned,  the  experimental 
attenuations  were  calculated  according  to  Equation  2.1 

v 

=  E^0  (1  +  F  ( T )  -2-  cos  6)  .  (2.1) 

This  equation  has  the  same  form  as  the  first  order 

Legendre  polynomial  expansion 

1  A 

l  Ak  Pk  (cos  0)  =  Aq  (1  +  cos  6)  (2.11) 


where  A,  are  the  coefficients  and  P,  (cos  0)  are  the 
k  k 

Legendre  polynomials.  Hence  the  average  y-ray  energies, 

E  ,  were  least  square  fitted  to  the  first  order  polynomial 

series  and  the  experimental  attenuations  were  determined 

from  the  ratios  of  the  coefficients.  A, /A  ,  once  the  v 

1  o  o 

was  specified. 


The  mean  initial  recoil  velocity,  v  ,  was  deter¬ 
mined  from  kinematics,  after  the  energy  loss  of  the 


.  i  s  i  3U(  ■  .■ f?nx  o 


J  u,  ,  \  A  v  .  i  IS  •  C  C .  3 
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incident  beam  through  the  target  was  taken  into  account 
(i.e.,  incident  beam  energy  was  taken  at  center  of 
target).  For  the  n-y  coincidence  runs,  the  finite 
angle  of  the  neutron  detector  was  also  taken  into 
account  by  using  the  mean  angle  subtended  by  the  neutron 
detector.  For  the  singles  experiment,  the  mean  initial 

recoil  velocities,  v  ,  were  determined  from  the  formula 

'  o' 

(Wa  63) : 


v  =  v 


+1 

/  a  + 

-i 


e  cos  6  )  o(0  )  d  cos  0 

cm  cm  cm 


cm 


(2.12) 


/  a (0  )  d  cos  0 

cm  cm 


in  which 

v  =  velocity  of  the  center-of-mass  in  the 
cm  1 

laboratory  system 

e  =  ratio  of  velocity  of  the  center-of-mass 
in  the  laboratory  system  to  the  velocity 
of  the  recoiling  nuclei  in  the  center-of- 
mass  system 

0  =  the  angle  of  the  outgoing  particles  (i.e. 

cm 

neutrons)  in  the  center-of-mass  system 

o(0  ) =dif ferential  cross  section  of  the  neutrons 

cm 

in  the  center-of-mass  system. 


The  neutron  differential  cross  sections  were  indepen¬ 
dently  measured  for  levels  up  to  3274  keV  excitation 
using  time-of-f light  technique  (Chapter  4). 


' 


' 
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At  a  bombarding  energy  of  5.70  MeV,  the 
37  37 

Cl(p,n)  Ar  reaction  is  assumed  to  proceed  pre¬ 
dominantly  through  the  compound  nucleus  mechanism. 

The  neutron  differential  cross  section  in  the  center- 
of-mass  system  will  therefore  exhibit  symmetry  about 
90°  (Section  4.1).  This  symmetry  property  applied  to 
Equation  2.12  leads  to  the  result  that  v  =  v  .  In 
other  words,  a  symmetric  angular  distribution  has  the 
same  influence  on  the  average  initial  recoil  velocity 
as  an  isotropic  angular  distribution.  The  average 
initial  recoil  velocities  determined  by  the  two  methods 
(i.e.  by  weighting  with  the  neutron  differential  cross 
section  or  by  a  straight  average)  agreed  to  within  one 
per  cent. 

The  expected  F(x)  vs  x  curves  were  calculated 
using  the  computer  code,  OXDS,  (Ro  69)  which  used  the 
formulae  of  Lindhard  (Li  61  and  Li  63)  to  calculate 
electronic  and  nuclear  stopping.  Large  angle  scattering 
was  also  taken  into  account.  A  typical  plot  of  the 
F(x)  vs  x  curves  is  shown  in  Figure  2.3. 

Since  the  target  was  not  infinitely  thick,  some 
of  the  recoiling  ions  would  escape  and  stop  in  the  Au 
backing.  In  such  cases  the  attenuation  factor,  F(x), 
was  adjusted  for  the  effect  of  the  recoil  stopping  in 
the  Au  backing.  The  target  was  considered  to  be  divided 
into  N  layers  such  that  the  attenuation  factor  for  the 
ith  layer  may  be  written  as 


■■ 


. 


O’t  !  T,  J'  4  •  -  -  I  *  v 
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Figure  2.3  F(t)  vs  t  curves  for  the  2795  keV  level.  0n  is  the  angle  of 

the  outgoing  neutron  and  vQ  is  the  initial  recoil  velocity. 

The  two  curves  represent  the  coincidence  (0  =  170°)  and  the 

singles  (0  =  ave)  cases. 


23 


Fi  (  t)  =  FU(t)  +  F2i(T) 


(2.13) 


where  1  and  2  denote  target  and  backing  respectively 
and  they  are  calculated  according  to  Equation  2.4  with 
the  limits  of  integration  bounded  by  the  thickness  of 
the  layer.  Hence  the  composite  attenuation  factor  is 
simply  the  sum: 


F(t) 


(2.14) 


For  some  levels  where  cascade  feeding  from  higher 
energy  levels  occurs  additional  corrections  to  the 
attenuation  factors  were  made.  These  corrections  took 
into  account  the  direct  population  and  lifetimes  of  the 
levels  involved.  The  method  of  determining  the  correc¬ 
tion  is  presented  in  Appendix  A. 


2 . 4  Results  and  Discussion 
2.4.1  Level  energies 

Energy  level  positions  were  inferred  from  y_ray 
energies  determined  by  the  fitting  procedure  described 
in  the  previous  section.  These  energies  were  determined 
independently  from  the  singles  and  coincidence  data  and 
were  in  good  agreement  with  each  other.  The  results  are 
summarized  in  Table  2.1  together  with  other  available 


determinations . 
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The  presence  of  the  1574  keV  y  ray  in  both  the 
singles  and  coincidence  data  indicates  that  it  origi- 
nated  from  a  state  in  the  Ar.  Its  energy,  however, 
did  not  fit  into  the  decay  scheme  known  at  that  time. 
From  the  yield  curve  study,  its  threshold  corresponded 
closely  to  the  level  at  3.17  MeV  but  a  14  keV  discre¬ 
pancy  would  then  exist  if  it  were  indeed  the  3.17  ->1.611 
MeV  transition.  It  was,  therefore,  postulated  that 
the  3.17  MeV  level  was  a  doublet  of  3171  and  3185  keV 
where  the  former  would  decay  by  the  ground  state  transi¬ 
tion  and  the  latter  would  decay  100%  to  the  1611  keV 
level.  To  obtain  direct  confirmation  of  the  presence  of 
the  proposed  3185  keV  level,  a  neutron  time-of-f light 
spectrum  was  accumulated  at  an  incident  proton  energy  of 
5.90  MeV.  The  experimental  set-up  which  is  described  in 
Chapter  4  was  arranged  to  provide  optimum  conditions  for 
the  separation  of  the  doublet.  A  neutron  group  corres¬ 
ponding  to  the  proposed  level  is  indeed  present  in  the 
time-of-f light  spectrum  (Figure  4.2)  giving  direct  con¬ 
firmation  of  the  3185  keV  level. 

The  3185  keV  level  has  since  been  inferred  by 
other  y-ray  experiments  (Table  2.1)  but  to  date  no  other 
direct  confirmation  by  particle  identification  has  been 
reported. 

The  present  energy  level  determinations  are  in 
agreement  (within  experimental  uncertainties)  with  other 
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results  with  the  exception  of  those  of  Ivascu  et  al 

(Iv  71) .  This  discrepancy  may  be  attributed  to  the  low 

3 

yield  of  the  coincidence  experiment  using  the  (  He,p) 
reaction. 

2.4.2  Lifetimes 

The  liftimes  of  levels  up  to  an  excitation  of 
2795  keV  were  investigated  using  the  coincidence  tech¬ 
nique  and  up  to  an  excitation  of  3274  keV  using  the 
singles  technique.  A  typical  y-ray  spectrum  in  coinci¬ 
dence  with  neutrons  is  shown  in  Figure  2.4.  Gamma-ray 
energy  shifts  corresponding  to  0  between  0°  and  120° 
for  all  the  levels  whose  lifetimes  were  measured  are 
presented  in  Figure  2.5.  The  corrected  attenuation 
factors  and  inferred  lifetimes  for  both  the  coincidence 
and  singles  experiments  are  presented  in  Table  2.2.  A 
summary  of  all  published  lifetimes  are  compared  to  the 
present  results  in  Table  2.3. 

The  1410  keV  level  The  1410  keV  level  was  partially 
populated  (^20%)  by  cascade  decays  from  the  3517  keV 
level  (Iv  71) .  The  attenuation  factors  and  lifetimes 
listed  in  Table  2.2  have  been  corrected  for  this  effect 
and  they  also  include  a  10%  correction  due  to  recoils 
decaying  in  the  target  backing.  The  final  lifetime. 


x=  943_^qq  fs,  is  in  good  agreement  with  the  result  of 
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a)  the  1410  keV  state 


.5  Plots  of  energy  vs  cos  0  for  y-ray  transi¬ 
tions  in  ^Ar  observed  at  Ep  =  5.70  MeV. 
y  rays  detected  in  coincidence  with  neutrons 
(COIN)  and  in  singles  (SINGLES)  are  shown 
for  comparison.  The  dashed  lines  represent 
full  energy  shifts  while  the  solid  lines 
correspond  to  the  least  square  fits  of  the 
experimental  shifts. 
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c)  The  2488  and  2795  keV  states 
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d)  The  3171,  3185  and  3274  keV  states 
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Attenuation  Factors  and  Lifetimes  of  Ar 


Transition 

Attenuation 

Factor  F(t) 

Mean  Lifetimes 

E . 

l 

(keV) 

Ef 

(keV) 

Coin 

Singles 

Coin 

(fs) 

Singles 

(fs) 

1410 

0 

0.17±0.02 

0.1110.02 

945+2°° 

-120 

940+^00 
*  -160 

1611 

0 

<0  .05 

>3000 

2217 

0 

0.2610.03 

0.2010.04 

o  o 

O  00 

1 — 1 

+  1 

o 

CM 

LO 

r  nn+HO 
500-  80 

2488 

0 

0.1910.05 

0.1610.04 

700+460 
'  -260 

640-200 

2795 

0 

0.9110.07 

0.8810.12 

2  2+  ^  ^ 

-  17 

2i+  26 

-  18 

3171 

0 

0.6110.11 

o  o 
co 

+  1 

l n 

00 

3185 

1611 

0.3110.05 

285+  ™ 

—  _)  D 

3274 

0 

0.7510.12 

+  29 

4  7 

-  25 

1611 

0.7210.10 

+  29 
52 
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t  =  1080_j^q  fs  obtained  by  Ragan  et  al  (Ra  71) .  How¬ 
ever,  agreement  with  other  lifetime  determinations 
using  similar  (p,n)  reaction  varies  from  fair  (Lu  72 
and  Kr  73)  to  very  poor  (Ca  71) . 


The  2217  keV  level  The  lifetimes  of  the  2217  level 

obtained  by  the  singles  and  coincidence  methods  show 

+  75 

good  agreement  and  the  average  value,  t  =  510  fs,  is 

—  b  U 

also  in  good  agreement  with  the  results  of  x  =  600±180  fs 

and  x  =  540±145  fs  determined  by  Ivascu  et  al  (Iv  71)  and 

Kruzek  et  al  (Kr  73)  respectively.  Fair  agreement  is 

achieved  in  comparison  with  the  value,  x=  380±50  fs, 

obtained  by  Ragan  et  al  (Ra  71) .  In  the  latter  work,  the 

analysis  of  the  2217  +  g.s.  transition  was  complicated  by 

the  presence  of  a  2232  keV  y  ray  due  to  the  first  excited 
32 

state  in  S  which  may  account  for  the  discrepancy 
between  the  values.  The  lifetime,  x  =910±180  fs,  reported 
by  Luketina  et  al  (Lu  72)  is,  however,  considerably  longer 
than  that  of  all  other  determinations. 


The  2488  keV  level  The  attenuation  factor  for  the  2488 
-+  g.s.  transition  was  modified  to  account  for  recoils 
decaying  in  the  target  backing  and  it  amounted  to  a 
correction  of  about  10%  in  the  lifetime.  The  resulting 
lifetime,  x  =  660^^q  fs  is  in  agreement  with  the  life¬ 
times  of  775+1qA  fs  and  755±170  fs  determined  by  Ragan 
et  al  (Ra  71)  and  Kruzek  et  al  (Kr  73)  respectively. 
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Fair  agreement  is  obtained  in  comparison  with  t  =  400±130 
fs  of  Ivascu  et  al  (Iv  71).  The  other  lifetimes  reported 
by  Caraca  et  al  (Ca  71)  and  Luketina  et  al  (Lu  72)  of 
t  =  88±24  fs  and  t =  1040±190  fs  respectively,  are  in 
very  poor  agreement  with  the  present  determination. 

The  2795  keV  level  The  2795  keV  level  has  the  shortest 
observed  lifetime  in  ^Ar.  Our  value  of  t  =  22^^  fs #• 
extracted  from  the  average  of  the  singles  and  coincidence 
measurements  of  21  and  23  fs  respectively,  has  been 
corrected  for  cascade  feeding  from  the  3274  keV  level. 

The  branching  ratio  presented  by  Ivascu  et  al  (Iv  71) 
was  used.  In  general,  all  reported  lifetimes  are  in  good 
agreement. 

The  3171  keV  level  The  100%  ground  state  transition  of 

the  3171  keV  level  was  weakly  excited  in  the  coincidence 

spectra  and  hence  its  lifetime  was  determined  from  the 

singles  spectra  only.  Separation  of  the  3163  keV  y  rays 
8  5  35 

from  the  Cl(p,p'y)  Cl  reaction  was  necessary  and 
resulted  in  a  slight  increase  in  the  errors  of  the  life- 
time,  t  =  85  ^  fs,  which  is  in  good  agreement  with 
T  =  100±20  fs  and  t  =  81±15  fs  of  Luketina  et  al  (Lu  72) 
and  Kruzek  et  al  (Kr  73)  respectively;  both  of  these 
determinations  have  used  the  (p,n)  reaction.  The  only 
other  reported  lifetime  for  this  level  is  that  by  Ivascu 
et  al  (Iv  71)  of  t  =  30+15  fs  which  is  in  considerable 
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discrepancy  with  the  present  measurement. 


The  3185  keV  level  A  spectrum  of  the  3185-)-  1611  keV 

transition  accumulated  at  the  three  angles  is  shown  in 

Figure  2.6  together  with  the  neighboring  1611 +g.s.  and 

3274  ->1611  keV  transitions.  The  resulting  lifetime  of 
+  70 

t  =  285_j-q  fs  is  in  excellent  agreement  with  t  =  300±60  fs 
and  310±80  fs  reported  by  Luketina  et  al  (Lu  72)  and 
Kruzek  et  al  (Kr  73)  respectively.  The  lifetime,  t  = 
950±450  fs,  inferred  by  Ivascu  et  al  (Iv  71)  is  more  than 
three  times  longer  than  the  present  value  and  this  dis¬ 
agreement  may  again  be  attributed  to  the  bw  experimental 
yield  that  was  mentioned  in  Section  2.4.1. 


The  32  74  keV  level  The  3274  +  g.s.  and  3274  +1611  keV 

+29  +29 

transition  yielded  lifetimes  of  t  =  47_25  an<^  ^2-21 

respectively  for  the  3274  keV  level.  Our  adopted  life- 

+  21 

time,  t  =  50  ^  fs,  is  in  excellent  agreement  with  all 
other  reported  values. 


An  overall  comparison  of  the  several  measurements 

37  . 

of  the  lifetimes  of  Ar  by  DSAM  shows  large  variation 

for  the  levels  below  and  including  the  2488  keV  level. 

An  explanation  of  this  variation  in  the  mean  lifetimes 
may  be  due  to  the  fact  that  for  different  experiments, 
the  recoiling  nuclei  slowed  down  in  different  materials. 
It  has  been  reported  (Br  73)  that  the  lifetimes  deter¬ 
mined  by  DSAM  vary  with  atomic  number  cf  the  slowing  down 
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3  7  3  7 

Figure  2.6  Direct  y-ray  spectra  from  the  'Cl(p,ny)  Ar 

reaction  exhibiting  Doppler-shif ts  of  the  y 
rays  originating  from  the  1611,  3185  and 
3274  keV  states. 


material.  The  lifetimes  of  the  higher  levels  show,  how¬ 
ever,  better  agreement  with  each  other,  in  particular 

3  7  3  7 

those  determinations  made  with  the  Cl(p,ny)  Ar 
reaction.  Except  for  the  1410  keV  level,  the  report 
of  Kruzek  et  al  (Kr  73)  appears  to  be  in  accord  with 
the  present  work. 

3  7 

The  mean  lifetimes  of  Ar  compiled  recently  by 
Endt  and  van  der  Leun  (En  73)  has  been  included  in  the 
last  column  of  Table  2.3  for  comparison.  Our  lifetime 
values  are  in  excellent  agreement  with  those  of  the 
compilation.  One  must,  however,  point  out  that  the 
present  results  were  included  in  obtaining  the  average 
lifetimes  reported  in  the  compilation. 


. 


CHAPTER  3 


GAMMA- RAY  ANGULAR  DISTRIBUTION  OF  37Ar 

The  study  of  y-ray  angular  distributions  from 
aligned  nuclear  states  is  a  useful  tool  for  extract¬ 
ing  information  of  spins  and  decay  modes  of  these 
states.  To  define  nuclear  alignment  one  considers  a 
reaction  of  the  type  A(a,b)B*.  If  neither  the  beam  of 
particles,  a,  nor  the  target  A  is  polarized  then  the 
y- radiation  from  the  residual  state,  B* ,  will  exhibit 
axial  symmetry  about  the  beam  direction.  Furthermore, 
if  the  excited  state,  |J2M2>,  of  B*  has  definite  parity 
then  reflection  about  a  plane  normal  to  the  beam  direc¬ 
tion  would  also  be  symmetric.  Under  these  conditions 
the  state  |J2M2>  is  said  to  be  aligned  i.e.  P(M2)=  P(-M2) 
where  PIM2)  is  the  population  of  the  magnetic  substate 


3 . 1  Condition  of  Nuclear  Alignment 

There  are  several  possible  ways  to  obtain  align¬ 
ment  in  the  nucleus  of  decaying  state,  one  of  which  is 
the  A(p,ny)B*  reaction  near  threshold.  At  low  bombarding 
energy  the  (p,n)  reaction  proceeds  mainly  through  the 
compound  nucleus  (see  Figure  3.1). 
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Compound  Nucleus 


Figure  3.1 


The  compound  nuclear  model  describes  the  reaction 
as  proceeding  in  two  stages;  the  formation  of  the  com¬ 
pound  nuclear  states,  |J^M^>,  by  the  projectile  and  the 
target  nucleus  and  the  eventual  decay  of  |J^M^>  into 
two  reaction  products  (in  this  case  the  neutron  and  the 
residual  nucleus  B*).  If  the  incident  proton  direction 
is  chosen  as  the  axis  of  quantization  then  the  maximum 
spin  projection  along  this  axis  for  any  compound  nuclear 
state  is 

| Mx  |  <;  Jq  +  s1  (3.1) 

where  Jq  is  the  target  spin  and  s^  is  the  spin  of  proton; 
the  projection  of  the  orbital  angular  momentum  of  the 
proton  will  not  contribute  because  the  orbital  angular 
momentum  is  perpendicular  to  the  quantization  axis. 

Now  if  the  reaction  proceeds  near  threshold  the 
outgoing  neutron  will  have  low  orbital  angular  momentum 
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and  hence  one  may  assume  that  the  1^  =  0  partial  waves 
will  predominate.  Consequently,  the  maximum  spin  pro¬ 
jection  of  the  residual  nucleus  B*  is 

I  M2  I  s  Jo  +  SJ  +  sv  •  (3.2) 

Thus  for  an  unpolarized  system,  states  with 
J ^  >  (J  +  1)  in  the  residual  nucleus  B*  will  be  strongly 
aligned  because  Equation  3.2  precludes  the  magnetic 
substates  greater  than  (J  +  1)  from  being  populated. 

A  nuclecir  state  is  said  to  be  strongly  aligned  if  its 
highest  magnetic  substate  ( |M? |  =  J2)  not  populated. 

This  results  from  the  fact  that  the  asymmetry  of  the 
angular  distribution  is  most  sensitive  to  the  lack  of 
population  of  the  highest  magnetic  substate. 

3.1.1  Gamma-ray  decay  of  aligned  state 

The  theoretical  derivation  of  the  angular  dis¬ 
tribution  of  y-ray  decay  from  an  aligned  state  has  been 
treated  in  detail  by  Litherland  and  Ferguson  (Li  61) , 
Poletti  and  Warburton  (Po  65)  and  Rose  and  Brink  (Ro  67). 
The  outline  presented  below  follows  the  notation  and  sign 
conventions  of  Rose  and  Brink. 

Consider  the  y-ray  transition  of  the  state  |J2M?> 
to  the  state  |J,.M3>.  The  probability  amplitude,  AM  ^  (£)  , 
for  emission  of  a  photon  along  k  with  circular  polariza¬ 


tion  q  is  given  by 


. 


,  . 
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am  m  (k)  ="<2?ii)  l 
2  3  LMtt 


qT^J2M2 


TLM>IJ3M3><4R) 


(3.3) 


where  T  is  the  electromagnetic  operator  with  n  =  0  or 
1  for  electric  and  magnetic  interaction  respectively. 

( R )  is  the  rotation  matrix  and  R  =  (aBy)  describes 
the  rotation  that  takes  the  quantization  or  beam  axis  to 
the  direction  of  k. 

When  the  orientation  of  the  spin,  J^,  is  not 
observed  the  transition  probability  for  emission  must 
be  summed  over  ,  namely 


l 

M. 


A«A(k) 


(3.4) 


If  the  radiating  system  is  cylindrically  symmetric,  the 
total  radiative  probability  for  photons  along  k  is 
obtained  by  weighting  the  contribution  from  each  substate 
M2  by  its  population  parameter,  P(M2) 


l 


M 


2 


p(m2) 


A 


m2m3 


(3.5) 


Using  the  normalization  condition 
l  P(M2)  =1  (3.6) 

m2 

the  Wigner-Eckart  theorem 

<J2M2|TLM>|J3M3>  =  (-1)2L(J3LM3m|j2M2)<J2I|t^>||J3> 

(3.7) 


fi.t  fs  j*  .  c> '  U  3 o  , 


. 


43 


and  the  properties  of  rotation  matrices  and  vector 
coupling,  the  expression  for  the  total  radiative  pro¬ 
bability  or  angular  distribution  formula  becomes 


k  q+J-.-J^+L'-L-K 

W(J2+J3fkq)=  I  Bk(J2)  (P  (cos  0)  (-) 

KLL  '  7T7T 


X  (2J2  +  1)  2(LL'q-q|KO)W(J2J2LL' ;KJ3) 


X  ’<J2MT<->||a3><J2||T<r>||J3> 


(3.8) 


with 


BK(J2) 


=  1 


M. 


p (m2)  (-) 


^2  ^2 


(2J2 +  1)  (J2J2M2-M2 |K0)  (3.9) 


VJ2)  =  1 


in  which  ( J-^M^J2M2  |  JM)  is  a  Clebsch-Gordan  coefficient, 
W(JiJiLL,;KJ2)  is  a  Racah  coefficient  and  <  |  |  T^tt>|  |j^  > 
is  a  reduced  matrix  element. 


L  "~2 


If  polarization  is  not  measured  then  one  sums 
over  the  circular  polarization  quantum  number,  q, 


W(J2 ->  J3;k)  =  W(J2 J3;kl)  +  W(J2 J3;k-1)  .  (3.10) 


Equation  3.8  can  be  written  in  a  more  compact  form  if 


one  introduces  the  mixing  ratio 

,  <  TT  > 


<7T> 


<J2 I lTL 


J3>/(2L+  1) 


3, 

"2 


<J2|  |  T<^>  |  |  J3>/(2L+  1) 


(3.11) 
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where-  L,  tt  is  the  lowest  multipolarity  and  the  coeffi¬ 
cient 


Rk (LL5 j2j3) =  (-) 


l+J^-J  +L* -L--K 


{  (2J2+1)  (2L+1)  (2L  1  +1)  }' 


x  (LL'l-l|KO)W(J0J2LL' ;KJ3) 


(3,12) 


Hence  the  formula  of  angular  distribution  for  y  rays  in 
terms  of  the  mixing  ratio  when  circular  polarization  is 
not  observed  is 


W(J2^J3;k)  =  W(0)  ~  .1  .  •tBK(J2)RK(LL,J2J3)PK(cos  6) 


LL  1  TT7T  1  K 


x  [i -  (-i)L+L,«+i',-K]6yyy>}/2 1  .  o.isj 

11  L  Ltt  " 

The  multipolarities  (L,  L*)  must,  however,  satisfy 
the  triangular  relationship, 


J2  —  J31  <  L,L’  <  (J2+J3) 


(3.14) 


In  all  but  a  few  rare  cases,  only  the  two  lowest  multi¬ 
polarities  are  important  as  contributions  (y-ray  tran¬ 
sition  probability  decreases  rapidly  with  increasing 
multipole  order) .  Hence  one  needs  to  consider  only  the 
lowest  two  multipolarities,  L  and  L* ,  whereby  L=  (j^  - 
and  L 1  =  L  +  l.  It  follows  that  for  transitions  between 
definite  parity  states,  the  quantity  [1-  (-)Tj+L  h 714-77  ~ AV] 
in  Equation  3.13  will  vanish  when  K  is  odd.  Moreover, 
one  arrives  at  the  same  result  that  the  odd  K  terms 


"  'y 
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disappear  if  the  decaying  state  | J2M2>  is  aligned.  This 
manifests  itself  from  the  definition  of  the  B„(J0) 
coefficients  (Equation  3.9)  that  they  vanish  for  odd  K 
values.  Therefore  the  y-ray  angular  distribution  from 
an'  aligned  state  when  only  the  lowest  two  multipolari¬ 
ties  are  considered  is  simply: 

W(0)  =  l  Bk(J2)Pk(cos  0) 

K  even 

R„(LLJ0J,,)  +  26  R„  (LL'  J0JO  +  6  2R^  (L  '  L  '  J0  J_  )  ' 
x  - 2_3 - K - 23 - K - 2_3_  (3>15) 

1  +  62  J 

3.1.2  Application  of  the  statistical  compound  nuclear 
model 

In  evaluating  the  right  hand  side  of  the  Equation 
3.15  one  must  first  calculate  the  BK(J2)  defined  earlier 

Jp-M2  i, 

Bk(J2)  =  £  P(M2)(-)  (2J2  +  1)  '  (J2J2M2-M2|  KO)  .  (3.9) 

m2 

The  B  ( J0 )  coefficients  describe  the  nuclear 
K  Z 

alignment  and  depend  on  the  population  parameter, 

P(M2),  of  the  magnetic  substates  of  the  decaying  state 
| J2M2> .  When  the  magnetic  substates  are  equally  popula¬ 
ted  i.e.  P(M2)  =  P  (M2)  for  M2  ^  m'2  ,  the  BK(J2)  coefficients 
will  vanish  indicating  that  the  y-ray  angular  distribu¬ 
tion  from  such  an  oriented  state  will  be  isotropic. 
Consequently,  BK(J2)  may  be  regarded  as  orientation 
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coefficients  describing  the  orientation  of  the  magnetic 
substate  and  hence  the  ensuing  y-ray  angular  distribu¬ 
tion.  The  manner  in  which  the  magnetic  substates  are 
populated,  however,  requires  the  knowledge  of  the 
nuclear  forces  involved  in  the  reaction  mechanism. 

Since  nuclear  forces  are  not  exactly  known  one  must 
rely  on  some  nuclear  model. 

The  statistical  compound  nuclear  (SCN)  model 
assumes  that  the  numerous  compound  nuclear  states 
overlap  to  a  great  extent  and  the  interferences  between 
different  spin  states  cancel  out.  In  essence,  the  for¬ 
mation  and  decay  of  the  compound  nucleus  through  its 
various  channels  may  be  described  by  energy  averaged 
transmission  coefficients.  Consequently,  the  population 
parameters  can  be  determined  from  these  transmission 
coefficients  in  terms  of  the  penetrability  term 


T 


Vi’(E,) 


l  T£"V(E,,) 


£"i" 


(3.16) 


where  T  . (E)  is  an  energy  dependent  transmission  co- 
£  j 

efficient  with  orbital  and  total  angular  momenta  £  and 
j  respectively.  The  unprimed  and  primed  quantities 
represent  the  incoming  and  outgoing  channels  respectively 
and  the  summation  extends  over  all  outgoing  channels. 

In  this  notation,  the  term  T^ (E)  in  Equation  3.16 


4  7 


describes  the  formation  while  the  quantity  in  the 
bracket  that  of  the  decay  of  the  compound  nucleus. 

The  transmission  coefficients  can  be  obtained  by  a 
Hauser-Feshbach  calculation. 

In  terms  of  the  penetrability,  t,  the  y-ray 
angular  distribution  formula  has  been  derived  by 
Sheldon  and  Van  Patter  (Sh  66) .  For  the  case  of  y-ray 
decay  from  an  aligned  state  and  considering  only  the 
lowest  two  multipoles  the  formula  reduces  to. 


W  ( 0 )  = 


(2Sj,  +  1)  (2Jq  +  1)  4 


l  <-) 


o  3  J2  2 


3i  3ok 
J1J2  even 


(2i1  +  1)  (2j1  +  1)  (2J1+  1)3s(2J2  +  l)k 


(£10£10 |kO)W(L1L1j1j1,Ks1)W(J1J1j1j1;KJo) 
W(J1J1J2J2;Kj2)  t  PR  (cos  0) 

'r„  (LLJ-J.J  +  26  R„  (LL'  J0J0)  +  62RT,  (L  '  L  '  J_  J0  ) 

i\  Z  J  J\  Z  J  K.  Z  J 


with  summation  over  all  incoming  and  decaying  channels 
and  even  K.  The  actual  calculation  of  Equation  3.17 
is  carried  out  by  the  computer  program,  'AK'  (Gr  71). 

In  order  to  apply  the  statistical  analysis  to  the 
experimental  results  one  must  ensure  that  a  large  number 
of  compound  nuclear  states  are  excited.  This  is  usually 
achieved  by  using  a  very  thick  target  or  a  broad  incident 


. 
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beam  (i.e.  beam  with  poor  energy  resolution).  The  thick 
target  method  has  been  chosen  for  this  investigation. 

The  experimental  y-ray  angular  distribution  can 
be  expanded  in  a  Legendre  polynomial  series, 

W(e)expt  =  ^  AKQKPK(cos  6)  (3.18) 

^  K  even 

in  which  A„  is  the  expansion  coefficient,  P  (cos  0)  the 
•ft  k 

Legendre  polynomial  and  the  attenuation  factor  des- 

J\ 

cribing  the  finite  size  correction  for  the  y-ray 
detector. 

Thus  by  fitting  the  theoretical  and  experimental 
angular  distributions  (Equations  3.17  and  3.18  respec¬ 
tively)  by  a  least  square  method  (Section  3.3.1)  the 
spin  of  the  decaying  state,  ,  may  be  determined  pro¬ 
vided  the  spin  of  the  final  state,  ,  is  known  or  vice 
versa. 

3 . 2  Experimental  Arrangement 

Angular  distributions  of  the  y  rays  following  the 
3  7  3  7 

Cl (p , ny )  Ar  reaction  were  measured  at  several  proton 
energies  between  3.25  and  5.15  MeV.  The  energies  were 
chosen  to  be  just  above  the  threshold  of  the  levels 
being  investigated  in  order  to  achieve  the  aforementioned 
alignment  of  the  decaying  states.  The  y  rays  were 

3 

detected  by  a  33  cm  Ge(Li)  crystal  mounted  on  a  rotating 


J  .  t  •  •  -if.  J 


I 


T.z  -  t  r;  j . :  ..  j 


trolley  at  a  distance  of  20  cm  from  the  target.  Inten¬ 
sities  were  measured  at  0°,  30°,  45° ,  60°  and  90°  to 

the  beam  axis.  Each  measurement  was  repeated  at  least 

3 

once.  A  43  cm  Ge(Li)  crystal  was  used  in  the  angular 
distribution  measurement  at  proton  energy  of  5.15  MeV. 
Both  detectors  provided  resolutions  of  typically  3.5 
keV  FWHM  for  the  1.332  MeV  ^Co  full  energy  peak. 

The  signals  extracted  from  the  Ge(Li)  detector 
were  fed  into  a  pre-amplifier  and  amplifier  system  and 
then  sent  to  a  single  channel  analyzer  (SCA) .  A  window 
was  set  on  the  SCA  to  select  the  energy  region  of  the 
spectrum  of  interest.  The  output  of  the  SCA  was  used 
to  gate  the  analog-to-digital  converter  which  was  inter¬ 
faced  to  a  Honeywell  516  computer. 

The  PbC^  target  described  in  the  lifetime  study 
(Section  2.2)  was  used  as  well  as  a  thick  target.  The 
thick  target  was  made  by  mixing  a  small  amount  of 
enriched  37C1  (^85%)  in  the  form  of  KC1  into  a  slurry 

of  polyurethane  and  benzene  and  depositing  the  mixture 
onto  a  0.5  mm  thick  tantalum  backing.  The  target  was 
placed  on  a  1.3  mm  thick  aluminum  plate  mounted  at  30° 
to  the  beam  direction  (Figure  3.2). 

Throughout  the  experiment,  the  target  was  cooled 
by  directing  a  jet  of  air  at  the  outside  surface  of  the 
aluminum  plate.  The  beam  was  well  collimated  to  pro¬ 
duce  beam  spots  of  less  than  2  mm  in  diameter  on  the 
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Figure  3.2  Target  arrangement  for  y-ray  angular  distribution 


targets.  With  this  arrangement,  the  counting  rate  in 
the  detector  was  maintained  around  8000  cps  with  beam 
currents  of  0.2  and  1.0  yA  for  the  thick  and  evaporated 
targets  respectively.  No  target  deterioration  was 
observed. 

An  initial  run  at  bombarding  energy  of  3.25  MeV 

3 

was  made  with  a  6  cm  Ge(Li)  crystal  fixed  at  90°  as 

monitor.  The  angular  distribution  of  the  1410  keV  y 

3  7 

ray  from  the  first  excited  state  of  Ar  was  found  to 
be  isotropic  as  expected  for  the  y-ray  transmission 
from  a  J  =  1/2  state  (Ro  65a)  (Figure  3.3).  Hence  in 
subsequent  runs,  the  1410  keV  peak  was  used  as  an 
internal  monitor  and  all  angular  distributions  were 
normalized  to  it.  The  existence  of  an  internal  monitor 
also  eliminated  correction  of  dead  time  and  anisotropie 
arising  from  the  correlation  table. 

A  yield  function  was  also  obtained  with  the 
evaporated  target.  The  beam  energy  was  varied  from 
3.00-5.40  MeV  with  the  Ge(Li)  detector  fixed  at  55°. 
This  angle  was  chosen  since  the  second  order  Legendre 
polynomial  vanishes  (P2 (0=  55°)  =0)  and  the  A4  coeffi¬ 
cient  of  the  distribution,  W(0)  =  Aq  +  A2P2  +  A4P4' 
is  usually  small.  Thus  the  yield  intensities  measured 
would  be  least  affected  by  the  angular  distributions. 
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Figure  3.3  The  angular  distribution  of  the  1410 

keV  transition  to  the  ground  state. 
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3 . 3  Analysis 

Gamma-ray  intensities  were  obtained  by  summing 
over  the  full  energy  peaks  after  the  background  had 
been  subtracted.  The  summing  process  was  done  directly 
from  the  computer  display  and  was  checked  later  by  a 
computer  code  'SUMMIT'  (Gr  71a).  In  cases  where  un¬ 
resolved  peaks  occurred,  the  sums  were  extracted  by 
fitting  them  to  Gaussian  shapes.  The  angular  distri¬ 
butions  were  then  normalized  to  the  1410  keV  y  ray 
before  being  fitted  to  an  even-order  Legendre  polynomial 
series  and  compared  to  predictions  of  the  statistical 
model  of  Sheldon  (Sh  66)  for  the  various  possible  spin 
and  mixing  ratio  combinations. 

The  statistical  model  calculations  were  carried 
out  by  the  computer  code  'AK'  (Gr  71).  All  kinemati¬ 
cally  allowed  neutron  and  proton  exit  channels  as  well 
as  the  four  strongest  alpha  exit  channels  were  included 
in  the  calculation.  The  necessary  transmission  coeffi¬ 
cients  for  these  channels  were  obtained  from  an  optical 
model  calculation  using  the  computer  program  'HAUSER' 

(Da  69) .  The  optical  parameters  used  were  those  of 
Rosen  (Ro  65)  for  protons,  Perey  and  Buck  (Pe  62)  for 
neutrons  and  Davison  (Da  69)  for  alphas.  A  second  set 
of  transmission  coefficients  was  generated  using  the 
code  'ABACUS'  (Au  64)  with  the  same  optical  parameters. 
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Although  the  absolute  values  of  the  two  sets  of  trans¬ 
mission  coefficients  differed  by  as  much  as  a  factor 
of  two,  the  ratios,  T^^/T^^,  differed  by  less  than 
30%  and  the  resulting  statistical  calculations  yielded 
A2/Aq  coefficients  that  differ  by  10%.  Since  the 
experimental  A9/A  coefficients  (Table  3.1)  have  uncer- 
tainties  of  this  magnitude,  the  SCN  analysis  using 
either  set  of  transmission  coefficients  would  give 
essentially  the  same  result.  This  is  expected.  Pilt 
(Pi  69)  has  shown  that  the  usefulness  of  the  statistical 
calculation  results  from  its  insensitivity  to  the  trans¬ 
mission  coefficients. 


3.3.1  Chi-square  analysis 
o 

A  X  subroutine  contained  in  the  'AK'  program 
provided  a  test  of  goodness  of  fit  between  the  experi¬ 
mental  and  predicted  distributions  for  various  spin  and 

2 

multipole  combinations.  The  quantity,  x  /  is  defined  as 


x2  =  ^-tt  l 

n  1  i=l 


n  [W(6i)  -  Y  (  0  i  )  ] 


(3.19) 


[a(0i)] 


where  W(0i)  and  Y^)  are  respectively  the  experimental 
and  predicted  y-ray  yields  at  angle  0^.  a(0^)  is  the 

statistical  uncertainty  in  W(0i)  and  n  is  the  number  of 
degrees  of  freedom. 
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Normally,  a  0.1%  confidence  limit  is  set  on  the 

2 

X  implying  that  all  spin  assignments  that  fall  outside 
the  limit  may  be  rejected  with  only  one  chance  in  a 
thousand  of  rejecting  the  correct  spin  assignment. 

Thus  for  some  angular  distributions  unique  spin  assign¬ 
ments  may  be  made. 

3.3.2  Branching  ratios 

The  y-ray  branching  ratios  were  determined  from 
the  angular  distributions  using  the  zero-order  term  of 
the  Legendre  polynomial  expansion  and  previously 
measured  efficiencies  of  the  Ge(Li)  detectors  (Ch  70). 
They  are  presented  in  Figure  3.4  together  with  the  spin 
and  parity  assignments  extracted  from  this  experiment. 

3.4  Results  and  Discussion 

Angular  distributions  of  y-ray  decays  for  all 

37 

levels  up  to  3274  keV  excitation  in  Ar  were  accumu¬ 
lated  at  incident  energies  of  3.25,  3.40,  4.05  and  4.30 
MeV  with  the  KC1  target  and  at  3.35,  4.55  and  5.15  MeV 
with  the  PbCl2  target  which  yielded  a  much  cleaner 
spectrum.  A  typical  y-ray  spectrum  observed  at  5.15 
MeV  proton  energy  is  shown  in  Figure  3.5.  In  fact, 
contributions  of  y  rays  resulting  from  proton  induced 

reactions  on  potassium  made  it  very  difficult  to  extract 

^  37 

the  yields  of  the  higher  energy  y  rays  from  Ar. 
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Figure  3 . 4 


Energy  levels  and  decay  scheme  of  Ar . 


-  1 191 
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3  7 

Several  targets  of  natural  and  enriched  Cl  were 
tried  before  the  PbC^  target  was  chosen. 

Legendre  polynomial  expansion  coefficients 
for  all  measured  angular  distributions  are  listed  in 
Table  3.1  and  mixing  ratio  determinations  are  summarized 
in  Tables  3.2  and  3.3.  The  errors  in  the  mixing  ratios 
were  calculated  using  the  method  described  by  Cline  and 
Lesser  (Cl  70)  in  which  the  errors  were  assumed  to 
follow  a  statistical  F  distribution.  A  10%  confidence 
limit  of  the  F  distribution  was  used  in  determining  the 
errors . 

The  electromagnetic  transition  strengths  of  all 
possible  transitions  were  calculated  using  the  average 
mixing  ratios  and  the  lifetimes  determined  in  Chapter  2. 
Some  mixing  ratio,  spin  and  parity  combinations  were 
ruled  out  based  on  the  upper  limits  for  y-ray  transition 
strengths  recommended  by  Endt  and  van  der  Leun  (En  74) 
listed  in  Table  3.4.  These  will  be  discussed  below  for 
each  level. 

3.4.1  The  1410  keV  level 

The  angular  distribution  of  the  y-ray  transition 

from  the  1410  keV  level  whose  spin  and  parity  were  known 

(J71  =  l/2+)  was  isotropic  to  within  3%  at  E  =  3.25  MeV 

P 

(Figure  3.3).  The  yields  were  corrected  for  absorption 
in  the  target  backing  and  retaining  plate  and  normalized 
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Table  3.1 

Legendre  Polynomial  Coefficients  for  Gamma-ray  Angular 

3  7 

Distributions  of  Ar 


Transition 

(keV) 

Bombarding  energy 
(MeV) 

a9/a 

l  o 

A. /A 

4'  o 

1410  ->  g .  s  . 

3.25 

-0.01±0.03 

0.0010.03 

1611  ->  g .  s  . 

3.35 

0 . 24  ±0 . 02 

-0.1010.02 

3.40 

0 . 18±0 . 03 

-0.0710.03 

4.05 

0.1810.01 

-0.0310.01 

4.30 

0.1610.01 

-0.0310.01 

4.55 

0.1710.01 

-0.0810.02 

2217  ->  g .  s  . 

4.05 

0.2910.02 

-0.0110.02 

4.30 

0.23+0.01 

-0.0H0.01 

- 

4.55 

0.2210.01 

-0.0410.01 

2488  ->  g .  s  . 

4.55 

0.0510.01 

0.0010.01 

5.15 

0.0510.01 

0.02+0.01 

2  795  ->  g .  s  . 

4.55 

-0.01+0.03 

0.0710.03 

5.15 

0.0110.01 

0.0210.01 

3171  ->  g.  s . 

5.15 

-0.3610.08 

-0.0410.08 

3185  ->1611 

5.15 

0.4710.01 

0.02+0.02 

32  74  ->  g .  s  . 

5.15 

-0.0810.01 

-0.0H0.01 

re 
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Table  3.2 

Mixing  Ratios  of  the  1611  keV  Level 


E 

P 

(MeV) 

S  C  N  Analysis 

arctan  6  ( ° ) 

Model  Independent  Analysis 

arctan  6  (°) 

3.35 

8  ±  4 

8  ±  3 

3.40 

12  ±  4 

13  ±  5 

4.05 

7  ±  4 

10  ±  3 

4.30 

8  ±  3 

8  ±  3 

4.55 

8  ±  4 

10  ±  3 

. 


Table  3.3 

3  7  2 

Summary  of  Mixing  Ratios  in  Ar  for  x  <0.1%  Confidence 
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Table  3.4 


Recommended  Upper  Limits  for  Gamma-ray  Transition  Strengths 


Transition 


Upper  Limits  (W.u.) 


El  0.03 

E2  100 


E3 


100 


Ml 


10 


M2 

M3 


3 

15 


2) 


1)  Reference  En  74 


2)  Reference  Ne  74 


\o  •  ot ; . .  }  ».* 


to  the  1410  keV  peak  from  the  monitor  detector. 


3.4.2  The  1611  keV  level 

The  y-ray  angular  distributions  of  the  1611  keV 
level  to  the  ground  state  transition  were  measured  at 
Ep  =  3.35,  3.40,  4.05,  4.30  and  4.55  MeV.  The  SCN 
model  analysis  resulted  in  a  unique  spin  assignment  of 
7/2  (see  Figure  3.6  of  the  angular  distribution  at  E^  = 
3.35  MeV  and  its  fit  to  the  SCN  model  prediction) 

which  is  in  agreement  with  previous  assignment  of  J  = 

7/2  (Ra  71  and  Ta  71) . 

The  y-ray  yield  curve  recorded  from  3.00  to 
5.40  MeV,  however,  shows  fluctuations  (Figure  3.7) 
indicating  that  some  non-statistical  contributions  were 
present.  Since  the  application  of  the  SCN  model  in  the 
analysis  depends  on  the  statistical  averaging  of  all 
partial  waves,  these  fluctuations  reflect  enhanced  con¬ 
tributions  from  some  dominant  partial  waves.  This  may 
affect  the  population  of  the  magnetic  substates  of  the 
decaying  level  and  can  result  in  incorrect  mixing  ratio 
and  spin  assignments.  Consequently,  in  addition  to  the 
SCN  analysis,  all  the  1611  keV  angular  distributions  were 
subjected  to  a  model  independent  study. 


Model  independent  calculation 

In  the  SCN  calculation  the  population  parameters, 
P(M2),  were  determined  from  the  penetration  factor,  t. 


V  .  •  •'  ;  ns  ,  . 


64 


*o  — 

C4 


O 

O- 


o 

CO 


O 

r> 

i 


o 

o 

I 


X 


CO 


c 

o 

•- 

u 

o 


S1NOOD 


Figure  3.6  The  angular  distribution  and  x  plot  for  the  1611  keV  transition 

to  the  ground  state. 
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Figure  3.7  Excitation  function  of  the  1611  keV  y  rays 

observed  at  0  =  90°  for  proton  beam  energies 
from  3.30  to  5.40  MeV.  Arrows  indicate 
energies  at  which  the  angular  distributions 
were  measured. 
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The  model  independent  calculation,  however, allows  the 
population  parameters  to  vary  freely  in  search  for 
best  fit  with  experimental  data.  For  each  set  of  P (M^ ) , 
the  coefficients  were  calculated  and  with  the 

BK^2^'  *'*le  an9ular  distributions  were  obtained  for 
different  mixing  ratios.  The  experimental  angular  dis¬ 
tributions  were  then  least  square  fitted  to  the  calcu- 

2 

lated  ones  using  the  x  analysis  described  in  Section 

3.3.1.  A  computer  code  'MAGPOP1  (Wo  72a  and  Appendix  B) 

was  written  to  facilitate  the  search  procedure.  The 

phase  convention  of  Rose  and  Brink  (Ro  67)  for  the 

mixing  ratio  has  been  adopted  throughout  the  analysis. 

Some  Bt.(J0)  coefficients  determined  were  not 

unique  combinations  of  population  parameters;  rather 

several  sets  of  population  parameters  yield  similar 

Bt.(J0)  values.  Nevertheless,  for  each  set  of  angular 
K  A 

2 

distributions  only  two  mixing  ratio  values  have  x  that 

fall  below  the  0.1%  confidence  limit.  These  mixing  ratios 

are  presented  together  with  the  ones  determined  by  the 

SCN  model  calculations  in  Table  3.2.  The  mixing  ratios 

are  in  good  agreement  for  all  beam  energies. 

If  one  adopts  the  lifetime  of  6.3±0.2  ns  for  the 

1611  keV  level  (En  73),  the  spin  and  parity  assignment 

of  7/2+  can  be  ruled  out  as  it  leads  to  a  M3  transition 

+  1  9 

strength  >45  W.u.  The  larger  mixing  ratio,  6  =  4.0_1*1, 
for  the  7/2”  assignment  is  also  ruled  out  as  it  would 
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require  an  E3  transition  strength >  107  W.u..  The 
accepted  value  of  the  mixing  ratio,  6  =  0.16 ±0.04, 
is  in  good  agreement  with  the  results  of  Champlin  et 
al  (Ch  71)  of  6  =  0.15±0.04,  Taras  et  al  (Ta  71)  of 
6  =  0.22±0.11  and  Ragan  et  al  (Ra  71)  of  6  =  0.09±0.03. 

The  consistency  of  the  mixing  ratio  determinations 
by  the  SCN  and  model  independent  methods  and  their  agree¬ 
ment  with  previous  results  indicate  that  the  yield 
fluctuations  do  not  cause  any  appreciable  error  in  the 
SCN  model  assignments.  Therefore  the  SCN  model  has 
been  used  for  the  analysis  of  the  other  levels. 

3.4.3  The  2217  keV  level 

The  2217  keV  level  decays  to  the  ground  state  only 
and  its  angular  distributions  yield  a  spin  assignment  of 
J  =  5/2  or  7/2  (Figure  3.8)  which  is  in  agreement  with 
the  results  of  Naude  et  al  (Na  68)  and  Ragan  et  al 
(Ra  71)  . 

Negative  parities  may  be  ruled  out  for  both  spins 
(5/2  and  7/2)  because  their  transition  strengths  exceed 
the  upper  limits  quoted  in  Table  3.4.  In  the  case  where 
J  =  5/2,  negative  parity  would  lead  to  a  M2  transition 
strength >  51  W.u.  and  where  J  =  7/2,  the  smaller  mixing 
ratio,  6  =  -0.03±0.04,  would  yield  a  M2  transition 
strength  >117  W.u.,  while  the  larger  mixing  ratio, 

6  >  s.l*1^  would  yield  an  E3  strength  >  1.1  x  105  W.u.. 
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Figure  3.8  The  angular  distribution  and  x  plot  for  the  2217  keV  transition 

to  the  ground  state. 
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To  be  consistent  with  measured  lifetime,  the 
spin  and  parity  for  this  level  is  limited  to  J  =  5/2+ 
or  7/2+  which  is  in  agreement  with  previous  assignments 
of  J  =  5/2  or  7/2  by  Naude  et  al  (Na  68)  and  Champlin 
et  al  (Ch  71)  and  J  =  7/2+  by  Ragan  et  al  (Ra  71) . 

Of  the  two  mixing  ratios  determined  for  the 
transition  ( 7/2+ ->  3/2+)  ,  the  larger  one  is  eliminated 

fl 

on  the  basis  of  a  M3  transition  strength  of  >  4.2  x  10 
W.u..  The  smaller  value  (6  =  -0.03±0.04),  however,  is 
in  good  agreement  with  the  values  of  6  =  -0.03  ±  0.04 
determined  by  Ragan  et  al  and  that  of  6  =  0.05±0.03  by 
Champlin  et  al  (Ch  71) . 

3.4.4  The  2488  and  2795  keV  levels 

The  ground  state  transitions  of  both  the  2488 
and  2795  keV  levels  could  be  fitted  with  either  J=  3/2 
or  5/2  (Figures  3.9  and  3.10).  These  levels  have 

ir  — 

previously  been  assigned  the  values  of  J  =  3/2  and 
5/2+  respectively  (En  73) .  Accepting  these  as  the 
correct  spin  and  parity  assignments,  the  smaller  mixing 
ratio  6  =  -1.6+?*j  for  the  transition  from  the  2488  keV 
level  may  be  ruled  out  as  it  would  lead  to  a  M2  tran¬ 
sition  strength  >24  W.u..  The  average  mixing  ratio  of 
5=  -0.21±0.13  or  >4.7  for  the  transition  from  the  2795 
keV  level  is  in  agreement  with  the  values  of  6  =  0.00±0.15 
or  8 . 0  ±1 . 5  obtained  by  Champlin  et  al  (Ch  71). 
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Figure  3.9  The  angular  distribution  and  x  plot  for  the  2488  keV  transition 

to  the  ground  state. 
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Figure  3.10  The  angular  distribution  and  x  plot  for  the  2795  keV  transition 

to  the  ground  state. 


3.4.5 


The  3171  keV  level 


The  3163  and  3171  keV  y-ray  peaks  from  the 
35  37 

Cl  and  Ar  respectively  were  not  resolved  completely 

.  35 

Since  the  3163  keV  level  of  Cl  has  a  relatively  long 

lifetime,  t  ^50  ps  (En  73)  its  peak  would  not  be  Dopple 
shifted.  A  Gaussian  shape  was  therefore  assumed  and 
fitted  to  it.  The  intensity  of  the  3171  keV  peak  was 
obtained  by  subtracting  the  3163  keV  peak  from  the  com¬ 
posite  peaks  (Figure  3.11).  The  resulting  angular  dis- 

2 

tribution  and  its  associated  x  plot  (Figure  3.12) 

7T 

yielded  a  unique  spin  assignment  of  J  =  5/2.  If  J  = 

5/2  is  assumed  for  this  level,  the  mixing  ratio,  6  = 

+n  8  2 

0.78  0*45/  would  yield  a  M2  transition  strength  >24  W.u 

7T  -f- 

making  J  =  5/2  the  accepted  assignment.  These  spin 
and  parity  and  mixing  ratio  determinations  are  consis¬ 
tent  with  those  reported  by  Champlin  et  al  (Ch  71) . 

3.4.6  The  3185  keV  level 

The  3185  keV  level  does  not  decay  to  the  ground 
state  but  only  to  the  1611  keV  level  whose  spin  and 
parity  assignment  is  J77  =  7/2”,  suggesting  that  it  has 
a  high  spin.  The  angular  distribution  of  the  1574 
(3185+1611)  keV  y  ray  indeed  could  be  fitted  only  for 
J  =  9/2  or  11/2  (Figure  3.13). 

Positive  parity  may  be  ruled  out  as  it  would 
require  M2  transition  strength  >374  and  >940  W.u.  for 
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Figure  3.11  The  full  energy  peaks  of  the  3163  and  3171 

keV  y  rays  of  the  ^r>Cl  and  ^Ar  respectively. 
The  dots  are  experimental  points.  The  solid 
curves  represent  the  Gaussian  fits  of  the 
3163  keV  y  rays  with  linear  backgrounds  which 
were  removed  to  give  intensities  of  the  3171 
keV  y  rays. 
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J  =  9/2  and  11/2  respectively.  The  J  =  11/2  assign¬ 
ment  is  also  ruled  out  as  it  would  require  M3  transi- 

g 

tion  strength  >5.9  x 10  W.u..  Hence  a  unique  assignment 

of  J  =  9/2  has  been  adopted  for  this  level. 

The  mixing  ratio  that  corresponds  to  this 
assignment  is  -0.58±0.09.  Since  no  previous  deter¬ 
mination  of  spin,  parity  and  mixing  ratio  were  reported, 
no  comparison  with  other  work  is  possible. 

3.4.7  The  3274  keV  level 

The  3274  keV  level  decays  to  both  the  ground 
state  and  the  1611  keV  level.  The  angular  distribution 
of  the  ground  state  transition  could  be  fitted  with 
either  J  =  3/2  or  5/2  (Figure  3.14).  The  1663  (3274  + 
1611)  keV  transition,  however,  could  be  fitted  with 
J  =  5/2  or  7/2  (Figure  3.15).  Thus  only  J  =  5/2  assign¬ 
ment  is  compatible  with  both  angular  distribution  results 
and  it  is  also  in  agreement  with  previous  assignment  of 
J  =  5/2  (Iv  71  and  Lu  72).  The  allowed  mixing  ratio 

for  the  ground  state  transition,  6  =  -0.03±0.14  or 
+  5  5 

4.0_^'^,  are  in  reasonable  agreement  with  the  values, 

6  =  0.20±0.25  or  2.7±0.5  determined  by  Champlin  et  al 
(Ch  71) . 

As  mentioned  earlier,  during  the  course  of  this 
experiment  several  investigations  have  been  reported  to 
consolidate  the  electromagnetic  properties  of  the  low- 
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Figure  3.14  The  angular  distribution  and  x  plot  for  the  3274  keV  transition 

to  the  ground  state. 
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Figure  3.15  The  angular  distribution  and  x  plot  for  the  1663  keV  transition 

to  the  1611  keV  state. 


37 

lying  levels  in  Ar.  As  a  result  most  of  the  spins 
and  parities  of  these  states  are  known.  Our  spin  and 
parity  assignments  are  consistent  with  the  existing 
information  as  is  evident  in  the  discussion  above. 

The  present  results  also  contributed  to  the  spin  and 
parity  values  established  in  the  recent  summary  of 
Endt  and  van  der  Leun  (En  73) ;  in  particular  those  of 
the  3171,  3185  and  3274  keV  levels. 
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CHAPTER  4 


NEUTRON  TIME-OF-FLIGHT  STUDIES  OF  THE  3 7C1 (p , n ) 3 ?Ar 

REACTION 

In  the  previous  chapter,  the  statistical  compound 
nuclear  (SCN)  model  was  applied  in  the  y-ray  distribution 
analysis.  The  model  independent  analysis  of  the  1611  keV 
level  has  shown  that  the  SCN  model  predictions  were  con¬ 
sistent.  Nevertheless,  its  application  was  based  upon 

37  37 

the  assumptions  that  the  Cl(p,n)  Ar  reaction  proceeds 
via  the  compound  nuclear  mechanism  and  that  a  large 
number  of  levels  are  excited  in  the  compound  nucleus. 

A  study  of  the  neutrons  following  the  reaction  could 
justify  these  assumptions. 

Neutron  angular  distributions  will  exhibit  symme¬ 
try  about  90°  in  the  center-of-mass  system  if  the  reaction 
proceeds  through  the  compound  nucleus  (En  66) .  Hence  a 
search  for  this  symmetry  will  indicate  the  validity  of 
the  first  assumption. 

The  differential  cross  sections  may  be  calculated 
using  the  statistical  model  (Hauser-Feshbach  (Ha  52  and 
La  58))  to  provide  comparison  with  measured  distributions. 
This  would  offer  a  test  for  the  potentials  and  parameters 
used  in  the  calculation  as  well  as  a  test  for  the  appli¬ 
cability  of  the  statistical  theory. 
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In  addition  to  the  justification  of  the  SCN 

model,  the  neutron  angular  distributions  would  also 

enable  one  to  determine  the  average  recoil  velocity  of  the 

3  7 

residual  nucleus  ~  Ar  knowledge  of'  which  is  essential 

in  the  DSAM  singles  analysis  where  the  outgoing  neutrons 

are  not  observed  (Chapter  1) .  Moreover,  a  neutron 

3  7 

spectrum  would  show  the  levels  in  Ar  that  are  directly 
excited  and  thereby  aid  in  the  identification  of  the  y 
rays . 

4 . 1  Compound  Nucleus  and  Statistical  Theory 

Nuclear  reaction  can  be  classified  into  compound 

nuclear  reactions  or  direct  interactions  according  to  a 

certain  time  scale,  the  nuclear  period.  This  is  defined 

as  the  time  taken  for  the  average  bound  nucleon  to  cross 

-22 

the  nucleus  and  is  typically  in  the  order  of  10  “  s 
(Pr  62) .  Direct  interaction  occurs  instantaneously  and 
has  a  reaction  time  in  a  'nuclear  period'  region.  On 
the  other  hand,  one  pictures  the  formation  of  the  compound 
nucleus  as  a  process  in  which  the  energy  of  the  incident 
projectile  is  randomly  distributed  among  the  nucleons 
in  the  target.  A  particle  within  the  compound  system 
may  gain  sufficient  energy  through  successive  collisions 
to  escape,  hence  the  decay  of  the  compound  nucleus. 

Since  this  process  requires  many  nuclear  periods,  one 
assumes  that  the  compound  nucleus  has  no  recollection 


ft 


. 


of  its  formation  when  it  decays.  Consequently,  the  for¬ 
mation  and  decay  can  be  considered  as  two  independent 
processes . 
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At  low  energies,  the  cross  section  exhibits  narrow 
isolated  resonances  which  correspond  to  the  levels  in  the 
compound  nucleus.  The  widths,  r,  of  these  resonances 
which  vary  from  around  0.03  to  2.0  eV,  equivalent  to 
lifetimes  (x  =  h/T)  of  2.2  *10  ^  to  3.3  x 10  ^  s,  are 
characteristic  of  the  compound  nucleus  according  to  the 
'nuclear  period'  time  scale.  Therefore  the  average  com¬ 
pound  nuclear  reaction  cross  section  over  any  resonance 

7T 

structure  of  spin  and  parity,  J  ,  can  be  written  as 


-Jtt  Jit  ,  . 

a  ,  =  a  (a) 
aa '  comp 


.Jtt 
’a ' 
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-Jtt 

vj  ■■ 
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(4.1) 


T  <IT 

where  a  (a)  is  the  cross  section  for  the  formation  of 
comp 

the  compound  nucleus  through  channel  a  and  the  factor 

T  /jt  "T  or 

[G , /\  G  ,,]  is  the  branching  ratio  of  its  decay  through 
a  it  cx 

a" 

channel  a' .  The  summation  over  a"  extends  over  all  decay¬ 
ing  channels. 

Since  compound  nuclear  reactions  conserve  particle 
flux  and  are  invariant  under  time  reversal,  it  follows 
that  any  reaction  and  its  reverse  reaction  will  satisfy 
the  reciprocity  theorem. 
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in  which 


( 2  J  +  1) 

(21  +  1)  ( 2 i  +  1) 


(4.3) 


is  the  statistical  spin  factor.  J  is  the  spin  of  the 
compound  nucleus.  I  and  i  are  the  spins  of  the  parti¬ 
cles  .in  channel  a  and  k  is  the  wave  number  for  the 

a 

relative  motion  of  these  particles.  The  superscript,  '  , 

denotes  the  corresponding  quantities  for  the  reverse  re- 

action  whose  cross  section  is  described  by  a  .  .  From 

a'  a 

Equations  4.1  and  4.2  one  obtains  the  relation, 


GJlJ 

a 


a 


,J7T 

comp 


(a) 


(4.4) 


As  the  mass  number  and  the  excitation  energy  of 
the  compound  nucleus  increase,  the  number  of  resonances 
excited  increases  and  the  average  spacing,  D,  between  the 
resonances  decreases.  If  the  experimental  condition  is 
such  that  the  spread  of  the  beam  energy,  AE,  encompasses 
many  of  the  resonances  (AE  >>  D) ,  the  average  cross 
section  must  take  into  consideration  transitions  through 
all  these  resonances  of  different  J77  simultaneously. 


namely , 


'  « 


r  -'O  I  :•  3  <-•5.  '  !  I  O.  1098 


At  some  energies  the  average  widths  of  the 
resonances  will  exceed  the  average  spacing  between 
them  (T  >  D) .  The  resulting  overlap  of  resonances 
introduces  interference  terms  in  the  cross  section. 

In  order  to  deal  with  these  interfering  terms,  one 
invokes  the  statistical  theory  which  assumes  that  there 
are  many  resonances  within  the  energy  interval  and  that 
their  widths  and  separation  energies  are  randomly  dis¬ 
tributed.  As  a  result  the  interference  will  average  out 
to  zero  and  the  interference  between  the  levels  of  the 
compound  nucleus  may  be  neglected. 

Also,  it  has  been  shown  that  the  various  incoming 
and  outoing  partial  waves  do  not  interfere  among  them¬ 
selves  (B1  52  and  Pi  69) ;  a  consequence  of  which  is  that 
the  average  cross  section  can  now  be  described  by  energy 
averaged  transmission  coefficients  which  are  related  to 
the  optical  model  phase  shifts,  6  Q,  according  to, 

Oi  X/ 


T^(a)  =  1 


(4.6) 


The  optical  model  views  the  formation  of  the  com¬ 
pound  nucleus  as  an  absorption  of  the  incident  particle 
by  the  target.  Hence  the  compound  nucleus  cross  section 
is  analogous  to  the  optical-model  absorption  cross  sec- 

o  ,  (a)  =  I  (a) 

abs  ^  comp 

Jtt 


tion. 


(4.7) 
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In  terms  of  the  complex  phase  shifts,  the  absorption 
cross  section  is 


u 

g 

r 

2i  6 

'v 

°abs<“>  =  I  *rf  l 

1  - 

^  a£ 

e 

J  7T  k  2  0 

a 

> 

in  which  the  channel  spin  formalism  has  been  adopted,  i.e. 


I  =  intrinsic  spin  of  target 
i  =  intrinsic  spin  of  projectile 

*4* 

£  H  orbital  angular  momentum  of  the  pair 
s  =  I  +  1  =  channel  spin 

"4"  “4*  -4 

J  =  s  +  £  =  total  angular  momentum. 

From  Equations  4.4,  4.6,  4.7  and  4.8  it  follows  that 

=  IT  l  T^a)  .  (4.9) 

s£ 


The  combination  of  Equations  4.1,  4.8  and  4.9 
yields  the  average  compound  nuclear  reaction  cross 
section  in  terms  of  the  transmission  coefficients  (Vo 


68) 


a  ,  =  £  g  l  T0 (a) 

aa '  ,_2  a  Ln  £ 

S  £ 


k~  Jtt 

a 


l  T  (a') 
s'£'  * 

l  T  »(a") 

s"£" 


(4.10) 


The  unprimed  and  primed  quantities  refer  to  the  incoming 
and  outgoing  channels  respectively.  Again  the  sum  in  the 
denominator  extends  over  all  channels  through  which  the 
compound  nucleus  decays. 


r  J  >  il  ri--  >"i  FT  *  i  X 


,qa  jian.ti.ini  S 


.•  %(  my^n^mc  m  •*  ■  I  '  Ir  i-  l ir< 


nioa  I  •nni-fi  * 


J&n:: 


a-  o  >o  :'  r  .  o  m  ’-jo.  ,■  ■  ;\  .  ^ 


srii  rir. iiiw  ipuoirij  amts  rfo  Ic  •  /o  a  >nsix^  to: an  no .'•->& 


-3louzi  bf^o-imco 


In  the  channel  spin  notations,  the  transmission 

coefficient  is  independent  of  J  and  s.  This  will  no 

longer  be  true  if  spin-orbit  coupling  is  present  and 

one  must  replace  each  l  T^(a)  in  Equation  4.10  by 

l  T 0 .  (a)  where 
j£  ^ 

j  =  i  +  £  (4.11) 


and  (a)  is  related  to  the  optical  model  phase  shift, 

60  .  by 
£  ja  2 


T£.(a)  =  1  - 


2i  S  „  . 


£  ja 


(4.12) 


One  can  alternatively  define  transmission  co¬ 
efficient,  TS£j(a)»  to  depend  on  s  and  J  instead  of  j. 
The  choice  of  s  or  j  as  quantum  numbers  dictates  the 
angular  momentum  coupling  scheme  used.  The  two  trans¬ 
mission  coefficients  are  related  by 


T 


s£J 


(a)  =  l  (2  j  +  1)  (2s  +  1)  Vv  ( £1 J i  ;  j s )  T  .  (a) 

j  J 


(4.13) 


in  which  W(£IJi;js)  is  a  Racah  coefficient.  With  this 
choice  of  transmission  coefficient,  one  simply  adds  the 
subscripts,  s  and  J,  to  each  transmission  coefficient  in 
Equation  4.10  to  include  spin  orbit  coupling. 

Based  on  the  approximation  that  no  interference 
occurs  between  all  incoming  and  outgoing  partial  waves, 
the  SCN  model  can  be  extended  to  describe  differential 
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cross  sections.  The  extension,  however,  concerns  only 
the  geometry  and  not  the  nuclear  physics  of  the  re¬ 
action.  It  involves  lengthy  calculation  of  Racah  algebra 
and  has  been  treated  by  many  authors  (B1  52,  Sh  66  and 
Vo  68) .  The  average  differential  cross  section  as 
derived  by  Vogt  (Vo  68)  is 


dQ, 


-±s  I 

4k  Jtt 


1 

(21+  1)  ( 2 i  +  1) 


l 

sis 1  l  ’ 


fTj,  (a)  Tr  (a') 


l  T  „ (a") 
^s"£"a" 


X  Z  UJ£J;sL)Z  U'JA'JJS'L)  (-)S~S'  PT  (cos  6  ) 

JLf 

(4.14) 

where 

Z  (£J£'  J'  ;sL)  =  (i)^  36  30  {  (2£+l)  (2£'+l)  (2J+1)  (2J'+1)  } 

x  W(£J£' J' ;sL) (£0£'0 |L0)  .  (4.15) 


The  lack  of  interference  has  a  direct  consequence 
on  the  shape  of  the  average  angular  distributions;  mak¬ 
ing  them  symmetric  about  90°.  This  property  of  Equation 
4.14  is  manifested  through  the  coefficient  Z(£J£'J';sL) 
which  vanishes  unless  L  is  even. 


4 . 2  Experimental  Method 

37  37 

To  study  the  neutrons  from  the  Cl(p,n)  Ar  re¬ 
action,  the  time-of-f light  technique  was  used.  A  pulsed 
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proton  beam  of  10  ns  duration  with  repetition  rate  of 
1  mHz  was  obtained  from  the  van  de  Graaff  accelerator 
and  compressed  to  ^0.5  ns  by  a  Mobley  magnet  bunching 
system  (Me  67) . 

The  electronics  associated  with  the  neutron 
detection  and  timing  is  shown  in  Figure  4.1.  A  cylin¬ 
drical  capacitive  pickoff  in  front  of  the  target  signals 
the  arrival  of  the  protons.  These  signals  are,  however, 
delayed  and  used  to  stop  the  time-to-amplitude  converter 
(TAC  1) .  Neutrons  are  detected  by  a  9  cm  diameter  by 
1.9  cm  thick  NE213  liquid  scillator  optically  coupled  to 
an  RCA  8575  photo-multiplier.  The  fast  output  of  the 
detector  is  sent  to  a  constant  fraction  pulse  height 
trigger  (CFPHT)  which  provides  the  (TAC  1)  with  the 
starting  pulses.  The  output  of  the  (TAC  1)  is  a  time 
spectrum  and  is  fed  into  an  analog-to-digital  converter 
(ADC)  interfaced  to  the  on-line  computer  (Honeywell  516). 

Gamma  rays  resulting  from  the  reaction  are  also 
detected  by  the  liquid  scintillator.  In  order  to  reduce 
the  y-ray  background,  a  n-y  discriminator  is  set  up. 

This  discrimination  is  based  on  the  different  rise  times 
of  the  neutron  and  y-ray  pulses  and  is  similar  to  the 
arrangement  described  earlier  in  Section  2.2.  Namely, 
the  fast  and  slow  signals  from  the  detector  are  used  as 
starting  and  stopping  pulses  for  a  second  time-to- 
amplitude  converter  (TAC  2) .  The  different  rise  times 
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Figure  4.1  A  block  diagram  of  the  neutron  time-of-f light  electronics. 
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of  the  neutrons  and  y  rays  yield  two  peaks  in  the  output 
of  (TAC  2)  .  A  window  is  placed  on  the  peak  corresponding 
to  y  rays  in  the  single  channel  analyzer  (SCA)  whose 
output  is  used  to  route  the  (ADC)  and  separate  the  neu¬ 
trons  and  y  rays  into  different  regions  of  the  computer 
memory. 

The  objectives  of  the  neutron  study  as  mentioned 

in  the  beginning  of  this  chapter  are  to  identify  the 

37 

excited  states  of  Ar  and  to  measure  the  neutron  angu¬ 
lar  distributions  leading  to  these  states.  Although  the 
experimental  arrangements  are  basically  the  same  for  both 
investigations,  other  conditions  differ  and  will  be  dis¬ 
cussed  separately. 

4.2.1  Neutron  identification 

37 

In  y-ray  studies  a  new  level  in  Ar  at  3185  keV 
was  proposed  (Section  2.4.1).  To  confirm  its  existence, 
good  experimental  resolution  was  required  to  separate 
it  from  the  3171  keV  level  and  the  conditions  were 
arranged  accordingly.  The  liquid  scintillator  detector 
was  placed  at  a  distance  of  6.32  m  from  the  target  at  0°. 

A  5.90  MeV  beam  of  protons  was  used  to  bombard  a  PbCl2 

3  7  2 

(78%  enriched  Cl)  target  of  0.1  mg/cm  thickness 

evaporated  on  0.5  mm  gold  backing.  A  total  charge  of 

50,000  yC  was  collected  on  the  target  over  a  period  of 

10  h  to  yield  the  spectrum  shown  in  Figure  4.2. 
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Figure  4.2  Neutron  time-of -flight  spectrum  of  the  37C1 (p,n) 37Ar  reaction  observed  with 

a  flight  path  of  6.3  m  at  Ep=5.90  MeV.  A  total  charge  of  50,000  yC  was 
collected  on  the  target.  The  neutron  groups  are  indicated  by  their  energies 
as  determined  from  this  experiment. 
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4.2.2  Neutron  angular  distribution 


Angular  distributions  of  neutron  groups  populat- 

.  3  7 

ing  all  levels  up  to  3274  keV  in  the  Ar  nucleus  were 

measured  at  a  proton  energy  of  5.75  MeV.  The  target 

2 

was  similar  to  the  one  mentioned  above  but  of  1.0  mg/cm 
thickness.  Two  NE213  detectors  were  mounted  on  a  cart  5 
apart  at  a  distance  of  4  m  from  the  target.  The  cart 
was  rotated  in  15  degree  intervals  from  0°  to  145°  with 
repeat  runs  taken  at  30  degree  intervals.  A  third 
detector,  the  monitor  counter,  was  fixed  at  2  m  in  the 
opposite  forward  quadrant.  A  total  charge  of  2000  yC 
was  collected  on  the  target  for  each  position  of  the 
neutron  detector.  The  target  was  air  cooled  and  the 
neutron  yield  was  monitored  throughout  the  experiment. 

No  target  deterioration  was  observed. 

4 . 3  Analysis 

The  peak  areas  of  all  neutron  groups  were  ex¬ 
tracted  with  the  aid  of  the  display  system  of  the  on¬ 
line  computer  and  checked  later  with  the  code  'Summit' 
(Gr  71a) .  The  resulting  sums  were  normalized  with 
respect  to  the  monitor  counts  and  also  corrected  for 
(ADC)  deadtimes  and  conversion  from  laboratory  to  center 
of-mass  angles. 
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In  order  to  calculate  the  absolute  cross  section 
the  efficiencies  of  the  neutron  detectors  were  required. 
The  detector  efficiencies  versus  neutron  energies  were 
obtained  by  a  computer  code  'Effi'  (Gr  67)  (Figure  4.3). 
Since  the  efficiency  curves  were  particularly  sensitive 
for  low  energy  neutrons  the  detector  thresholds  had  been 
determined  carefully  before  and  after  the  runs.  Never¬ 
theless,  for  neutrons  below  1.4  MeV  in  energy  the  error 
in  efficiency  was  determined  by  assuming  a  50  keV  uncer¬ 
tainty  in  neutron  energy  and  for  the  higher  energy 
neutrons,  a  5%  efficiency  error  was  imposed. 

The  absolute  differential  cross  section  was  then 
calculated  according  to  the  formula: 


dq  _  A  N 

dft  ~  N  n  Ax  ft 


(4.16) 


in  which  AN  =  number  of  detector  neutrons 

N  E  number  of  incident  protons 
n  E  number  of  target  nuclei  per  unit  volume 
Ax  E  thickness  of  target 
ft  E  solid  angle  subtended  by  detector. 


Experimentally,  N  is  determined  from  the  total  charge 
collected  on  the  target  and  AN  is  determined  from  the 
area  under  the  neutron  peak  after  all  corrections  have 


been  taken  into  account. 
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Figure  4.3  A  detector  efficiency  vs  neutron  energy  curve. 


4.3.1  Hauser-Feshbach  calculation 


The  theoretical  differential  cross  sections  were 
determined  by  the  statistical  model  (Hauser-Feshbach) 
calculation  discussed  in  Section  4.1.  A  computer  code 
written  by  Smith  (Sm  65)  and  modified  by  Davison  (Da  69) 
was  used  to  facilitate  the  calculations.  All  energeti¬ 
cally  allowed  proton  induced  reaction  (i.e.  (p,p')(p*n) 

and  (p,a))  were  included  as  exit  channels. 

The  transmission  coefficients  necessary  for  the 
computation  of  the  average  differential  cross  section 
were  obtained  from  an  optical  model  routine  within  the 
code  using  a  potential  in  the  form 


U  <r)  =  Vc(r)  -  Vof(r,ro,ao) 
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iW 


r . 
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m  c 
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rdr 


(f (r , r 


o 


(4.17) 


where  V  (r)  is  the  Coulomb  potential  of  a  uniformly 
c 

1/3 

charged  sphere  of  radius  r  A  /  .  V  ,  W  and  V  are 

o  o  s 

the  depths  of  the  real,  imaginary  and  spin-orbit  poten¬ 
tial  wells  respectively.  m^  is  the  mass  of  the  pion 

and  f(r,r  ,a  )  is  the  Wood-Saxon  form  factor, 
o  o 

f(r,rQ,ao)  =  {1  +  exp  [  (r  -  rQAly/3) /aQ]  }  .  (4.18) 


in  which  r  and  a  are  the  radius  and  diffuseness  para- 
o  o 

meters  respectively.  The  subscripts  o  and  i  denote 
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respectively  the  quantities  contained  in  the  real  and 
imaginary  terms  involved.  One  notes  that  the  deriva¬ 
tive  form  of  the  form  factor  in  Equation  4.17  produces 
surface  absorption. 

The  optical  model  parameters  that  were  used  in 
the  calculations  are  tabulated  in  Table  4.1. 

4 . 4  Results  and  Discussion 

The  neutron  identification  investigation  using 
2 

the  thin  0.1  mg/cm  target  at  a  proton  bombarding  energy 

of  5.90  MeV  shows  neutron  groups  leading  to  the  excited 
37 

states  in  Ar  (Figure  4.2).  It  is  evident  that  a  doublet 
corresponding  to  the  3171  and  3185  keV  levels  exists. 

This  is  a  direct  confirmation  of  the  3185  keV  level  that 
was  postulated  from  the  y-ray  experiments  (Chapters  2 
and  3) . 

Angular  distributions  of  all  neutron  groups  popu- 

37 

lating  states  below  3274  keV  excitation  in  Ar  have  been 

measured  at  an  incident  energy  of  5.75  MeV.  The  thicker 
2 

(1  mg/cm  )  target  and  shorter  flight  distance,  however, 
did  not  allow  the  3171  and  3185  keV  states  to  be  resolved. 
Therefore  their  differential  cross  sections  were  not 
calculated.  The  differential  cross  sections  of  all  the 
other  states  are  presented  in  Figure  4.4  (a  to  g) . 

The  theoretical  (Hauser-Feshbach)  differential 
cross  sections  for  various  spin  values  are  shown  by  the 
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Table  4.1 


Optical  Model 


Proton 
(Ro  65) 


vo 

(MeV) 

53.8 

w 

(MeV) 

7.5 

V 

s 

(MeV) 

5.5 

ro 

(fm) 

1.25 

r  . 

1 

(fm) 

1.25 

a° 

(fm) 

0.65 

a . 

l 

(fm) 

0.70 

Parameters 


Neutron 
(Pe  62) 

Alpha 
(Da  69) 

41.35 

185.0 

3.95 

20.0 

7.2 

0.0 

1.32 

1.50 

1.25 

1.70 

0.62 

0.80 

0.65 

0.60 
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(  q  ui)  u p/jdp 


Figure  4.4  b)  The  1410  keV  state 
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(  q  w)  UP/-°P 


Figure  4.4  c)  The  1611  keV  state 
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Figure  4.4  d)  The  2217  keV  state 
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(  qw)  UP/^P 


Figure  4.4  e)  The  2488  keV  state 
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(  q  uj  )  UP/^P 


Figure  4.4  f)  The  2795  keV  state 
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Figure  4.4  g)  The  3274  keV  state 
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solid  and  dashed  lines;  the  solid  lines  represent 
calculations  with  the  accepted  spins,  i.e.  predicted 
by  y-ray  work  (Chapter  3). 

The  agreement  in  the  shape  of  the  angular  dis¬ 
tributions  for  the  measured  and  calculated  differential 
cross  sections  is  fairly  good.  Symmetry  about  90° ,  a 
consequence  of  compound  nuclear  reaction,  is  also 
evident  in  the  differential  cross  sections.  To  provide 
a  more  quantitative  analysis,  the  experimental  differen¬ 
tial  cross  sections  were  fitted  to  the  Legendre  polynomial 
series.  Table  4.2  shows  the  Legendre  coefficients  and 
the  chi-squares  of  their  fits  for  the  different  neutron 
groups. 

The  odd  order  Legendre  coefficients  should 
approach  zero  if  the  reaction  proceeded  through  the 
compound  nuclear  mechanism.  One  observes  that  all  the 
odd  order  coefficients  do  indeed  approach  zero  (within 
errors)  except  for  the  3rd  and  1st  order  coefficients 
for  the  neutron  angular  distributions  leading  to  the 
1410  and  1611  keV  states  respectively.  These  deviations 
of  the  differential  cross  sections  from  symmetry  about 
90°  lead  one  to  question  the  assumption  that  the 
37C1  (p,n) 37Ar  reaction  at  bombarding  energies  below 
6  MeV  is  predominantly  a  compound  nuclear  reaction. 

Direct  reactions  which  usually  exhibit  forward  peaked 
differential  cross  sections  may  have  contributed  to  the 
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reaction  and  caused  the  asymmetry. 

A  literature  survey  of  the  direct  (p,n)  reactions 
shows,  however,  that  direct  reaction  cross  sections  are 
negligible  at  the  present  bombarding  energies  of  <6  MeV 
and  are  most  unlikely  to  cause  the  asymmetry  in  the 
measured  differential  cross  sections.  The  arguments 
against  direct  reaction  contributions  are: 

a)  Direct  (p,n)  reactions  are  described  as  charge- 
exchange  reactions.  This  occurs  most  readily  in  re¬ 
actions  involving  isobaric  analog  states  (states  with 
identical  spin,  J,r,  and  isospin,  T)  because  the  initial 

and  final  states  of  the  reactions  are  basically  the 

37  37  37 

same.  In  the  Cl(p,n)  Ar  reaction,  the  ~  Cl  target 

IT  + 

has  a  ground  state  spin  and  isospin  of  J  ,  T  =  3/2  ,  3/2 

3  7 

and  hence  its  isobaric  analog  state  in  Ar  must  have 

the  same  (J^T).  The  lowest  T  =  3/2  state  in  ^Ar  is 

at  4.989  MeV  excitation  (Bu  68)  and  it  has  been  con- 

37 

firmed  as  the  analog  to  the  ground  state  of  Cl  (FI  68) . 
Since  this  is  the  lowest  analog  state  and  it  cannot  be 
excited  with  the  present  bombarding  energies,  no  direct 
analog  reaction  can  contribute  to  our  differential  cross 
sections . 

b)  It  is  also  possible  to  get  charge-exchange  by 
knocking  a  neutron  out  with  a  proton.  However,  this  type 
of  direct  reaction  does  not  require  the  initial  and  final 
states  to  have  the  same  (Jit,T).  In  terms  of  the  shell 


t  i 
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model  description,  the  proton  may  enter  any  orbit  after 
knocking  out  a  neutron  from  any  orbit.  In  contrast, 
the  direct  analog  reaction  leaves  the  nucleons  of  the 
residual  nucleus  in  the  same  shell  configuration  as 
that  of  the  target  nucleus.  Hence  one  would  expect 
the  knockout  reaction  to  proceed  with  a  much  smaller 
probability  (smaller  cross  section)  than  the  analog 
reaction. 

Most  direct  (p,n)  reactions  have  been  observed 

at  proton  energies  above  10  MeV  (Ma  69  and  Fi  74). 

With  proton  energies  from  7  to  13.5  MeV,  Anderson  et  al 

2  7  2  7 

(An  69)  found  no  forward  peaking  in  the  Al(p,n)  Si 
reaction  cross  sections  below  10  MeV  incident  energies. 
Hence  even  if  knockout  should  dominate  the  direct  (p,n) 
reaction,  one  would  expect  it  to  contribute  insignifi¬ 
cantly  to  our  differential  cross  sections  (proton  incident 
energies  <6  MeV) . 

c)  Results  of  direct  (p,n)  reactions  leading  to 

analog  states  at  Ep=^18  MeV  yield  an  average  neutron 
differential  cross  section  of  1  mb/sr  at  0°  and  an 
average  do (0° ) /da (120° )  ratio  of  four  (Ma  69).  From 
theoretical  treatment  of  direct  knockout  reaction, 

Atkinson  and  Madsen  (At  68)  have  shown  that  its  differen¬ 
tial  cross  section  formula  has  the  same  form  as  that  of 
the  direct  analog  reaction.  Thus  assuming  the  same  for¬ 
ward  backward  ratio  in  the  differential  cross  section 
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for  the  knockout  reaction,  it  would  be  difficult  to 
account  for  the  asymmetry  in  the  observed  differential 
cross  sections  (because  the  knockout  cross  section  is 
much  smaller  than  that  of  the  analog  reaction  as  men¬ 
tioned  earlier) . 

d)  Finally,  if  direct  reaction  is  indeed  present, 

one  would  expect  it  to  be  present  in  all  the  differen¬ 
tial  cross  sections,  especially  that  of  the  ground 
state.  This  is,  however,  not  the  case  as  the  neutron 
differential  cross  section  leading  to  the  ground  state 
exhibits  symmetry  about  90°  (Figure  4.4a). 

The  absolute  values  of  the  present  experimental 
differential  cross  sections  were,  however,  lower  than 
those  calculated  using  the  Hauser-Feshbach  code  and  a 
normalization  factor  of  1.25  was  required.  This  discre¬ 
pancy  could  have  arisen  from  the  detector  efficiency 
calculations  and/or  the  target  thickness.  In  any  case, 
the  general  accuracy  of  the  Hauser-Feshbach  calculation 
is  not  better  than  30%  because  of  uncertainties  in  the 
value  of  the  transmission  coefficients.  Taking  these 
errors  into  consideration,  the  magnitude  of  the  calcu¬ 
lated  and  measured  differential  cross  sections  are  in 


reasonable  agreement. 
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One  concludes  from  these  results  that  the 
37  37 

Cl(p,n)  Ar  reaction  at  incident  energies  below 
6  MeV  proceeds  mainly  through  the  compound  nucleus 
and  that  the  statistical  model  is  indeed  applicable 
in  this  case. 
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CHAPTER  5 


SUMMARY  OF  37Ar 

The  results  of  the  various  investigations  on  the 
37 

Ar  nucleus  are  summarized  in  this  chapter  and  compared 
with  theoretical  calculations. 

The  neutron  experiment  has  identified  all  the 
low-lying  levels  in  37Ar  up  to  3274  keV  (Section  4.2.1) 
and  their  energies  have  been  determined  from  y  rays 
detected  in  singles  and  in  coincidence  with  neutrons 
(Section  2.4.1).  The  spin  and  parity  assignments  for 
these  levels  were  derived  from  the  y-ray  angular  dis¬ 
tributions  and  transition  strengths  based  on  Weisskopf 
estimates  discussed  in  Section  3.4.  Although  no  unique 
spins  could  be  extracted  from  the  neutron  angular  dis¬ 
tributions,  comparisons  of  these  distributions  with  the 
statistical  model  calculations  were  consistent  with  the 
spins  determined  from  the  y-ray  work.  The  neutron  angu¬ 
lar  distributions  do  show,  however,  that  in  the  case  of 
1410  keV  level,  the  spin,  J  =  1/2  is  the  most  probable 
assignment.  Similarly,  the  spin,  J  =  7/2  is  more  proba¬ 
ble  than  J  =  5/2  for  the  2217  keV  level. 

3  7 

The  electromagnetic  transition  strengths  in  Ar 
below  3274  keV  excitation  have  been  calculated  from  the 
lifetimes  (Section  2.4.2)  and  mixing  ratios  (Section  3.4) 
and  are  tabulated  in  Table  5.1.  Spin  and  parity 
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assignments  are  based  on  the  combined  results  of  the 
present  study  except  for  the  2488  and  2795  keV  levels 
for  which  the  previously  established  values  (En  73) 
have  been  adopted. 

Several  shell  model  calculations  on  the  mass 

A  =  37  nuclei  have  been  published.  The  first  major 

study  was  undertaken  by  Glaudemans  (G1  64).  They  used 
2  8 

Si  as  the  core  and  confined  the  remaining  nucleons  to 
the  2s-jy2  anc^  1^3/2  orhits*  More  recent  calculations 
by  Dieperink  and  Brussaard  (Di  69)  and  Wildenthal  et  al 
(Wi  71)  have  enlarged  the  configuration  space  by  allow¬ 
ing  the  excitation  of  nucleons  from  the  ld^.^  to  the 
2Sjy2  anc*  1^3/2  orhits.  Essentially,  this  established 
■^0  as  the  inert  core  and  expanded  the  vector  space  to 
include  the  entire  2s-ld  shell.  The  difference  between 
the  latter  two  calculations  lies  in  the  Hamiltonian  used 
In  the  shell  model  calculation  of  Wildenthal  et 

al  (Wi  71) ,  the  transition  strengths  between  positive- 

3  7 

parity  state  in  Ar  were  also  determined  and  have  been 

included  in  Table  5.1  for  comparison. 

The  above  shell  model  calculations  yield  positive 

parity  states  only.  To  obtain  negative-parity  states, 

32 

Erne  (Er  66)  performed  a  calculation  using  S  as  the 
core  with  one  nucleon  in  the  If 7^2  subshell  and  the  othe 


extra  core,  nucleons  in  the  Id 


3/2 


subshell . 
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The  results  of  these  calculations  together  with 
the  properties  measured  in  this  investigation  are  listed 
in  Figure  5.1. 


All  the  Shell  model  calculations  have  yielded 

J  =  3/2  for  the  ground  state  of  ;Ar;  corresponding 

to  the  last  (nineteenth  neutron)  unpaired  neutron  being 

in  the  ld^^?  orbit.  The  1410  keV  state  with  J1'  =  l/2+ 

has  been  identified  as  a  s  state  with  the  configura- 

3  6 

tion  (2s-jy9  )  ~  ^  ^^5/2  ^  0  *  Since  Ernd's  calculation 
assumed  the  2s^y2  orbit  to  be  closed,  no  s^y2  state  was 
possible.  The  calculated  energies  for  this  level  vary 
from  1.43  to  1.55  MeV  and  are  in  good  agreement  with  the 
measured  value  of  1410  keV. 

The  unique  assignment  of  J  =  7/2  for  the  1611 
keV  state  (Section  3.4.2)  restricts  comparison  to  the 
2.07  MeV  level  in  Ernd’s  report  (Er  66).  According  to 


the  model  pr  ‘diction  this  level  has  the  configuration 
(ld3/2)  ^  df  y/2 )  7/2  and  contains  m°st  of  the  single  par¬ 
ticle  f j, 9  strength. 


The  J1T  -  1/2'  state  at  2217  keV  excitation  has 
been  reproduced  well  by  all  the  calculations.  On  the 
contrary,  the  J7r  =  3/2~  assignment  for  the  24  88  keV 
level  cannot  be  accounted  by  any  of  the  aforementioned 
shell  model  calculations.  This  is  due  perhaps  to  the 
exclusion  of  the  2p3^2  subshell  in  the  calculations. 
Stripping  reactions  (Ro  65a  and  Ch  71)  have  indicated 
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that  this  state  possesses  most  of  the  2p  ,  single 

J/ z 

particle  strength.  The  3/2  state  at  3.50  MeV  exci¬ 
tation  predicted  by  Erne  (Er  66)  has  the  configuration 


^ld3/2^  2  ^lf 7/2* 7/2 # 

The  2795  keV  level  corresponds  to  the  lowest 
calculated  J77  =  5/2+  state.  With  the  exception  of 
Glaudemans '  (G1  64)  prediction  at  3.98  MeV,  all  other 

calculations  reproduced  this  level  in  the  correct  energy 
region.  The  3171  and  3274  keV  levels  with  J77  =  5/2  + 
and  5/2”  respectively  may  be  compared  to  the  second 

7T  + 

J  =  5/2  state  calculated  by  Dieperink  and  Brussaard 
(Di  69)  and  Wildenthal  et  al  (Wi  71) . 


Erne  (Er  66)  predicted  a  J71  =  9/2  state  at  3.87 
MeV  excitation  which  may  describe  the  observed  3185  keV 
level  whose  spin  and  parity  assignments  were  uniquely 
determined  to  be  9/2  (Section  3.4.6). 

One  observes  from  this  comparison  that  the  agree¬ 
ment  between  experimental  and  theoretical  results 
improves  when  the  calculation  involved  a  larger  vector 
space.  Nevertheless  the  overall  agreement  is  far  from 
good.  With  the  exception  of  Erne's  calculation,  the 
exclusion  of  the  negative  parity  orbits  lf^2  and  ^3/2 
may  account  for  the  discrepancy.  This  is  particularly 
true  since  the  1611  keV  state  has  been  identified  with 
the  'single  particle'  f^2  level.  Moreover,  Wildenthal 
et  al  (Wi  71)  pointed  out  that  exclusion  of  negative- 
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parity  states  does  not  affect  the  calculation  of  positive- 
parity  states  since  negative-  and  positive-parity  confi¬ 
gurations  do  not  mix.  But  the  results  are  affected  by 
the  exclusion  of  positive-parity  states  due  to  excitation 
of  an  even  number  of  particles  from  the  sd  to  the  fp 
shell . 

Finally  one  must  mention  that  these  shell-model 

calculations  were  performed  when  experimental  information 
3  7 

on  Ar  was  lacking.  The  results  of  recent  investiga¬ 
tions  and  this  experiment  established  most  of  the  pro- 

3  7 

perties  of  the  low-lying  levels  in  Ar.  This  will 
provide  a  better  comparison  for  future  theoretical  cal¬ 
culations  . 
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CHAPTER  6 


A  STUDY  OF  THE  STATISTICAL  MODEL 
29  29 

WITH  THE  Si (p , p ' y )  Si  REACTION 

6 . 1  Motivation 

The  previous  chapters  have  shown  that  the  statis¬ 
tical  compound  nuclear  (SCN)  model  was  successful  in  its 

37  3  7 

description  of  the  Cl(p,n)  Ar  reaction.  The  assump¬ 
tion  that  many  levels  in  the  compound  nucleus  were 
excited  was  inferred  from  the  use  of  thick  targets. 
Agreement  of  the  experimental  results  with  model  cal¬ 
culations  suggests  that  probably  sufficient  levels  were 
excited.  Nevertheless,  it  would  be  interesting  to  inves¬ 
tigate  how  many  levels  are  required  to  justify  the 
statistical  approach. 

One  recalls  (Section  4.1)  that  as  the  mass  of  the 
nucleus  and  also  the  excitation  energy  increase,  the 
number  of  compound  nuclear  levels  increase.  In  order 
to  study  the  effect  of  the  number  of  compound  nuclear 
levels  on  the  statistical  theory,  one  must  have  an 
experimental  condition  whereby  one  can  excite  a  few 
or  many  levels  in  the  compound  nucleus  by  varying  the 
target  thickness. 

An  appropriate  mass  region  to  conduct  this  inves¬ 
tigation  is  in  the  mid  s-d  shell  (A  ^30)  where  studies 
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on  structure  and  resonance  have  reported  (Sh  70  and 

Ru  72)  level  densities  suitable  for  the  experimental 

condition  mentioned  above.  In  particular,  the  proton 

29 

capture  and  inelastic  scattering  yields  on  Si  nucleus 

have  exhibited  considerable  fluctuation  suggesting  low 

level  density  (De  70  and  Di  71) .  In  addition  to  its  low 

29 

level  density,  the  Si  nucleus  has  several  desirable 
properties  which  would  make  it  ideal  for  this  investi¬ 
gation. 

29 

Firstly,  the  Si  nucleus  has  a  ground  state  spin 

7T  -f- 

and  parity  (J  =  1/2  )  which  would  enhance  nuclear 

alignment  (Section  3.1.1).  Secondly, at  proton  incident 

energies  below  4.8  MeV,  only  the  (p,p')  reaction  is 

energetically  possible  (all  other  proton  induced  reactions 
29 

on  Si  have  larger  negative  Q-values) .  Lastly,  the  level 

energies  and  electromagnetic  properties  of  the  low-lying 
29 

levels  in  Si  are  well  known  to  afford  comparison. 

Particle-gamma  angular  correlation  studies  of 

Becker  et  al  (Be  67)  confirmed  the  spin  and  y-decay 

29 

modes  of  the  low-lying  levels  in  the  Si  nucleus 
(Figure  6.1).  Lifetimes  of  these  states  have  also  been 
subjected  to  thorough  investigations  by  many  authors 
using  different  reactions  (Ba  68a,  De  70  and  Me  71) . 
Detailed  stripping  (El  72)  and  two  nucleon  transfer 
(Cu  69,  Me  71  and  Na  72)  studies  which  included  useful 
optical  model  parameters  from  DWBA  analysis  have  been 
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reported. 

29 

The  theoretical  aspects  of  the  Si  nucleus  have 
also  received  considerable  attention  because  it  is  in 
the  mass  region  where  nuclear  deformation  changes  from 
prolate  to  oblate  (Br  57) .  The  low-lying  positive- 
parity  states  have  been  described  in  terms  of  the  Nilsson 
model  (Ni  55)  as  members  of  two  mixed  rotational  bands 
with  oblate  deformation  (Br  57  and  Hi  69).  Single  par¬ 
ticle  interpretation  had  been  unsatisfactory  (Ba  68) 
until  the  advent  of  large  model  space  calculations. 

Recent  shell  model  calculations  have  given  a  fair  overall 
description  (Vo  72  and  Wi  73) . 

In  view  of  these  desirable  features,  the 
29  29 

Si(p,p*Y)  Si  reaction  was  chosen  for  this  investi¬ 
gation  to  determine  when  sufficient  averaging  over  com¬ 
pound  states  is  achieved  such  that  the  experimental 
results  can  be  compared  to  statistical  theory  predictions. 

6 . 2  Experiment 

Gamma-ray  angular  distributions  following  the 

n  q  2  9 

Si(p,p*y)  Si  reaction  were  measured  with  a  52  cc  Ge(Li) 
detector  at  0°,  31°,  55°  and  90°.  Each  angle  was  repea¬ 
ted.  A  second  Ge(Li)  detector  of  48  cc  volume  was  fixed 
at  90°  in  the  opposite  forward  quadrant  as  a  monitor. 

The  detector  distance  from  the  target  were  12  and  8  cm 


c  v)  f  ,a !  b 


pc 


or 


VI  w  .■  * .■'“’•Mr  d  .1  °C  y<  DiJ 
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respectively.  Lead  disks  of  3  mm  thickness  were  placed 
in  front  of  both  detectors  to  absorb  the  low  energy  y 
rays. 

2 

The  target  consisted  of  94  yg/cm  of  SiC>2  (61% 

29 

enriched  Si)  evaporated  onto  a  gold  backing  of  0.13 
mm  thickness.  It  was  mounted  at  30°  to  the  beam  direc¬ 
tion  and  air  cooled  in  the  arrangement  described  in 
Section  3.2. 

The  energy  of  the  proton  beam  obtained  from  the 
van  de  Graaff  accelerator  was  varied  from  3.00  to  3.60 
MeV  in  10  keV  intervals.  To  ensure  energy  stability 
and  accuracy,  extreme  care  was  taken  in  increasing  beam 
energy  without  changing  other  conditions.  Throughout 
the  experiment,  the  target  current  was  maintained  at 
^0.7  yA  and  monitored  by  an  audio  alarm  system  which 
would  sound  if  the  target  current  should  fall  below  a 
preset  level.  When  a  run  was  interrupted  the  last  run 
was  repeated  before  resuming  the  experiment.  The  repro¬ 
ducibility  was  within  the  statistical  errors  and  no 
renormalization  was  required. 

Following  the  above  series  of  runs,  the  target 

was  replaced  by  a  thicker  one.  This  target  consisted 

29 

of  90%  enriched  Si02  mixed  in  a  slurry  of  polyurethane 

and  benzene  smeared  onto  a  gold  backing.  Its  estimated 

2 

thickness  was  ^5  mg/cm  .  The  experiment  was  repeated  in 
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100  keV  intervals  with  the  angular  distributions  taken 
at  the  same  angles  mentioned  earlier  and  the  addition 
of  a  point  at  71°. 

The  electronics  and  computer  configuration  asso¬ 
ciated  with  the  detection  and  data  collection  were 
similar  to  that  previously  described  (Section  3.2). 

Intensities  of  the  1273  and  2028  keV  y-ray  tran- 
sitions  to  the  ground  state  in  Si  were  extracted  from 
the  full  energy  peaks  by  the  same  procedure  as  described 
in  the  earlier  chapters.  They  were  then  normalized  with 
respect  to  the  1273  keV  peaks  in  the  monitor  detector 
and  corrected  for  absorption  in  the  target  backing  and 
retaining  plate.  Dead  times  in  the  analog-to-digital 
converters  were  also  compensated  for  although  they  accoun 
ted  for  less  than  one  percent  error. 

6 . 3  Results  and  Discussion 
6.3.1  Thin  target 

Preliminary  survey  runs  with  the  thin  target 
(94  yg/cm  )  had  indicated  that  the  results  disagreed 
with  the  SCN  model  predictions.  This  discrepancy  may 
be  attributed  to  insufficient  statistical  averaging  as 
a  result  of  too  few  levels  being  excited  in  the  compound 
nucleus.  Hence  the  experimental  y-ray  angular  distri¬ 
butions  measured  with  the  thin  target  were  not  subjected 
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to  the  normal  statistical  analysis  (Section  3.3). 

Instead  they  were  least  square  fitted  to  the  even 
Legendre  polynomial  series,  W(0)  =  Aq  +  A2^2P2 +  A4^4P4 
(notations  similar  to  Equation  3.18).  The  coefficients 
extracted  from  the  fits  were  expressed  as  a  function  of 
the  bombarding  energies. 

Figure  6.2  show  such  a  plot  for  the  1273  and  2028 
29 

keV  y  rays  of  Si.  The  zeroth  order  coefficients,  Aq, 
represent  the  yield  function  and  the  higher  order  co¬ 
efficients,  A~/A  and  A„/A  ,  the  distribution  functions. 
Rather  than  complicate  the  graphs  only  typical  error-bars 
are  shown.  The  distribution  coefficients  obtained  from 
the  SCN  calculations  are  shown  by  the  dotted  line  for 
comparison.  These  theoretical  coefficients  are  calculated 
with  the  adopted  spins  and  mixing  ratios  (En  73) .  The 
transmission  coefficients  required  for  the  statistical 
calculations  were  obtained  from  the  code  'Hauser*  (Da  69) 
using  the  optical  model  parameters  from  the  report  of 
El-Naiem  et  al  (El  72).  They  are  listed  in  Table  6.1. 

Before  analyzing  the  results  further  one  questions 
if  the  resonance-like  peaks  in  the  yield  curves  are 
authentic,  i.e.  do  they  arise  from  the  reaction  and  not 
from  fluctuation  in  the  experimental  conditions?  The 

yield  of  the  1273  keV  y  ray  had  been  measured  by  Koen 

29  29 

(Ko  72)  using  the  same  Si(p,p'y)  Si  reaction  with  a 
2  29 

5  yg/cm  enriched  Si  target.  His  result  shows 
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Figure  6.2  Observed  Legendre  polynomial  coefficients  vs  energy  for  the  1273 

and  2028  keV  y  ravs. 
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Table  6.1 
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Optical  Model  Parameters  for  the  Si(p,py)  S 


vo 

(MeV) 

50.31 

w 

(MeV) 

7.95 

V 

s 

(MeV) 

6.2 

ro 

(fm) 

1.25 

r . 

l 

(fm) 

1.25 

ao 

(fm) 

0.64 

a . 

(fm) 

0.74 

Reaction 


1)  Notations  similar  to  those  in  Table  4.1 


, 


12  7 


resonance-like  peaks  that  occur  at  the  same  energies  as 
our  peaks.  This  confirms  that  the  peaks  observed  in 
Figure  6.2  are  due  to  resonance  structures  in  the  reac- 
t  ion „ 

The  fluctuations  of  the  measured  angular  distri¬ 
bution  coefficients  about  the  SCN  line  is  expected  since 
th  e  thin  target  was  chosen  specifically  to  excite  only 
a  few  levels  in  the  compound  nucleus.  Perhaps  the  pro¬ 
perties  of  these  levels  do  influence  the  fluctuations 
in  the  coefficients  and  therefore  an  attempt  was  made  to 
see  if  some  correlation  does  exist. 

The  statistical  model  assumes  that  compound  nuclear 
levels  of  all  spin  and  parity  combinations  are  excited 
provided  they  satisfy  the  conservation  laws.  The  calcu¬ 
lation  then  sums  (neglecting  any  possible  interference) 
over  partial  wave  contributions  leading  to  and  from  all 
these  states.  When  the  compound  nuclear  level  density 
is  low,  only  a  few  levels  are  excited  or  in  the  extreme 
case  perhaps  only  a  single  level  is  excited.  In  these 
cases  a  CN  model  calculation  may  still  provide  a  reason¬ 
able  description  of  the  reaction. 

To  study  the  effect  of  low  level  density  on  the 
SCN  calculation,  one  restricts  the  summation  over  all 
incoming'  and  outgoing  partial  waves  in  Equation  3.17  to 
those  for  a  specific  compound  nuclear  spin  and  parity, 

J^.  This  was  accomplished  by  modifying  the  summation 
routine  in  the  code  *AK’  (Gr  71)  and  resulted  in  a  set 
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of'  theoretical  coefficients  for  distinct  values. 

They  are  presented  graphically  in  Figure  6.3  for  com¬ 
parison  with  Figure  6.2. 

The  range  of  angular  distribution  coefficients 
calculated  for  a  specific  encompasses  all  the  measured 
values.  This  indicates  that  the  CN  calculations  for  spe- 
cific  could  account  for  the  observed  fluctuations  in 
the  angular  distributions.  Comparison  of  the  experimen- 
tal  and  theoretical  coefficients  for  specific  may  also 
yield  spin  and  parity  information  on  possible  resonance 
peaks.  The  distribution  coefficients  of  both  the  1273 
and  2028  keV  y  rays  approach  zero  at  3.23  MeV  suggesting 
transition  through  J-j  =  0  states  in  the  compound  nucleus. 
Therefore  the  resonance  peak  at  this  energy  may  correspond 
to  a  state  of  =  0.  In  practice,  however,  incomplete 
averaging  of  the  interference  terms  could  result  in  dis¬ 
agreement  with  the  simple  case  of  the  single  CN  resonance. 
S imulation  of  Thicker  Targets  by  Yield  Averaging 

So  far  the  experimental  evidence  is  consistent 
with  the  hypothesis  that  excitation  of  insufficient  levels 
in  the  nucleus  causes  the  failure  of  SCN  predictions.  It 

I 

follows  then  that  the  inclusion  of  more  compound  nucleus 
levels  should  overcome  the  discrepancy.  In  order  to 
study  the  effect  of  averaging  over  these  resonances,  one 
combines  the  yields  from  two  adjacent  measurements  of 
successive  energies.  By  repeating  the  process  at  the 
respective  angles  one  simulates  an  angular  distribution 
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Figure  6.3  Legendre  polynomial  coefficients  calculated  with  the  modified 

SCN  model. 
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obtained  from  a  target  twice  the  original  thickness. 

Thus  a  new  set  of  Legendre  coefficients  may  be  ex¬ 
tracted  by  the  least  square  procedure  already  described. 

This  procedure  of  simulating  thicker  targets  is 
repeated  for  larger  and  larger  energy  intervals. 
Following  each  increment  of  target  thickness  (i.e. 
energy  bite) ,  the  least  square  fitting  is  performed. 
Hence  one  generates  a  whole  series  of  Legendre  coeffi¬ 
cients  which  correspond  to  gradual  increases  of  target 
thickness.  A  code  'WAVY'  (Wo  74  and  Appendix  C)  was 
written  to  carry  out  the  computation.  The  results  are 
shown  in  Figure  6.4  (a  to  d)  for  the  energy  intervals  of 
50,  100,  200  and  400  keV  respectively. 

For  comparison,  the  coefficients  determined  with 
the  thick  target  are  presented  in  the  same  manner  in 
Figure  6.5.  Although  some  of  the  experimental  coeffi¬ 
cients  in  Figure  6.5  do  not  fall  (within  experimental 
errors)  on  the  dashed  line,  they  are  accepted  in  the  SCN 
model  analysis  (Section  6.3.2)  because  of  the  uncertain- 

i 

ties  in  the  mixing  ratios.  Since  uncertainties  are 
are  associated  with  mixing  ratio  determinations  one 
must  take  these  into  consideration  when  comparing  with 
the  theoretical  values. 

These  graphs  show  that  as  the  energy  interval  is 
increased,  the  fluctuation  in  the  coefficients  are 
smoothed  out.  When  the  energy  interval  reaches  400  keV 
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Figure  6.4  Angular  distribution  coefficients  obtained  from  averaging  proce 
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Figure  6.4  b)  Energy  interval  of  100  keV 
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Figure  6.4  d)  Energy  interval  of  400  keV 
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Figure  6.5  Angular  distribution  coefficients  obtained  using  the  thick  target. 
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(Figure  6.4d)  the  coefficients  for  both  the  1273  keV 
and  2028  keV  y  rays  vary  smoothly  with  energy  and  they 
resemble  the  curves  for  the  thick  target  (Figure  6.5). 
Also  the  ^2^0  an<*  A4^Ao  coefficients  approach  the  values 
predicted  by  the  SCN  model;  implying  that  the  model  may 
be  successfully  applied  for  this  energy  interval  or 
greater. 

It  remains  then  to  show  that  the  normal  SCN  analy¬ 
sis  would  indeed  predict  the  correct  spins  and  mixing 
ratios  for  the  thick  target  case. 

6.3.2  Thick  target 

The  angular  distribution  of  the  1273  and  2028  keV 

29 

y-ray  transitions  in  Si  measured  with  the  thick  target 
at  different  incident  energies  are  shown  in  Figure  6.6. 
Cascade  feeding  from  the  2028  to  1273  keV  level  contri¬ 
buted  to  less  than  one  percent  in  the  yield  of  the  1273 
keV  y  ray  and  therefore  has  been  neglected. 

Assuming  that  a  sufficient  number  of  compound 
nuclear  states  have  been  excited  with  the  thick  target, 
the  y~ray  angular  distributions  were  compared  to  the 
SCN  model  predictions  in  the  manner  described  in  Section 
3.3.  The  analyses  yielded  unique  spin  assignments  of 
J  =  3/2  and  5/2  for  the  1273  and  2028  keV  levels  respec¬ 
tively  for  all  incident  energies.  These  spins  are  con¬ 
sistent  with  previous  assignments  (En  73).  Sample  chi- 
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Figure  6.6  Angular  distribution  of  the  1273  and  2028  keV 

Y  rays  measured  with  the  thick  target. 
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square  vs  mixing  ratio  plots  for  the  1273  and  2028  keV 
transitions  are  shown  respectively  in  Figures  6.7  and 
6.8  for  the  incident  energy  of  3.40  MeV ;  plots  at  other 
incident  energies  are  similar. 

The  experimental  Legendre  polynomial  coefficients 

for  all  energies  are  presented  in  Table  6.2.  The  mixing 

2 

ratios  corresponding  to  the  minimum  x  obtained  from 
comparison  with  SCN  model  calculations  are  also  tabula¬ 
ted.  The  larger  mixing  ratios  for  both  transitions  are 
not  included  since  they  have  been  ruled  out  in  previous 
investigations  (En  73) . 

The  mixing  ratios  determined  at  all  the  incident 
energies  agree  well 'with  each  other.  The  average  value 
of  6  =-0.13±0.03  for  the  1273  keV  transition  is  to  be 
compared  with  the  accepted  value  of  6  =  -0.206±0.011 
(En  73) .  For  the  2028  keV  transition,  the  average  6  = 
0,07±0.07  is  consistent  with  the  established  fact  that 
the  transition  is  pure  E2. 

6 „ 4  Conclusion 

The  thin  target  study  has  resulted  in  a  visual 
picture  of  the  statistical  process.  By  gradual  simula¬ 
tion  of  thicker  targets,  the  fluctuations  in  the  yield 
and  angular  distribution  coefficients  were  shown  to  be 
daimped  out.  This  is  analogous  to  increasing  the  number 
of  levels  for  statistical  averaging.  When  sufficient 
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2 

X  vs  mixing  ratio  for  the  2028  keV  y  rays 


Figure  6.8 
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Table  6.2 

29 

Angular  Distribution  Coefficients  and  Mixing  Ratios  of  Si 


Transi¬ 

tion 

Spin 

P  roton 
energy 

Legendre  Polynomial 
Coefficients 

Mix¬ 

ing 

ratio 

Min¬ 

imum 

E.  E. 

l  f 

(keV) 

J. 

l 

Jf 

E 

P 

(MeV) 

VA0 

VA0 

Arc- 
tan  6 
(°) 

2 

X 

1273  0 

3/2 

1/2 

3.10 

-0.026 ±0.009 

0.016±0.012 

-12±4 

1.04 

3.20 

-0 . 054  ±0 . 00  7 

0.00610.009 

-1013 

0.49 

3.30 

-0 . 033±0 . 007 

0.00410.008 

-1213 

0.29 

3.40 

-0 . 045±0 . 006 

0.01410.008 

-9±4 

2.36 

3.50 

-0.069±0.006 

0.00010.007 

-6±4 

1.16 

2028  0 

5/2 

1/2 

3.10 

0.34510.016 

-0.00110.021 

-417 

1.31 

3.20 

0.252±0.009 

-0.02210.012 

8±8 

3.96 

3.30 

0.26410.009 

-0 .02610.011 

6±4 

1.47 

3.40 

0.26510.008 

-0.05010.010 

4±4 

1.13 

3.50 

0.29610.009 

-0.05710.011 

4±4 

0.81 
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levels  have  been  included  in  the  statistical  treatment, 

29  29 

the  SCN  model  becomes  applicable.  In  the  Si(p,p'y)  Si 
case,  sufficient  levels  required  an  energy  width  corres¬ 
ponding  to  400  keV  in  the  compound  nucleus. 

Gilbert  et  al  (Gi  65)  reported  level  density 
50  MeV  ^  for  the  mass  A  =  30  region.  Using  their  level 
density  formula  and  the  level  information  of  20P  (En  73) , 
a  level  density  of  ^40  MeV  ^  was  obtained  for  the 
excitation  energy  considered.  Hence  the  energy  of  400 
keV  would  correspond  to  ^20  levels  (the  lower  limit  for 
successful  statistical  averaging) . 

If  one  repeats  the  calculation  for  the  compound 
3  8 

nucleus  of  Ar,  one  obtains  a  level  density  of  ^1000 

MeV~^.  Thus  to  excite  ^20  levels  in  the  compound  nucleus, 

one  would  require  a  target  thickness  of  only  ^20  keV. 

From  proton  range  calculations  (Ja  66),  the  evaporated 

target  (500  pg/cm2  PbCl2)  used  in  the  37C1 (p,ny ) 37Ar 

experiment  has  a  corresponding  thickness  of  ^25  keV. 

Hence  one  is  assured  that  sufficient  compound  nuclear 

levels  were  excited  to  justify  the  statistical  averaging 
3  7 

in  the  Ar  experiment. 

The  agreement  between  previously  established 

29 

results  and  the  SCN  model  predictions  in  both  the  Si 
3  7 

and  Ar  studies  establishes  the  wide  applicability  of 


the  SCN  model. 
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APPENDIX  A 


Modification  of  F(t)  due  to  Cascade  Feeding 

In  lifetime  investigations,  if  the  level  of 
interest  is  fed  not  only  by  direct  population  but  also 
by  cascade  y  rays  from  higher  energy  levels  then  the 
observed  F (t )  will  have  to  be  appropriately  compensated. 
To  simplify  the  calculations,  consider  a  case  where  only 
a  single  cascade  is  involved  (Figure  Al) . 


Compound 

nucleus 


Figure  Al 


n^  and  n^  represent  the  direct  population  of  the  states 
with  lifetimes  and  x2  respectively  and  they  are 
normalized,  (i.e.  n^  +  =  1) 

At  any  time,  t,  the  population  of  the  lower  state 
due  to  direct  contribution  can  be  written  as 


-t/x 


n,  (t )  , .  ,  =  n,  e 

lv  ' direct  1 


(A .  1 ) 


-t/T, 

where  e  is  the  decay  probability.  Similarly  for  the 


upper  state: 

n~  (t) 


-t/x 


0  w  , .  ,  =  n0  e 

2  direct  2 


(A. 2) 
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The  decay  rate  of  the  upper  state  which  equals 
the  cascade  feeding  rate  of  the  lower  state  is 


dt  [n2 (t) direct3  ~ 


1  t/T2 
t 2  n2  6 


- —  n0  (t )  , . 
t 2  2  direct 


(A.  3) 


From  the  laws  of  successive  decay  it  can  be  shown  that 
the  decay  rate  of  the  lower  state  due  to  cascade  feeding 


is 


dt  ^nl direct3  t n2^t3 direct  x^  nl^t3 cascade 


(A. 4) 


in  which  the  terms  represent  the  cascade  feeding  rate 
minus  its  decay  rate.  Equation  A. 4  can  be  integrated 
to  yield: 


n9  t9  -t/x 

n,  (t)  ,  =  -----  [e 

1  cascade  T2“ti 


(A. 5) 


Combining  above  with  Equation  A.l,  the  instantaneous 
population  of  the  lower  state  becomes: 


nl (t) total  nl (t) direct 


~t/Tl  T1 

ni  e  +  n2  ^x  2~t  i 


+  nl ^ cascade 

-t/T  -t/T 

)  [e  ^  -  e 


]  . 


(A. 6) 


The  attenuation  factor  F (x )  for  directly  popula¬ 
ted  level  with  lifetime,  x^,  and  initial  velocity,  vq , 
has  been  defined  (Section  2.1)  as: 


. 
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(T.) 


1 


V  T  . 


°  3 


oo 

Vi(t)  COS 

d 

0 


-t/T 

e 


j 


dt 


(A. 7) 


where  v^(t)  is  the  instantaneous  velocity  and  <j>^  (t)  is 
the  angle  between  the  initial  and  instantaneous  recoil 
directions.  By  abbreviating 


vi(t)  cos  ^(t)  =  Vi(t) 


we  obtain 


F.  ( t  . ) 
1  3 


Ou 


0 


-t/T  . 


V. (t)e 

l 


V.  (0)  T. 
13 


dt 


(A. 8) 


Now  if  one  combines  Equations  A. 6  and  A. 8  with 


the  definition  of  F(t)  one  gets 

-t/T,  T 


CO 

f 


obs 


0 


1  'I  't/T2  -t/TX 

[n1e  xV1(t)  +  n2  ^  _T  (e  -e  )  V±(t)  J 

2  1 


n1F1(T1)+  n2 


T2^*2  ^ t2 ^  -T1^2  ^  ^ 

T2  "  T1 


V,  (0) 

v1  (0) 


(A. 9) 


which  reduces  to 


F(Tl)obs  "  nlFl(Tl)  +  n2 


t2F  (t2)  "  TiF  <Ti) " 

T2  "  T1 


(A. 10) 


when  the  initial  velocities  are  the  same. 

To  a  first  approximation  the  cos- (jTTTtT  term  may 
be  neglected  and  the  instantaneous  velocity  can  be  written 

as 


1 


.•  .  '  .  . . 


(8. A) 
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v(t)  =  v  e  t//a 
o 


(A. 11) 


where  a  is  the  slowing  down  time.  Upon  integration,  the 
attenuation  factor,  F(x),  becomes 


F  (t) 


a 

a  +  t 


(A. 12) 


Substituting  Equation  A. 12  into  Equation  A. 10,  one 
arrives  at 


FTt) 


obs 


n1F(T1) 


n2F(T1) 


F(t2) 


(A. 13) 


In  this  simplified  form,  the  F(t^)  factor  for  the 
lower  state  may  be  extracted  if  one  knows  the  feeding 
fractions  (n^  and  n^)  and  the  FCt^)  of  the  upper  level. 
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APPENDIX  B 


MAGPOP  -  A  Program  to  Calculate 
Angular  Distribution  Coefficients 
of  y  Rays  in  Terms  of  Magnetic 
Population  Parameters. 


153 


154 


MAGPOP  t 

E.R.  Wong  and  B.O.  Aaquist 

Nuclear  Research  Centre 
University  of  Alberta 
Edmonton,  Alberta,  Canada 

Internal  Report  UAE-NPL-47 


August,  1972 


t 


Work  supported  in  part  by  the  Atomic  Energy  Board  of  Canada 


MAGPOP 


E.R.  Wong  and  B.O.  Aaquist 


ABSTRACT 

This  program  is  written  in  FORTRAN  ti  for  the  XDS  920 
computer  at  the  Nuclear  Research  Centre  of  the  University  of 
Alberta.  It  calculates  the  angular  distribution  coefficients, 
Ap,,  of  y-ray  decay  between  levels  with  spins  and  J 2  for  all 
possible  combinations  of  the  magnetic  population  parameter 
P(M  )  {  “Jj  <u  <.  Jj  }  and  mixing  ratio,  6.  A  least  square 
fit  of  the  calculated  coefficients  to  the  measured  y-ray 
angular  distribution  is  included. 


THEORY 


The  general  y-ray  angular  distribution  between  levels  with  spin  J 
and  J2  can  be  expressed  as 

W (0)  =  l  AkPk(cos  0) 

K  even 


•j  ***  ■■  -J*  '  >:  -  .  :  ' 


If  is  aligned  [ie.  P(Mj)  =  P(-M  )]  and  we  consider  only  the  lowest 
two  multipolarities  (usually  the  case),  L  and  L1  =  L  +  1  then  one  can  write 


„  1(ll  [  Rk(LLJ,J2)  +  26Rk(LL'J1J2)  +  <S2Rk(L-L'J1J2)  ]  * 

K  K<  l)x  [  (77T7) - J  <2> 

where  ^  (21  +  1 ) i  <J1|[TL.  I|J2> 

(2L'  +  1)4  <  Jill  tl  I  I  J2  > 

Rk(LL'JiJ2)  =  (-)'  +  Ji'J2  +  L''L'K  X  |(2Jl  +  I)  (2L  +  I)  (2L 1  +  I )] 
x  (LL ' 1-1  I KO)  W  (JjJjLL1 ;KJ2) 

BK(Ji>  =  l  PK(J,M1)  P  <M,> 

Mi 

PK(J lMi )  “  <2J1  +  (J1J,H1-M|  KO) 

6  is  the  mixing  ratio  of  the  higher  to  the  lower  multipole. 

R^'s  and  p^'s  are  angular  momentum  coupling  coefficients  and  have  been 
tabulated  (RO67) . 


The  reaction  dependent  term  is  the  P(Mj),  Frequently,  compound  Nuclear 
Statistical  Model  has  been  used  to  determine  it  in  terms  of  the  Hauser 
Feshbach  penetrability  factor,  T,  (SH66).  In  fact,  computer  codes  like 
Mandy  (SH69)  and  AK(GR71)  have  been  written  to  facilitate  the  calculation. 


*  Note:  phase  convention  of  Rose  and  Brink 


O 
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DESCRIPTION  OF  PROGRAM 

In  some  cases  the  condition  of  using  the  C.N.  Statistical  Model 
may  not  be  fully  justified  and  one  may  wish  to  perform  a  model  independent 
analysis.  This  program  has  been  written  to  achieve  this  model  independency. 
For  a  given  y-ray  decay  (ie.  -+  J2)  it  varies  the  P ( M i ) 1 s  to  calculate 
the  B^'s  and  A^'s  for  any  given  6  or  range  of  6.  It  also  does  a  least 
square  fit  of  the  theoretical  angular  distribution  resulting  from  these 
calculated  A^'s  to  the  experimental  measurements. 

The  program  allows  the  user  to  specify: 

i)  range  of  lowest  P ( M a )  ie.  P(M1=^)or  PCM^O) 

ii)  steps  in  which  the  P(M})  are  to  be  varied 

i i i )  range  of  5 

iv)  steps  to  vary  6 

v)  option  to  output  all  calculated  vales  or  those  that  fall 
below  certain  X2 (to  be  input  by  user) 

vi)  option  to  input  fixed  values  of  P(Mj) 

PROGRAM  INPUT 

Program  input  cards  are  of  the  form  (see  G R 7 1 ) 

ICONTR,  (DUM(I),  1  =  1,7)  FORMAT  ( 1 2 ,  8X,  7F10.0) 

where  ICONTR  dictates  the  type  of  input  card  or  those  that  follow. 

The  program  will  remember  all  quantities  and  perform  calculations  with 
them  until  they  are  changed.  Default  values  has  also  been  incorporated 
into  the  program  such  that  they  wi 1 1  be  used  in  case  the  user  did  not 
input  them. 
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I CONTR  DUM(I)  VARIABLE 
I  = 


COMMENT 


-1 


PROGRAM  EXITS 


0 

1  1 


2  1 


3  1 

2 

3 

k 

5 

6 

7 


NRHO 


NRK 


PM  I  N 

PMAX 

PST ( 1 ) 
PST (2) 
PST (3) 
PST(4) 
PST (5) 


The  following  card  is  title  card 
(FORMAT  20AA) 


Number  of  RHO  coefficients 

RHO  coefficients  follow  on  the  next 
card  or  cards  (FORMAT  8F10.0)  in  the 
order  P2(J1  ~  Jj),  P2(J1  "  Ji  +  !)••• 

P 2  ( J i , J i )  for  ascending  K. 

Number  of  R^  coefficients 

R  coefficients  follow  on  next  and  or 
cards  (FORMAT  8F10.0)  in  the  order 
Rk(LLJ1J2)>  Rk(LL'J1J2),  Rk(L'L'J1J2) 

for  ascending  K. 


smallest  value  of  lowest  magnetic  population 
parameter 

largest  value  of  lowest  magnetic  population 
parameter 


value  to  step  lowest  population 
parameter 


value  to  step  2nd  lowest  population 
parameter 

value  to  step  3rd  lowest  population 
parameter 

value  to  step  **rth  lowest  population 
parameter 

value  to  step  5th  lowest  population 
parameter 


program 
will  read 
(NRHO | 2) 
of  these. 
Extra  cards 
to  input 
these  may  be 
used  if 
necessary . 


etc . 


* 


■ 


. 


ICONTR  DUM(l) 


VARIABLE 

NAME 


COMMENT 


4 

1 

ANORM 

NORMALIZATION  FACTOR  FOR  EXPERIMENTAL  ANGULAR 
DISTRIBUTION  (ie.  A  ) 

2 

SENSIT 

The  program  will  print  values  for  which  the 

X2  falls  below  this. 

Calculation  begins  after  this  card. 

5 

1 

NANG 

Number  of  experimental  angles  (NANG  <.  20) 
Angles  are  input  on  following  card  or  cards 
(FORMAT  8 F 1 0 . 0 ) 

6 

Y  ( 1 ) 

(1=1 , NANG) 

The  next  card  accepts  the  experimental  angular 
distribution  in  the  order  of  the  angles  input 
earl ier  (FORMAT  8F10.0) 

7 

- 

Ml) 

The  next  card  accepts  exp.  dist.  error 
(FORMAT  8F10.0) 

8 

1 

DM  IN 

MINIMUM  DELTA  values 

2 

DMAX 

MAXIMUM  DELTA  values 

3 

DSTEP 

value  to  step  the  delta  values 

9 

1 

2 

Q2 

Attenuation  coefficients 

10 

1 

1  OPT  ION 

=1  program  will  print  output  only  if  the 

X2  is  below  the  sensitivity 

5^1  program  will  print  all  output 

11 

1 

NPOP 

number  of  fixed  P(M1)  values  to  follow  on 

the  next  card  (FORMAT  8F10.0) 
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DEFAULT  VALUES 


Default  values  are  used  unless  overriden  by  input. 


DESCRIPTION 

VARIABLE 

Mixing  ratio,  6 

DMIN  =  0.0 

DMAX  =0.0 

DSTEP  =  30.0 

Magnet i c 

Popu 1  at i on 
parameter 

PM  1 N  =  0.0 

PMAX  =0.5 

PST ( 1 )  -  0.1  for  1  =  1,10 

Attenuat i on 
coef f i ci ent 

Q2  =  1.0 

QA  =  1.0 

Experimental 

angles 

NANG  =5.0 

THETH(l)  =  0.0 

"  (2)  =  31.0 
"  (3)  =  55.-0 

"  (A)  =  71.0 

"  (5)  =  90.0 

Output 
opt  ion 

1  OPT  ION  =1.0 

SENSIT  =  A. 6 

COMMENTS 

The  program  is  limited  to  calculation  for  initial  state  of  spin, 

Jj  <.  19/2  and  for  expansion  of  angular  distribution  up  to  the  order  K 

This  latter  condition  restricts  unique  population  combination  for  each 

value  of  B  for  J  5/2. 

K.  * 
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WAVY 

E.R.  Wong 


Abstract 

A  FORTRAN  IV  program  written  for  expanding  experimental  angular 
distributions  into  Legendre  polynomial  series  of  any  order  (<  15)  in 
all  or  even  orders  by  the  least  square  method.  It  also  allows  one 
to  average  simultaneously  several  (  <20)  excitation  functions  over 
any  range  of  energy  interval.  It  is  the  intention  of  this  code  to 
use  both  these  functions  for  the  analysis  of  yield  curves  taken  at 
several  angles.  However,  each  of  the  aforementioned  subroutines  may 
be  executed  separately. 

Program  Description 

Differential  cross  sections  (or  angular  distributions)  are  fre¬ 
quently  expanded  into  the  Legendre  polynomial  series 

w(ei)  =  l  akPk  (cos  ei) 

•where 

W(0.)  =  intensities  at  angle  6^ 

P.  ( e - )  =  Legendre  polynomial  of  order  k 

K 

ak 


coefficient  of  order  k 


The  program  uses  the  method  of  least  square  to  fit  an  experi¬ 
mental  angular  distribution  to  these  polynomials.  As  the  procedure 
is  well-known  and  has  been  described  in  detail  elsewhere  *)  only 
a  brief  outline  will  be  given  here. 

If  the  measured  value  is  represented  by  V ( 6 ^ )  then  the  method 
of  least  square  specifies  that  this  be  achieved  by  chosing  a^  so 
that 

a  =  I  w.  [m  (e.)-v  (e.)]2 

is  a  minimum,  is  the  weighting  factor  and  generally  expressed 

in  terms  of  standard  deviation,  o.,  by  w.  =  —  . 

1  J  l  a. 

Since  the  errors  of  measurement  (i.e.  standard  deviation,  a. ) 
are  known,  then  A^-p  has  a  x2  distribution  with  (n  -m)  degrees  of 
freedom  where  n  is  the  number  of  observations  and  m  the  number 
of  fitted  parameters. 

Given  a  set  of  observations  (angular  distributions)  the  program 
will  determine  the  a^  coefficients.  The  maximum  order,  kmax,  of 
these  coefficients  is  governed  by  the  user.  However,  kmax  must  not 
exceed  the  number  of  observations  or  15.  An  option  is  also  provided 
in  cases  where  only  even  order  Legendre  polynomial  expansion  is 
required,  (e.g.  reactions  that  exhibit  symmetry  about  90°).  Aside 
from  printing  these  a^  coefficients  and  their  associated  errors,  the 
program  will  also  give  the  x2  Per  degrees  of  freedom  for  the  fit. 

The  other  part  of  the  program  concerns  the  averaging  process 
whereby  an  excitation  function  recorded  with  good  energy  resolution 
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may  be  averaged  to  simulate  yield  with  poorer  resolution.  In  particu¬ 
lar  by  averaging  two  adjacent  yield  points  one  would  simulate  an 
excitation  obtained  with  a  target  twice  as  thick.  Consequently, 
increasing  the  averaging  interval  would  correspond  to  thickening  the 
target.  Of  course  one  may,  on  the  other  hand,  regard  it  as  deterior¬ 
ation  of  beam  energy  resolution.  The  program  allows  the  user  to 
impose  lower  and  upper  limits  as  well  as  the  increment  of  the  averaging 
process. 

The  averaging  procedure  may  be  combined  with  the  least  square 
fitting  procedure  to  analyse  yield  curves  obtained  at  several  angles 
simultaneously.  Thus  the  program  will  perform  the  least  square  fit 
after  each  averaging  interval. 

The  program  consists  of  a  mainline  which  controls  the  input  and 
output  as  well  as  the  calling  of  subroutines.  It  also  contains  the 
option  of  writing  the  output  into  a  disk  file  beside  the  normal 
printout.  This  latter  option  is  handy  when  one  wishes  to  plot  the 
yield  curves.  The  format  with  which  the  output  is  written  onto  the 
file  is  compatible  with  the  input  format  of  the  plotting  code, 
"AUTOPLOTTER"  14 ). 


Input 

A  control,  ICONTR,  is  used  to  distinguish  between  the  different 
input  quantities  and  dummy  variables  are  then  used  to  accept  all 
input  2»3).  The  input  cards  are  of  the  form 


ICONTR,  (DUM(I),  I  1,7)  FORMAT  ( 1 2 , FI 8 . 0 , 6F1 0 . 0 ) 


•  ■ 
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Once  a  quantity  has  been  entered  it  will  be  used  until  changed  with  the 
exception  of  ICONTR  =  2  card  which  must  be  included  for  each  set  of 
averaging  calculation,  otherwise  the  default  value  will  be  used.  The 
codes  for  control  and  dumny  variables  are  tabulated  in  Table  I  and  the 
default  values  for  some  variables  in  Table  II. 

As  mentioned  earlier  the  two  parts  of  the  program  may  be  executed 
independently;  the  key  control  cards  for  calling  the  respective  sub¬ 
routines  are  ICONTR  =  2  and  ICONTR  =  5.  If  ICONTR  =  2  card  is  omitted, 
then  the  program  will  not  call  the  averaging  routine.  On  the  other  hand 
if  on  the  ICONTR  =  5  card,  the  quantity,  NANG  =  1,  then  no  least  square 
fitting  will  be  performed. 

In  the  averaging  process,  the  number  of  yield  points  are  determined 
from  the  maximum,  minimum  and  increment  and  hence  the  ICONTR  =  2  card 
must  precede  ICONTR  =  4,6  and  7  cards.  Yields  for  each  energy  are  then 
expected  to  be  entered  on  separate  cards  in  order  of  increasing  energy 
following  ICONTR  =  6  card.  The  normalization  factors  (ICONTR  =  4)  and  yiel 
errors  (ICONTR  =  7)  are  input  in  the  similar  manner.  All  quantities  in 
general  may  be  input  in  any  sequence. 

The  program  will  start  executing  after  the  ICONTR  =  7  series  of 
cards  have  been  read  and  will  continue  calculating  for  any  number  of 
sets  of  data  until  it  encounters  ICONTR  =  -1  or  other  terminating 


commands. 
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TABLE  I 


ICONTR 

DUM(I) 

I  => 

VARIABLE 

COMMENT 

-1 

Stop 

0 

TITLE(I) 

Title  read  on  the  following  card 
(FORMAT  20A4) 

1 

1 

KB 

Maximum  order  of  Legendre  polynomial  to 
be  fitted  (KB  <  15) 

2 

NOPT 

0  only  even  orders  fitted 

1  all  orders  fitted 

3 

NPLOT 

0  no  output  for  plot 

1  output  for  plot  without  errors 

>1  output  for  plot  with  errors 

2* 

1 

ENMIN 

Minimum  energy 

2 

ENMAX 

Maximum  energy 

- 

3 

ENDEL 

Energy  increment 

4 

XINTMI 

Minimum  energy  interval  to  be  averaged 

5 

XINTMA 

Maximum  energy  interval  to  be  averaged 

6 

XINTDE 

Energy  interval  increment 

3** 

Q(I) 

Attenuation  coefficients  follow  on  the 
next  card  or  cards  (FORMAT  8F10.0) 

4 

TEMP(I) 

Normalization  factors  follow  on  the 
next  card  or  cards  (FORMAT  8F10.0) 

5 

1 

NANG 
(  <  20) 

Number  of  experimental  angles 

6 

YY(I) 

Experimental  angular  distributions  read 
in  the  order  of  angles  on  the  next  card 
or  cards 

7 

YEY(I) 

Experimental  distribution  errors  cor¬ 
responding  to  YY(I)  on  the  following 
card  or  cards 

Program 

starts  execution  after  ICONTR  =  7  series  of  cards  are  read. 

*  This  card  must  precede  ICONTR  =  4,  6  and  7  cards  if  it  is  used. 

**  This  card  must  follow  ICONTR  =  1  card  if  it  is  used. 
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TABLE  II 


DEFAULT 

VARIABLE 

DESCRIPTION 

3 

KB 

Maximum  number  of  order  of  the 

Legendre  polynomial 

0 

NOPT 

Option  to  fit  only  even  orders  of 
Legendre  polynomial 

0 

NPLOT 

Option  of  NO  output  for  plot 

1 

NANG 

Number  of  experimental  angles  =  1 

1  .0 

Q(l) 

Attenuation  coefficients  =  1.0 

1.0 

TEMP(I) 

Normalization  factors  =  1.0 
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Neutron  and  y-ray  angular  distributions,  and  y-ray  polarisation  correlations  have  been  measured  with 
the  reaction  40  Ar(p,  ny)40K.  Positive  parity  states  have  been  established  at  1,644  MeV  (0+),  1,959 
MeV  (2+ ),  2.261  MeV  (3+)  and  2.290  MeV  (1> ). 


Recent  interest  in  deformed  states  in  nuclei 
around  4^Ca  [1,2]  led  us  to  search  for  such 
states  in  40k,  where  they  would  be  character¬ 
ised  by  their  positive  parity.  The  39K(d,p)40K 
reaction  [3]  was  known  to  preferentially  excite 
the  single  p-h  (negative  parity)  states,  the  four 
lowest  states  at  0,  30,  800  and  891  keV  are 
fi  di*  states  whilst  those  at  2047,  2070,  2103  and 
262 ^  keV  arise  from  the  pj  dl^  configuration.  We 
have  studied  the  40Ar(p,n)^0K  reaction  which 
strongly  excites  nearly  all  the  low-lying  levels 
in  40k  including  the  four  lowest  levels  not  de¬ 
scribed  by  single  p-h  configurations.  These  are 
at  1644,  1959,  2261  and  2290  keV. 

The  lifetime  of  the  1644  keV  state  of  0.49  rs 
[4]  together  with  its  mode  of  decay  to  both  the 
800  keV  (2")  and  30  keV  (3")  states  indicate  a 
probable  assignment  of  0+.  On  the  other  hand,  a 
negative  parity  with  J  ^  2  is  favoured  by  a  recent 
DWBA  analysis  of  data  from  the  40Ar(T,t)40K 
reaction  [5].  We  have  measured  the  neutron  an- 
ular  distributions  from  the  five  lowest  levels  in 
oK  at  a  bombarding  energy  of  5.5  MeV  using  the 
pulsed  proton  beam  and  neutron  time -of -flight 
facility  at  the  University  of  Alberta.  The  only 
distribution  to  exhibit  any  marked  anisotropy 
was  that  for  neutrons  exciting  the  1644  keV  state. 
The  symmetry  about  90°  of  this  distribution, 
shown  in  fig.  1,  indicates  the  dominance  of  the 
compound  nuclear  mechanism.  The  theoretical 
distributions  shown  were  calculated  with  a 
Hauser -Feshbach  computer  program  [6]  which 
uses  the  average  potentials  of  Rosen  [7].  Small 
changes  in  the  potentials  did  not  affect  the  shape 

*  Work  supported  in  part  by  the  Atomic  Energy  Con¬ 
trol  Hoard  of  Canada. 

**  On  leave  of  absence  from  University  of  Liverpool, 
England. 


Fig.  1.  Neutron  angular  distribution  for  1644  keV 
state  in  the  reaction  40Ar(p,n)40K  at  a  proton  bom¬ 
barding  energy  of  5.5  MeV.  The  fitted  curves  are  de¬ 
rived  from  a  compound  nuclear  model  calculation  for 
spin  assignments  of  0,  1  and  5.  The  other  possible 
spin  values  of  2,  3  and  4  give  essentially  isotropic  dis¬ 
tributions. 

of  the  distributions  which  clearly  indicate  the 
1644  keV  state  has  spin  0.  The  negative  parity 
assignment  is  ruled  out  as  it  requires  an  M3  en¬ 
hancement  factor  of  104  for  the  1614  keV  y-ray. 

The  1959  keV  state  was  studied  at  a  proton 
energy  of  4.44  MeV,  90  keV  above  its  threshold 
excitation  energy.  We  measured  the  angular  dis¬ 
tributions  of  both  the  1159  keV  (—800  keV  2~ 
state)  and  1929  keV  (—30  keV  3“  state)  y-rays 
with  a  45  cm3  Ge(Li)  detector;  the  y-y  correla¬ 
tion  between  the  1159  keV  and  770  keV  y-rays 
with  Nal  crystals;  and  the  polarisation  direction 
correlation  of  the  1159  keV  y-ray  with  a  Comp¬ 
ton  polarimeter  with  the  45  cm3  Ge(Li)  counter 
as  the  scatter er  and  two  Nal  crystals  as  the  ana¬ 
lysers.  The  data  were  fitted  in  terms  of  the 
1959  keV  level’s  magnetic  substate  populations, 
which  were  limited  to  lie  within  values  calcu¬ 
lated  at  threshold  and  200  keV  above  threshold 
by  the  MANDY  program  of  Sheldon  and  Van  Pat¬ 
ter  [8],  which  is  based  on  the  compound  nuclear 
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model.  The  large  «2  Legendre  coefficients  of 
0.40  ±  0.01  and  -0.22  ±  0.02  for  the  angular  dis¬ 
tributions  of  the  1159  keV  and  1929  keV  y- rays, 
respectively,  could  only  be  fitted  by  a  spin  2 
assignment  to  the  level,  which  was  confirmed  by 
the  y-y  correlation.  The  polarisation  of  the 
1159  keV  y- ray  was  measured  at  angles  of  9  = 

=  90°,  45°  and  0°  relative  to  the  beam  direction, 
the  results  at  0°  being  used  to  correct  for  slight 
differences  in  the  efficiency  of  the  analysing 
crystals.  The  definition  of  polarisation  used  was 

b  1  N(0,  0  =  0)-  N(0,  <p  =  90) 

1  R  N{9,  0  =  0)  +  N(9,  0  =  90) 

where  N(9,  0)  is  the  number  of  counts  in  the 
photo-peak  of  the  sum  coincidence  spectra  at  an¬ 
gles  0,  0.  The  efficiency,  R,  of  the  polarimeter 
was  calculated  with  the  solid -angle  correction 
technique  suggested  by  Ferguson  [9]  and  checked 
by  the  results  for  the  2+  -*  0+  1460  keV  transi¬ 
tion  in  The  measured  polarisation  of 


Fig.  2.  A  plot  of  x2  versus  6  for  various  spin  and  par 
ity  assignments  to  the  1959  keV  state.  The  angular 
distributions  of  the  1159  keV  and  1929  keV  y-ray  and 
the  polarisation-direction  correlation  of  the  1159  keV 
y-ray  are  fitted  simultaneously. 


40 
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Fig.  3.  A  summary  of  the  spin  and  parity  assignments 
together  with  the  y-ray  branching  ratios  for  the  low- 
lying  states  in  40^ 


-0.76  ±  0.21  (the  efficiency  calculation  contrib¬ 
uting  most  of  the  error)  enabled  the  definite  as¬ 
signment  of  positive  parity  to  the  1.959  MeV 
state.  Thex^  plot  versus  6,  shown  in  fig.  2,  for 
various  spin  and  parity  assignments  indicates 
pure  El  character  for  the  1159  keV  y-ray.  The 
1929  keV  y-ray  was  determined  to  have  anM2/El 
mixing  ratio  of  0.10  ±  0.07. 

From  similar  techniques  the  angular  distri¬ 
butions  of  the  2231  keV  and  2261  keV  y-rays  in¬ 
dicate  that  the  2261  keV  state  has  a  spin  of  3 
whilst  the  polarisation  correlation  at  9  =  90°, 
yielding/)  =  -0.57  ±  0.28,  shows  that  the  state 
has  positive  parity.  The  2290  keV  state  observed 
in  the  39K(d,p)40K  reaction  [3]  has  been  shown 
to  be  a  doublet  of  energies  2290  and  2291  keV. 
The  2290  keV  state  decays  to  the  1644  keV  0+ 
state  (646  keV  y-ray)  and  the  800  keV  2"  state 
(1490  keV  y-ray)  whilst  the  2291  keV  state  de¬ 
cays  to  the  4'  ground  state.  The  angular  distri¬ 
bution  of  the  646  keV  y-ray  can  only  be  fitted 
with  a  spin  assignment  of  1  and  the  compound 
nuclear  calculation  also  gives  a  fit  three  times 
better  for  1+  than  for  1".  The  presence  of  the 
strong  decay  mode  to  the  0+  state  also  indicates 
the  state  has  positive  parity.  The  angular  dis¬ 
tribution  of  the  2291  keV  y-ray  with  an  a.’i  of 
0.55  ±  0.03  can  only  be  fitted  for  a  spin  of  4  to 
the  2291  keV  state. 

These  results,  summarised  in  fig.  3,  show 
that  at  least  four  positive  parity  states,  with 
possible  configuration  f?di2  and  with  spins  0+, 
1+,  2+  and  3+,  lie  at  the2 same  energy  in  40^  as 
the  single  p-h  pj  dl*  states. 
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Abstract:  The  mean  lifetime  of  levels  in  37Ar  below  3.30  MeV  excitation  were  measured  via  the 
37CI(p,  ny)37Ar  reaction  with  the  Doppler-shift  attenuation  method.  Gamma-ray  angular 
distributions  were  obtained  at  various  incident  energies  and  compared  with  the  predictions  of 
the  statistical  compound  nuclear  model.  The  angular  distribution  and  lifetime  results  yield  an 
assignment  of  Jn  ji(  +  )  for  the  3.171  McV  level  and  J n  !p  for  the  3.185  MeV  level.  The 
experimental  results  are  discussed  and  compared  with  shell-model  calculations. 


NUCLEAR  REACTION  37CI(p,  ny),  £’  =  3.25-5.90  McV;  measured 
o(£;  Ey,  0y),  Doppler-shift  attenuation.  37Ardedueed  levels, r,y-branching  ratios, 
Jn,  y-rrixing.  Enriched  targets;  Neutron  time-of-flight  technique;  Ge(Li)  detector. 


1.  Introduction 

Nuclei  near  doubly  closed  shells  are  of  considerable  interest  due  to  the  relatively 
simple  shell-model  configurations  expected,  and  experimental  data  can  be  readily 
compared  in  detail  with  theoretical  predictions.  However,  the  fact  that  the  simple 
shell  model  cannot  account  for  even  the  number  of  levels  in  37Ar  up  to  an  excitation 
of  4  MeV  [ref.  ’)],  suggests  the  presence  of  more  complex  configurations.  The  mea¬ 
surements  of  other  properties  such  as  the  electromagnetic  decay  and  lifetimes  should 
provide  useful  information  about  these  configurations. 

The  present  experiment  was  undertaken  to  provide  information  on  level  spins, 
multipole  mixing  ratios,  branching  ratios  and  lifetimes  for  37Ar  levels  up  to  an  exci¬ 
tation  of  3274  keV.  The  results  are  compared  with  previous  information,  where  avail¬ 
able.  The  existence  of  the  3185  keV  level  1  ■ 2)  is  directly  confirmed  by  a  neutron  time- 
of-flight  measurement.  The  lifetime  and  mixing  ratio  results  are  discussed  in  terms  of 
the  available  theoretical  predictions. 

2.  The  experiment 

2.1.  EXPERIMENTAL  ARRANGEMENT 

2.1.1.  Lifetime  determination.  Levels  in  3  ,Ar  up  to  an  excitation  of  3.28  MeV 
were  populated  via  the  37CI(p,  ny)37Ar  reaction  ( Q  =  -1.596  MeV)  using  a  5.70 

t  Present  address:  Nuclear  Reactions  Section,  Tata  Institute  of  Fundamental  Research,  Homi 
Bhabha  Road,  Colaba,  Bombay,  5  India. 

+t  This  work  has  been  supported  in  part  by  the  Atomic  Energy  Control  Board  of  Canada. 
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MeV  proton  beam  from  the  University  of  Alberta  6  MV  Van  de  Graaff  generator. 
The  target  consisted  of  500  gg/cm2  of  PbCI2  (enriched  to  87  %  in  37C1)  evaporated 
onto  a  0.5  mm  thick  Au  backing.  Gamma  rays  were  detected  both  in  coincidence  with 
neutrons  and  in  singles  by  a  45  cm3  Ge(Li)  detector  placed  10  cm  from  the  target. 
Neutrons  were  detected  in  a  5.1  cm  thick  x  (17.8  cm  OD,  2.5  cm  ID)  NE218  annular 
neutron  detector  subtending  angles  between  1 57°  and  177°  to  the  beam  axis  and  placed 
15.8  cm  from  the  target.  The  y-flux  in  the  neutron  detector  was  suppressed  by  a  1.95 
cm  thick  lead  shield.  Detected  y-rays  were  rejected  by  time-of-flight  gating  and  pulse- 
shape  discrimination.  Details  of  the  y-rejection  and  n-y  coincidence  circuits  are  given 
in  ref.  3). 

Gamma-ray  peaks  from  88 Y  and  60Co  sources  were  accumulated  during  all  runs. 
In  singles  runs  these  were  added  directly  to  the  spectra,  but  in  the  n-y  coincidence  runs 
additional  circuits  were  used  to  select  y-y  coincidences  between  the  Ge(Li)  detector 
and  a  7.5  cm  x  7.5  cm  Nal(Tl)  crystal  placed  nearby.  This  y-y  coincidence  source 
spectrum  was  accumulated  in  parallel  with  the  n-y  spectrum  and  provided  a  contin¬ 
uous  monitor  of  the  system  stability. 

2.1.2.  Angular  distribution  measurements.  Gamma-ray  distributions  for  levels  up 
to  an  excitation  of  3274  keV  were  measured  in  singles  at  several  proton  energies  be¬ 
tween  3.25  and  5.15  MeV.  The  Ge(Li)  detector  was  mounted  on  a  rotating  trolley  at 
a  distance  of  20  cm  from  the  target.  Measurements  were  carried  out  using  thick  tar¬ 
gets  [see  ref. 3)]  of  enriched 3  7CI  (»  85  %)  in  the  form  of  KC1  deposited  onto  a  0.05 
cm  thick  Ta  backing  as  well  as  the  evaporated  PbCl2  target.  An  additional  small 
(6  cm3)  Ge(Li)  detector  fixed  at  90°  was  used  in  the  initial  stage  of  the  experiment  to 
monitor  the  reaction  yield. 

2.2.  DATA  ANALYSIS 

2.2.1.  Lifetime  analysis.  Observed  energy  shifts  for  the  various  levels  were  ob¬ 
tained  from  the  centroid  position  of  the  appropriate  y-rays  by  a  least-squares  fitting 
procedure.  The  mean  initial  recoil  velocity  was  determined  from  kinematics,  taking 
into  account  the  energy  loss  of  the  incident  beam  through  the  target.  For  n-y  coin¬ 
cidence  runs,  the  finite  angle  of  the  neutron  detector  was  also  taken  into  account. 
For  the  singles  experiments,  the  influence  of  the  neutron  angular  distribution  on  the 
mean  initial  recoil  velocity  was  accounted  for  using  the  method  of  ref.  4).  Neutron 
angular  distributions  were  also  independently  measured  using  a  time-of-flight  tech¬ 
nique  [see  ref.  5)]  to  explicitly  calculate  the  mean  initial  recoil  velocity;  results  ob¬ 
tained  using  the  two  methods  agreed  in  all  cases  to  better  than  1  %. 

The  F( t)  curves  were  calculated  using  a  computer  code  6)  which  uses  the  formulae 
of  Lindhard  7)  to  calculate  electronic  and  nuclear  stopping.  There  is  evidence  8)  for 
systematic  deviation  from  the  Lindhard  predictions  for  electronic  stopping  of  up  to 
30  %,  and  in  fact  the  stopping  power  calculations  of  Firsov  9)  have  been  shown  in 
many  cases  10)  to  give  closer  agreement  with  experiment  than  those  of  Lindhard. 
However,  for  Ar  ions  stopping  in  Cl,  where  the  largest  deviation  from  the  Lindhard 
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prediction  might  be  expected  in  the  present  experiment,  the  two  predictions  agree  to 
within  1  %.  Consequently  the  Lindhard  values  have  been  used  throughout  and  a  20  % 
error  assigned  to  the  total  stopping  power. 

2.2.2.  Angular  distribution  analysis  and  branching  ratio  determination.  Gamma- 
ray  yields  were  obtained  by  summing  counts  over  the  full-energy  peak  of  the  appro¬ 
priate  y-rays  and  normalizing  to  the  yield  of  the  isotropic  decay  of  the  1410  keV  37Ar 
level  ( Jn  —  i  +  ).  Angular  distributions  were  fitted  to  an  even-order  Legendre  poly¬ 
nomial  series  and  also  compared  with  the  predictions  of  the  statistical  model  of  Shel¬ 
don  n).  In  addition,  comparisons  were  made  allowing  the  spin  substate  populations 
to  vary  in  order  to  evaluate  the  dependence  of  extracted  mixing  ratios  on  the  validity 
of  the  statistical  model  assumption  (subsect.  3.3.1).  The  phase  convention  of  Rose 
and  Brink  12)  for  mixing  ratios  was  used  throughout.  Mixing  ratios  errors  were  estim¬ 
ated  using  the  method  proposed  by  Cline  and  Lesser  13). 

Angular  distributions  predicted  by  the  statistical  model  were  calculated  and  com¬ 
pared  with  experiment  using  the  computer  code  AK  [ref.  14)].  All  kinematically 
allowed  neutron  exit  channels  as  well  as  the  proton  elastic  channel  and  the  four 
strongest  a-particle  exit  channels  were  included  in  the  calculation.  The  necessary 
transmission  coefficients  were  obtained15)  using  the  optical-model  parameters  of 
Rosen  16)  for  protons,  Perey  and  Buck  1 7)  for  neutrons  and  Davis  18)  for  a-particles. 
Transmission  coefficients  were  also  obtained  from  the  code  ABACUS  [ref.  I9)]  and 
yielded  essentially  the  same  results. 

Branching  ratios  were  determined  using  the  zero-order  Legendre  polynomial  co¬ 
efficients  and  the  previously  measured  detector  efficiencies  20). 


3.  Experimental  results 

3.1.  EXCITATION  ENERGIES  AND  DECAY  SCHEMES 

Excitation  energies  were  inferred  from  the  various  y-ray  energies  determined  by  the 
fitting  procedure  described  in  subsect.  2.2.1.  Gamma-ray  energies  were  determined 
independently  from  the  singles  and  the  coincidence  data,  and  found  to  be  in  good 
agreement.  The  excitation  energies  are  summarized  in  fig.  1,  together  with  the  branch¬ 
ing  ratio  information.  The  present  results  are  in  generally  good  agreement  with  those 
of  refs.  lf  2>  21  ~23). 

It  has  been  proposed  from  y-ray  studies  *' 2)  that  the  3.17  MeV  level  is  a  doublet, 
consisting  of  levels  of  3171  and  3185  keV  excitation.  To  obtain  direct  confirmation 
of  the  presence  of  the  3185  keV  level,  a  neutron  time-of-flight  spectrum  was  accumu¬ 
lated  using  a  pulsed  proton  beam  5 )  at  an  incident  proton  energy  of  5.90  MeV  and  a 
100  gg/cm2  thick  PbCI2  target  (enriched  to  87  %  37C1).  A  neutron  group  correspond¬ 
ing  to  the  proposed  level  can  be  clearly  seen  (fig.  2)  in  the  time-of-flight  spectrum, 
giving  direct  confirmation  of  the  existence  of  a  level  in  37Ar  at  an  excitation  of  3185 
keV. 
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Fig.  1.  Energy  level  and  decay  scheme  of  37  Ar  levels  up  to  an  excitation  of  3274  keV.  The  branching 
ratios  and  level  positions  are  those  determined  in  the  present  work;  Jn  assignments  combine 
previous  results  discussed  in  the  text  with  the  present  results. 


CHANNEL  NUMBER 

Fig.  2.  Time-of-flight  spectrum  of  neutrons  from  the  37CI(p,  n)37Ar  reaction  observed  for  a  flight 
path  of  6.37  m.  Neutrons  were  detected  at  0°  to  the  beam  axis  in  a  NE2I3  liquid  scintillator.  The 

incident  proton  energy  is  5.90  MeV. 
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3.2.  LIFETIMES 

The  lifetime  of  levels  up  to  an  excitation  of  2795  keV  were  investigated  using  the 
coincidence  technique  and  up  to  an  excitation  of  3274  keV  in  singles.  Gamma-ray 
energy  shifts  between  0°  and  120°  for  several  levels  are  indicated  in  fig.  3.  The  present 
lifetime  information  is  summarized  in  table  1  and  also  compared  with  previous  results. 
The  attenuation  factors  and  lifetimes  listed  for  the  singles  measurements  have  been 
corrected  for  cascade  feeding  and  all  the  results  include  a  small  correction  ( %  10  %) 
due  to  recoiling  37Ar  nuclei  which  decayed  in  the  Au  backing. 

The  lifetime  of  the  3185  keV  level  determined  in  the  present  experiment,  t  = 
285  ±60  fs,  is  to  be  compared  with  t  =  950  +  450  fs  obtained  by  Ivascu  et  al.  2). 


Fig.  3.  Plots  ofy-ray  energy  versus  cos  0V  for  transitions  in  37Ar.  The  solid  lines  correspond  to  a 
least-squares  fit  to  the  experimentally  observed  energies  and  the  dashed  lines  correspond  to  the  full 
kinematic  shifts.  The  plots  shown  are  for  coincidence  experiments. 
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A  spectrum  of  the  3185  -»  1611  keV  transition  is  shown  in  fig.  4  [cf.  fig.  2  of  ref.2)]. 


t 


*0 


Fig.  4.  Variation  of  energy  with  detection  angle  for  the  1574,  161 1  and  1662  keV  y-rays,  originating 

from  the  3185,  1611  and  3274  keV  levels  in  37 Ar. 

3.3.  ANGULAR  DISTRIBUTIONS 

Gamma-ray  angular  distribution  measurements  were  carried  out  for  levels  up  to 
3274  keV  excitation.  Angular  distributions  measured  at  Ep  =  3.35,  4.55  and  5.15 
MeV  were  done  with  the  PbCl2  target;  the  rest  were  done  with  the  KC1  target.  The 
angular  distribution  of  y-rays  from  the  1410  keV  level  ( J"  =  i  +  ),  used  to  normalize 
all  other  y-ray  yields,  was  independently  measured  at  a  beam  energy  of  3.25  MeV  to 
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be  isotropic  within  3  %  using  the  small  fixed  Ge(Li)  monitor.  Legendre  polynomial 
expansion  coefficients  for  all  measured  angular  distributions  are  listed  in  table  2; 
spin-  and  mixing-ratio  determinations  are  summarized  in  tables  3  and  4. 


Table  2 

Legendre  polynomial  coefficients  for  angular  distributions 


Transition 

Bombarding 

«2 

a* 

(keV) 

energy  (MeV) 

1410  ->  g.s. 

3.25 

-0.01  ±0.03 

0.00  ±0.03 

161 1  ->  g.s. 

3.35 

0.24  ±0.02 

— 0.10±0.02 

3.40 

0.1 8  ±0.03 

-0.07  ±0.03 

4.05 

0.1 8  ±0.01 

-0.03  ±0.01 

4.30 

0.1 6  ±0.01 

-0.03  ±0.01 

4.55 

0.1 7  ±0.01 

-0.08  ±0.02 

2217 ->  g.s. 

4.05 

0.44±0.02 

0.03  ±0.02 

4.30 

0.23  ±0.01 

-0.01  ±0.01 

4.55 

0.22  ±0.01 

— 0.04±0.01 

2488  ->  g.s. 

4.55 

0.05  ±0.01 

0.00±0.01 

5.15 

0.05  ±0.01 

0.02±0.01 

2795  ->  g.s. 

4.55 

-0.01  ±0.03 

0.07  ±0.03 

5.15 

0.01  ±0.01 

0.02  ±0.01 

3171  ->  g.s. 

5.15 

— 0.36±0.08 

— 0.04±0.08 

3185  -*  1611 

5.15 

0.47±0.01 

0.02  ±0.02 

3274  -*•  g.s. 

5.15 

-0.08  ±0.01 

-0.01  ±0.01 

Table  3 

Mixing  ratio  of  the  161 1  keV  ->  g.s. 

transition 

SCN  prediction 

Best-fit  population  distribution 

(MeV) 

determination 

arctan  <5(deg)  M 

*) 

arctan  d  (deg) 

P(M) 

3.35 

8±4 

4 

0.239 

8±3 

0.247 

l 

0.179 

0.173 

4 

0.082 

0.080 

l 

0 

0 

3.40 

12±4 

4 

0.238 

1 3  ±4 

0.233 

t 

0.179 

0.182 

4 

0.083 

0.085 

l 

0 

0 

4.05 

7±4 

l 

0.205 

1 0  ±  3 

0.210 

5 

0.165 

0.162 

i 

0.098 

0.097 

7 

2 

0.032 

0.031 

4.30 

8±3 

.1 

0.202 

8±3 

0.203 

3 

2 

0.163 

0.140 

l 

0.099 

0.115 

l 

0.036 

0.042 

4.55 

8  ±4 

4 

0.199 

1 0  ±  3 

0.214 

l 

0.162 

0.168 

i 

0.100 

0.085 

l 

0.039 

0.033 

')  Populations  are  listed  only  for  positive  substates. 
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Table  4 

Summary  of  spin  and  mixing-ratio  assignments  in  37Ar 


Transition 

Present  results 

Champlin 
ref.  ‘) 
d 

E,  E, 

J 

>) 

2217  ->  g.s. 

l 

-  0.03+0.04 

-0.03+0.04 

2488  ->  g.s. 

3 

2 

unrestricted 

2795  g.s. 

l 

- 0.21  ±0.13 

0.00+0.15 

3171  ->  g.s. 

1 

0  7+1.5 
”  -0.6 

0.3  +0.3 

3185  —  161  1 

2 

—0.58+0.09 

3274  -  g.s. 

-0.03+0.13,  4.0 +  5- 5 

0.20+0.25  or  2.7  +0.5 

3.3.1.  The  1611  keV  level.  Angular  distributions  of  the  1611  keV  ground  state 
transition  were  measured  at  Ep  =  3.35,  3.40,  4.05,  4.30  and  4.55  MeV.  A  y-ray  yield 
function  measurement  [see  ref.  22)],  however,  indicates  that  there  is  some  non-statis- 
tical  contribution  on  the  reaction  yield,  bringing  into  question  the  accuracy  of  statis¬ 
tical  compound  nuclear  (SCN)  model  predictions  in  the  present  experiment. 

The  presence  of  partial  wave  strength  in  addition  to  the  expected  statistical  strengths 
could  change  the  magnetic  substate  populations  of  a  decaying  level  so  that  interpreta¬ 
tion  of  an  angular  distribution  using  an  unmodified  SCN  model  calculation  would 
result  in  an  incorrect  mixing-ratio  assignment  and,  in  extreme  cases,  an  incorrect  spin 
assignment.  Consequently  all  the  1611  keV  angular  distributions  have  been  analysed 
both  with  SCN  predictions  and  by  allowing  magnetic  substate  populations  to  vary 
arbitrarily. 

The  1611  keV  y-ray  angular  distribution  analysis  resulted  in  a  unique  J  =  \  assign¬ 
ment;  this  is  in  good  agreement  with  the  previous  assignment 2,1 22)  of  Jn  =  -J~.  Also, 
the  mixing  ratios  at  all  beam  energies  are  in  good  agreement  (table  3).  The  accepted 
value  for  the  mixing  ratio  (S  =  0.16  +  0.04)  is  in  good  agreement  with  the  value  of 
Champlin  et  at.  *),  0.1 5 ±0.04,  Ragan  et  al.  21 ),  0.09 ±0.03  and  Taras  et  al.  22), 
0.22±0.11.  It  was  therefore  concluded  that  the  observed  yield  fluctuations  are  in¬ 
sufficient  to  cause  any  appreciable  error  due  to  using  SCN  model  predictions,  and 
accordingly  they  have  been  used  throughout  the  following  analysis. 

3.3.2.  The  2217,  2488,  2795  and  3274  keV  levels.  Analysis  of  the  measured  angular 
distributions  of  the  ground  state  decays  of  the  2217,  2488,  2795  and  3274  keV  levels 
yield  results  in  agreement  with  the  previous  assignments  of  ]  +  [ref.  2I)],  \  [ref.  *)], 
|+  [ref.  21)]  and  f  [ref.  2)],  respectively.  The  mixing  ratios  obtained  in  the  present 
work  are  also  in  good  agreement  with  the  previous  results  (see  table  4). 

3.3.3.  The  3171  and  3185  keV  levels.  The  yield  of  the  3171  keV  y-ray  was  obtained 
by  subtracting  the  contribution  due  to  the  35CI  3162  keV  y-ray  using  a  Gaussian  line- 
shape.  It  was  possible  to  make  a  unique  J  =  \  assignment  to  the  3171  keV  level.  Al¬ 
though  it  is  not  possible  to  completely  rule  out  negative  parity,  due  to  the  large  mixing- 
ratio  error,  positive  parity  is  more  likely. 
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The  angular  distribution  of  the  1574  keV  transition  from  the  3185  keV  level  (fig.  5) 
is  consistent  only  with  J  =  §  or  Jj1-.  Furthermore,  using  the  measured  lifetime  and 
mixing  ratios,  it  is  possible  to  assign  J*  =  to  the  3185  keV  level. 


ARCTAN  8 


Fig.  5.  Gamma-ray  angular  distribution  (insert)  and  % *  plots  for  the  1374  keV  transition  from  the 
3183  keV  level  to  the  1611  keV  level  measured  at  a  proton  energy  of  3.13  MeV.  The  solid  and 
dotted  lines  indicate  the  best  fit  for  various  spin  and  mixing-ratio  combinations. 


4.  Discussion 

Several  DSA  lifetime  determinations  exhibit  a  large  variation  between  different 
experiments.  In  particular,  in  37Ar  the  recent  results  of  Caraca  et  al.  23)  show  a  dif¬ 
ference  from  other  work  of  up  to  a  factor  of  10  (see  table  1 ),  which  cannot  be  account¬ 
ed  for  by  the  normally  expected  uncertainties.  They  concluded  that  the  only  possible 
reason  for  the  observed  discrepancy  is  that  the  nuclear  stopping-power  prediction  is 
strongly  underestimated  by  the  Lindhard  theory  for  light  ions  in  high-Z  solids  below 
certain  critical  velocities.  However,  it  is  not  clear  that  this  is  the  only  possible  explana- 
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tion.  It  is  instructive  to  look  at  the  difference  in  the  quoted  lifetimes  for  the  1410  and 
2488  keV  levels  between  the  results  of  Caraca  et  al.  and  both  those  of  Ragan  et  al. 
and  the  present  results.  In  each  case  a  time  difference  of  approximately  650  fs  is  indi¬ 
cated.  This  time  difference  would  correspond  to  the  recoiling  nuclei  travelling  ap¬ 
proximately  0.5  pm  (assuming  p  =  2.8  x  I0-3)  before  starting  the  slowing-down 
process  described  by  the  F(x)  calculation.  As  Caraca  et  al.  used  a  thin  target,  thick 
backing  configuration,  a  mean  target-backing  gap  of  about  0.5  pm  would  also  explain 
the  quoted  differences. 

Two  shell-model  transition-rate  calculations  for  positive-parity  states  in  37Ar  up  to 
2796  keV  excitation  have  been  reported  2I).  The  calculations,  of  Dieperink  and  of 
Wildenthal,  use  the  whole  2s-ld  shell  configuration  space. 

The  1410  keV  level  Ml  lifetime  predicted  by  Dieperink  and  Wildenthal  are  730  fs 
and  470  fs,  respectively.  Both  are  faster  than  the  present  experimental  value  of 
940l|oo  fs.  Similarly  a  lifetime  of  «  3.5  fs  is  predicted  for  the  2795  keV  level  by  both 
calculations,  compared  with  the  experimental  value  22 ±14  fs.  The  absolute  values  of 
the  mixing  ratio  prediction  of  Dieperink,  |<5|  =  0.07  and  Wildenthal,  |<5|  =  0.12  for 
the  ground  state  transition  are  to  be  compared  with  the  experimental  value  —0.21  ± 
0.13. 

The  E2  lifetime  of  the  2217  keV  level  was  calculated  by  Wildenthal  to  be  420  fs,  in 
reasonable  agreement  with  the  experimental  value  t  =  5 1 0 ±  70  fs.  This  level  was  not 
treated  by  Dieperink. 

The  161 1  keV  level  ( J"  =  }~)  has  most  of  the  lf^  strength  in  37Ar  and  is  most 
likely  of  a  single-particle  nature.  Harris  24)  has  calculated  the  theoretical  M2  lifetime 
and  E3/M2  mixing  ratio  using  a  1  f^  configuration.  The  predicted  lifetime  of  12  ns  is 
to  be  compared  with  the  measured  value  2S)  of  6.5 ±0.3  ns.  The  predicted  mixing 
ratio  is  0.087,  in  fair  agreement  with  the  experimental  value  of  0.1 6 ±0.04.  However 
it  should  be  pointed  out  that  correcting  the  discrepancy  with  the  M2  lifetime  would 
result  in  worsening  the  agreement  for  the  mixing  ratio. 

Note  added  in  proof:  The  lifetime  of  the  3185  keV  level,  among  others,  has 
recently  been  measured  by  Luketina  et  al.  26).  They  obtained  an  average  lifetime  of 
t  =  300±60  fs,  in  good  agreement  with  the  present  value. 
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Appendix  F 


Listing  of  ' OXDS1  -  a  program 
to  calculate  F(t) 
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r" 

OS  ALGOL  F 

SC 

SOURCE  PROGRAM 

SOURCE  STATEMENT 

> _ 10 .00  0. 

•  nrc;  ?h  • 

00000 
00  0  00 

coooo 

•REAL*  A1A,A2A.ERR1  ,  t*  ,  PH , V  M 1 N , CELT  A , E FJR? , 0 C • DUW . AC .  INT2.AV,ETAU,A1 
A2 , A3  .  A4  ,  VMAX  .  IM I  N  *  AO  »  A  A 1  , ACCC , PPP , A A A .BBB.FUDGNU . FUO C E L , NNE N « PATH 
RANGE ; 

♦ 

♦ 

00001 

0000  1 

coon? 

•  I NTEGER *  NSS  ,NEN,  I  ,NUCLE 4R, NR,  N,  TK» COUNT T , J,K . F 1  ,F2,F10.F , ZI  ,  II . 
R.NL  ,  PP.  STEPS  ,SUP  ,NVMAX,  I  KL ; 

•PEAL*  •  ARRAY*  0  ,PFO  ,  A  ,  Z  (  /  1  :  6/  )  ,  EXN  .  RES  ,  GOL  ,  COL  1(  /(•:  in/lt 

00002 

0000? 

OOOOJ 

E ,DE , TAU(  /o: 40/ )  , VMAXX ( / 1  150/ )  ,EXFC. ShAPE.L IM,T1ME,V  INT , 
oec.velefe(/-i:201/) ,cni ( /o : io, o: io/i ; 

•  PROCECURE*  V1R  ITE  (  OV  ,F  ,M  ;•  VALUE  •  OV  ,F  ,  N  ;•  INTEGER  •  OV  ,F  :  ’REAL  •  N ;  *COOF 

• ; 

0  00  OH 

OOCl  3 

0001  6 

•  PROCEOURE*  vi'R  ITE  TEXT  (  DV,  S  );  'VALUE  *  07  ;  *ST«  I  NG  •  s;  •  INTEGER  *DV  ;  ’CODE  *  *, 

•  INTEGER*  • FKOC  FCURE*  FCPMAT (S) ;  ‘STRING  *S;  'CODE  •  ; 

0  0  0  16 

00017 

000  17 

•  PROCEDURE*  C CP Y ; 

•  BEGIN* 

•  T  NTEGFR 'BEG, I  ,j;  •  INTEGER*  *ARFAY*CCFY(/1 130/)  ; 

0  0  019 

0  0  0?  1 

0  002? 

SEARCH  :  I  MSYMBOL  <  0  ,  •  (  •*')*,  REG  )  5  *  I  F  •  BE  G->=  1  *  THE  N**G0TU*SEARCH; 

i  :  =  o; 

I  npijt  :  i  :  =  if  i  ; 

00023 

00024 

0002^ 

I  NS Y MB CL  <0  , *  ( *  «  A  EC  CEE  GH I JKLMNOPGRSTUVWXY  ZC l 23 4 56 7 6 9 (  )/,,*•  ) •  , COPY!  /I  / )  )  ; 

•  IF  •  COPY  (/!/)-,=  1  •  THEN  »  ‘GOT  C*  I  NPL’J 

0002  5 

o  o  02  r> 

00025 

*  FUR • J := 1 «  STEP*  1*  UNT I L  *  I - l  *  DC • 

OUTSYMEOL  (  1  ,  M  ' * ADCDEF  GH I  JKLMNCFCRSTUVWXYZO  122456789(  )/.,*•)*, 
COPY!  / J /  )  )  ; 

00026 

00  02  7 

0  0  02  6' 

Y’RI  TE  TE  X  T (  I  ,*(•»(»  C •)••)')  ; 

*  END •  ; 

0  0  02  r. 

00029 

0  0  0  3  1 

•REAL*  *  PROCEDURE*  1. 1  NOEL  {  Zl  ,  Z2  ,  A2  ,  E  )  ; 

•VALUE*  Z1,Z2,A2,E;  'REAL* Zl , Z2 , A? ,E I 

•DEGIN*  *  CCMPENT  •  E  IN  ME V/NUCLE C N , L I NDEL=DE/OX  IN  MEV-CM2/M.G;  - 

000  31 

0  0  03? 
000  33 

LI  KpEL:=73. 0*  (  Zl  **  (  7/6)  )  7  2  2  *SCR  T  (  F  1/A2/I  (  Z  1  (  2/3  )  +  22 **  (  2/ .2  )  )  ** { 3/2 )  )  ; 

•END • ; 

00  0  53 

00035 

0  0  036 

•  RE  AL  *  •  PROCEDURE  *  MK  {  7. 1  ,  Z  2  *  A 1  ,  A2  »  V  )  ;  '  VALUE  *  Z'l  ,Z2»A  1  ,A2,V  5 

•HEAL*  Zl ,Z2 , A1 , A2,v; 

•BEGIN*  'PEAL*  f  P  ,K  ,R  ,  AL  ,  RL  ,  I  U,  1 1. ,  CF  ,  NL  ,  AES  ,  EE  ; 

00037 

0  0  036 

00036 

EE := ( 465 ,5*V**2  )/( 9 .0  •  4 )  ; 

EP:= {3. 2 4*1. 0'4*A1*A 2* EE )/(SCRT<21**0. 6667 +Z 2* *0.6 66 7)*Z1*Z2*{A1+A 
NL { RH/ { A2* 1 . Ob ) ) 7 1 .0*3; 

2  > )  ; 

00040 

0004  1 

00042 

AFS: =1 .096/( Z 1 **0.6667+Z2**C. 6667) J 

r:=ai /A2; 

al : = ( 1 +« ) /2;be :=( AL**2) /r; 

0  0  044 

00045 

0  004  & 

k:=o.o; 

•IF*  E  P<=  0. 38CC01  'THEN*  *  EEG I N  *  I L :  =  2 . 1 4 1 50  ; ’GOTO*  E  2  E  ;  •  END •  ; 
•IF*  FP<=3.  00001  'THEN*  *  H  E  G  I  N ’  1U:  =  3. 14159;  * GO  TO » SA ;  •  END '  ; 

dtsTAAL  sro«fcs 


jj u  :n  as.1-. 


SC  SOURCE  STATEMENT 


SOURCE 


PRO  GRAM 


0  005 1 
00053 
00054 
0  0  055 
00056 
00  05  7 
0  0057 
00058 

00  05  9  _ _____ _ 

00  0  ' 

0  0  'REAL*  •PROCEDURE*  RARCSIN(X);  'VALUE'X;  *REAL*X; 

0  0  063 ’  DEGIN*  *  fE  *  (  1-X**2  )  <  )  .0  »  -  \  Q  ♦THEN*  *£551  N  *  R  ARCS  IN  :  :_=_3 . 1  4  15  5/2; _ 

00064  'GCTCNJfl;  'END'  ;  R  A  R  CS  I  N  :  =  AR  C  T  AN  (  (  X/ S  C  K  T  (  1  -X  *  *  2  )  )  }  ; 

0.0 067  Njs:  'END*; 

0  0  368 _ _ _ _ _ . _ 

00C68  'REAL*  'PROCEDURE*  I  I l ( P , A L ,  AU ) ;  ’  V  ALU E • R » AL . AU ;  » RE AL * R .  AL . AU  ; 

00071  'BEGIN*  ‘REAL*  ALF.BE.PHL.FHU; 

00072 _ alf:  -(  i  »R)/2  ;  be: =alf**2/e  ; _ _ _ 

C0074  PHL: =RARCSI N ( (SI N( AL/2 J/SCRT ( BE ) ) ) ; 

0  0  075  PHU: ~h ARCS  INI  ( SIN t AU/2 i/SORT ( EE)  )  )  ; 

_0_0  C  76  IM:-<-i/(2*(SIN(AU/2))**?)+l/(2MSIN(AL/?))**2)) _ _ _ 

00076  - ( ALF  7 ( 2  *BF*SCF  T (R  *0C )  )  ) *(  { -CCS { PHU ) / (  ( S I N ( PHU  )  )  *  *2  )  ) 

0  O  C  7  7?  +LN( 1 /SIN  (PHU  )-CC5(  PHU  )/SIN(FHU)  )  ) 

00C76_  -LM  (CCS(FHU.-?)/SIN(  PHU/2  )  )  ) /S  CR  T  (  R  ♦  RE  ) _ _ _ _  _ _ 

0  0076  +  <  ALF/ <  2*DE*SGRT (R*BE  >  )) *(  ( -COS (PHL )/(  ( S IN (PHL)  )**2>  ) 

00  0  76  +l_M  1  /  SI  N  (  PHL  )  -  CCS  (  PHL  )/S  1M  PEL  )  )  ) 

0  0  0  76  fLM  (COS  (  PHL/ 2  )  /  SIM  PHL/?  )  ))/SCRT(R*9E); _ _ _ _ 

00  077  ‘•END’; 

0007  8 

00078  •  REAL'  *  PROCECUPE*  I  I  2  (  R  ,  A  L  .  AU  )  ;  »VALUF»  R  »  AL  ■»  AU  ;  'PEAL'  R.AL,AU; _ _ _ 

00  0  3 i  'BEGIN*  'REAL*  ALF  •  BE .PHL .PHU ; 

0CCL82  ALF  :  =(  1ER) /2;  3E:=ALF**2/G; 

00 064  PHU: ARCS  IN( S  IK ( AU/2 )/SQRT ( EE  )  )  ; _ 

OQO^b  PHL.  :  ARCS  I  N  (  (  S  I  N  (  AL/2  l/SCFT  (  EE  )  )  )  ; 

0  0  086  I  12 :=LM  SI N( AU/2) /S I M  AL/2  >  ) - (1 /SCRT { R  >  ) * (SORT  ( AES { BE- (S IN( AU/2  .C )>  **2>  ) 

00  086_. _ _ _ _ _ _ _ _ 

COO  36  -SCR  T ( AES (BE- ( 3IN( AL/2  .0  )  )  **2  )  )  ) 

00036  .  -( ALF /SCRT ( BE*R)  )*LN(  ( 1/SI N (PEU )-COS (PHU )/S IN( PHU )  )/ (  1/S  I N( PHL ) 

AQCAo  —CO  S  (  P  HI  )/SIN(PHL)  ));  _ , _ ’ _ 

00087  'END*; 

0  0  088 

00  0  86 _ '  R  FA  L*  'PRUCFD'JRF'  II3(R.AL);«  VAlUE'P.AU;  »REAL'R,AU; _ 

00091  'BEGIN'  'REAL*  ALF.DE.FHU; 

00002  ALF :=( 1 +R )/2 .0;  B E : - ALF* * 2/R ; 

0  0  094  PHU:  =  RARCSIN  (  (  S  I  N  (  AU/ 2 . 0  )  /  S  C  R  T  (  EF  )  )  )  ; _ _ _ 

0009  5  I  I  3I=-1  .  C/SI  N  (  AU/2.0  )  +  (  ALF  /  (  BE*SCR  T  (  R  )  )  )  * (SORT  (AES(BE-(SIN(AU/2.0))**2)> 

00095 

OOCR5  /SI N(AU/2.0  > )  +  PHU/SCRT (R  ) ; _ ■ _  _ 


RL  :  =  2*ARCT  AN(  (  3/S0RT( E F* *2 -« )  ) )  ;  I  U :  =  RL  ; 

sa:i l  :=2* arct  an ( ( c  . 3 q /sort ((ef* *2 -0.1444  ) ) ) ) ; 

'IF*  FP<=  urincci  1  T  HCH  >  'GOTfl  '  e>  :  _ _ _ 

K:=K+I  II  {  R.RL.3.  14  159)/(r.P*v2  )  ; 

b:k:=k  +  o.£T14*(  i  i  3  (  r  ,  iu)-i  1 3 ( k  .  il  )  )/ep-o.2*ii2(r.  il,  iu  ) ; 

B2B  :  *  IF*  EP**2  *1 . 0*4— 1>1.0»— 2Q  •  THEN  •  _ _ _ 

K;=K  +  0.87*(  I  1 3  (  F  ,  IL)-II3(R,2*ARCTAM  l/SCfiT(EF>**2*l.C*4-i  )  )  )  )/fp; 
mk:=(k*3.14i 59*nl*afs )/ 2 ; 

■end'  ; _ 


»to*£3  mi-565 


- 

. 
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00  096  • END  »  ; 

00-097 

0  C  0  )  7 _ '  REAL*  ♦ PROOF CURF *  U  C  E  CX  (  7  1  .  7  ?  .  A  1  ,  A  2  ♦  EE  )  !  1  V  ALU?  *  F  1-  .  7  1  ,  7  ?  ,  A  1  .  A  ->  : _ 

00099  'INTEGER*  7.1.Z2;  ‘REAL*  A1.A2.EE; 

0  0  X  04  *  HF  GIN*  'REAL*  EP.AE.IDJ 

00  JO  2 _  CP  :=  ( -i  .,?4  »  4»A  1  *A?*EE  )  /  (  SORT  (  (7  1  .66fc74Z?^»0 .6667)  )^?1«Z?MA  1 +  A  2  )  )  ; _ 

"  0  i  103  "  I  O:  =  (Zl-S'Z2*Ai*i).475>/(1.66:»A24(Al+A2)*S0RT<Zi*40.6be7  4Z24  4C  .666  7))  ; 

00104  'IF*  EPC14  • TEEN  *  'GOTO'RRR; 

_ 0  C  I  OS  "  IF  »Ff><  ’  >  'THEN*  »GCTC»  Ra; _ _ _ 

OOiOi  NDEPX  := { i  l D*LN(2 .259*EP)  )/2.0*EP )  ; 

00107  '  GOT  0  *  fim; 

©01  oh _ ra  :ae  :=sqrx (ep)  ; _ _____ _ _ _ 

00109  NDEDX 1=10*0 .25* ( (LN( 2.259*EP )* (EP- 1 4) ) / ( 2. C*SP) +( 13-EP) 

001  09  *<  1 . 97  4  AE*EXP<  -  l.82*AE')  +  0.9225*AE)  )  ; 

_ 0 0  1  1C-  «GCT>)»  FINN _ 

ooui  prr;cp:=sgrt  (ep)  ; 

0  0 1 1 2  AH : -{  *  IF  *EP> (2.2-SCRT  (6.574.1  )  )  ’  AND*  EF< ( 2  .2  +  SORT { 6 . 5/4 . 1  )  ) 

00 1 1 ? _ *  THEN*  (0.06  5—  0  »  C  4  1«(  ( EP-2 . 2 ) *  4 2 )  )  *  E  L  £E  *  0.0); _ 

00  1)  3  NDED X 1 -  C  ;  * ( < 1  • 97  *  EXP ( - 1  .3  24EP  )  *£P  +  0 .0225  4EP )  +  AE)  ; 

00114  FI NI  :  ’END*  J 

_ 0  CM  15 _ _ __ _ 

00115  'PROCEDURE*  P R  I  NT  ( ACO ♦ ARES . AOE G )  i  'VALUE*  ACO , AR ES» ADEG J •  I NTEGER • AOEG  ; 
00113  • ARRAY • ACO .ARES i 

00  1  19 _ •  DEG  IN  *  *  INTEGER  '  AT  ,  A)  .90  C _ _ _ _ _ 

00120  WP’l  TE  TEXT  {  1  ,*(••(•  CC*  >  *  CEGKEE  444**  CUE  FF  IC  IEN  T4PLUS*RES  I  DUE  *  (  •  CC  *  )  *  *  )  *  ) 

00121  R0:=FUPMAT(*  (•-D.DCC'-ND*)* )  J 

00122  WPITF  TEX  T(  1  .»(««(»  ;s  *)»*)*)  ; _ _ 

00123  *r-OR*AI  :  =0  •  STEP*  1  *  UNTIL*  ACEG  *CC*«6EGIN*WRITE<1  »  F  ?_  *  A  I  )  * 

00124  «FOR* A JI =0 • STEP*  1  * UNT IL *  A  I  • DO •* BEGIN • WR ITE  TE X T (  1  ,'(••(* 2S *)*•)•)  J 

0  C  1  ?  5 WRIT  E  (I  *  RO.  ACO  (/A  I  .  AJ/))1 _ 

00  126  *  END  *  1 

C0127  WRITE  7EXT< 1 ,*(••(• 2S »)••)*)  ;  h R I TE ( 1  . PO  ,  ARE S ( / A  1 / )  )  ; 

_ 00129  WRITE  TEX  f  <1  ,*(*»(  *C*  )««)»)  ; _ 

00130  » END* ; *ENC*UF  PRINT; 

0013? 

00  13? • RE  AL  *  1  PRO CEO  UR F •  E  X  P  1  NT  (LC.UF.TAU)  1  *  VALUE' LO. UP. T  AU  ;  'REAL'LC,  UP,  T AU! 

0.01  35  *  B  E  G  in’*  EXPIN  T  :=-TAU4{  EXP{  -LP/  T  AU)  -E  XP  (  -  LO/TAU)  )  ; 

GQ1.36  ‘.END*; 

0  01  37 _ _ _ _ _ _ _ _ _ _ _ _ 

00  13  7  *  Rfc  AL '  ' PROCEDURE • PGL Y(X.A.M.N); 

00  1  3  6  *  VALUE  *x  ,  n.  m;  »  RE  AL  *  x ;  *  arrava;  •  INTEGEP'N.M; 

C  0  1  4  2_  MT  E  G I N 1  •  IMTEGFR*  I  ;  »REAL*P;P:=A(/N/); _ 

00145  FOR •  I  : =N-1  • STEP*- 1 • UNT I L'  M’DC*F:=P4X+A{/I/)  ; 

ooi46  poly:=p 

_ 0  0  14  6  'END'  ; _ _ _ _ _ _ _ _ _ 

0014  7 

00147  • PROCEDURE  *C UR VEFIT(N)  POINTS:  (X.Y)NAXCE GREET  (DEG)COEFFIC IENTS:  ( C ) 

OP14Y  RFSIDUALS:  (RES) ;  ' VALUE ■ N .X ,V  .DEG ;  » I NTEGER'N . DEG  ;  » APR  A Y*  X  .  Y  .C  ,Pc s; 
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0  0  1  51 

0.0  153 

0.0  1  53 

•REG  IN • » I NTEGER*  1  ,  J  *P; ‘REAL • SUN  X,  SUMY  ; 

•array*rc/i:n/)  ,cc.d*m(/o:dec/)  ♦t(/o:ofg»i  :n/  )  ,3(/o:deg»o:deg/), 

PI  .  R0  (  /)  :  DEC.  /  )  ;  SUM  /  :  =  sumy  :  =o.  o  ; 

•C  O  l  5  5 

00  1  56 

001  57 

•FOR*  I  :  =  1  *  STEP*  1  ' UNT  I  L  *  N • CC •  *  E EG l N • SU N X : = SUM  X  +  X  C /  I  / )  ; 

SUM y:  =  sumy  +  y(/j/  ); 

•END* ; 

00156 

0  0155 
C01  6.0 

•FOR  *  P  :  =  0  •  TEP*  1  • UNT  I  L • CEG • DC*  • F  CR  *  J : =P+ 1  • STEP  *  1  •  UNT  IL  *DEG  *DO  * 

c(/p, j/>:=o.o;  res< /o/ ) : =o .o ; 

m(  /o/i  :  =  n ;  D(  /  o/):  3-smiwy  ; c(  /  o ,  o  /  >  :  =  c  r  ( /o  /  )  :  =  D(  /  o/ )  /mc  /o  / )  ? 

0  0  1  63 
00165 
CGI  66 

p(/o,o/)  :=R(  / 1 . 1/)  :  =  i ;  r  (/i  ,o/  ) :  =  -su  /n  ; 

•FOR  •  I  :=  1  •  STEP*  1  'UNTIU'N'OD'  'BEGIN*  s'i/O*  I/>  :=1  ; 

T(  /  l  ,  I  /  )  :-X(/I/)+R(/l  .0/1; 

00  167 

C  0  1  6  8 

0  016  0 

RES( /O/ ) : =RES( /O/ )+ ( Y (/ 1/ 1-C (/O .0/ )) **2 «END*  ; 

•rOR'PI— 1  'ST EP •  1 'UNT  IL*  DE  G • D  0 • 

•BEGlMM!fc'S(/p/)  :=d(/p/)  :  =  m(/p/):=p i (/f/): -fg(/p/ ):=o.o ; 

OO  160 

00  1  o9 

0.0170 

•FOR  *  I !  =  1 • STEP • 1  »UN  ri L'N'DO* 

•  BEG  I N*  M ( /P/ )  : - M  C/P/ >+T(/F,I/)4TC/P«I/); 

D ( /p / ) : =D ( sp /) +y( / 1/ i^T ( /p ,i/) ; 

QO  1  71 

0017! 

0  01  72 

•IF*  P-»=DEG  *  THEN  ’  •  BEG  1  N  •  RQ  {  /P/  )  :=20(/P/)Dt(/I/)*T(/P,  I/) 

*t(/p-i  ,i/);fi(/p/):=pk/f/)4xi/i/)=»t(/p,i/)*t{/p,i/) 

»  end* ; • end* ; 

00  1  7-4 
00175 

00  1  7  6 

c ; /p ,p/i :~cc(/n/ > : =d ( /p/ )/m  c /p/ ) ; 

•  I/--  •  P~.=OEG*  THEN*  *3EGIN«RC{/F/):=P0(/P/)/M(/P-l/)  ; 

P  I  {  /P/  i:  =p  n  /p/  )/m(  /p/  )  ; 

0  017  7 

00  1  77 

OO  1.76- 

•  FOR  *  I  := 1  • STEP  *  1  *  UNT I L  *  N • OC • 

t  c /p+i ,[/>:  =  (x(/i/)-pi i/p/ n*r(/p. i/)-kc(/p/)*t(/p-i, i/); 
a  c /p  +  i  ,o/ ) :=-P! (/p/»iH(/p,o/)-RO(/p/)*H( /«- 1 , 0/ ) ; 

001  79 

00  17° . 

00?  eo 

•FCVR  •  J  :  =  r>  •  STEP*  1  •  UNT  I  L*  P-2  *oc  • 

D(/P  M  , J+3  / )  :=R( /P»  J/ )-PI  C/P/ ) *E ( /p»  J  +  1/ )-EG { /P/ )*e(/P-l,J+l/i; 
ch/ph  ,p/|  :=B(7F,p-i/)-Pi(/P/)  ;e(/p  +  i.p+i/):  =  i*fno  •; 

QC  lW 

0  C  16  3 

0  C  i  6  3 

•  FOR • J 1  =  0 •  STOP  *  1  '  UNTIL* P— 1 '00*  C(/i’>J/)  I  C  (  /P  —  1  »  J  /  )  +  CC  i/p/  )*B(/P,J/>; 

•  FOR  *  1  : =1  •  STEP  *1  • UNTI L ’ N • DO • 

•BEGIN *H{/1/)  :=Y( /  1/ ) ; 

00 16  A 

0  0164 
001c  5 

•FOR • j:=0 • STEP*  1 • UNTIL  *P • CC 'R (/ 1/ ) :=R{/ 1/ )-C( /P*  u/) 4 
(  •  I F  *  J  =f-  »  THEN  •  1  *  ELSE  •  X  (  /  1/  )  *  3  J  1  ; 

RE  SC  /P/)  :=RES(/P/)+M/I/)vRC  /  1/  )  •  END*  J 

0  0  1 

C  DC  86 

OO  1  tVH 

•END’  ;  ’ END*  CUP VEF I T ; 

*  PROCEDURE*  ST  MRIJLE  (  4  »  R  ♦  F  ♦  INTEGRAL  ♦  C  «  ACC  *  F  A  I  L  )  ; 

00  J  09 
0C1  91 

cei4M 

•  VALUE •  O.O.ACC;  • RE AL •  A ,H ♦  JNTEGR AL , C , ACC ; 

•LAB EL ‘FAIL; 'REAL* ’PRCCFCURE'F; 

•  PEG  I N» •  J  NTFGER*  M  »  N . NN I  *  RE  AL*  C  »  H. CO , G1  »  G2  »  G3  »  G4 , AO , A l . A  2  »  RES  »  X  ; 

001  95 

0  01  99 

0  0  202 

•ARPAY'S! *  s  2  *  S  3 ( / 1 :i5/);res:  =  c:  =  o.o:m:=n :=?;NN:-ib; 
go:=fca) ;G2:=fc (ahu/2.0) ;ga:=fcb) ; 

AO:=CE-A)*(GO+4.04G21G4)/2.0; 

00203 
00206 
.Jn  oajj 

aa :d :=2.o**n;h :  =  <  b-a) / ( 4.o*d ) ; x:=a  +  ( 4.o*m  +  i . o > *h ; 

G i :=f ( X) ;G3:=f{x+2.04h); 

A1  ;-H'?(G0  +  4.0’>Gl  •»  G  2  )  IA?:=H4(G2  +  4.0^  G  3  +  G  4  )  J 

<T  *.  N  T  A  AL  S  IO*tS  Mr 
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0021  0 
Q021  0 

fi&zxz 

•  IF*  APS<  A1  +A2-A0  }>ACC/C  1 T  FEN  * 

•  DEG  1 n*m:=m+m;n:  =  k+ i ; 

•  IF  *  N >  N M  '  THEN*  '  GOTO*  FAll  ; 

0  0'2  1  5 
Q‘?  3*«> 
00^2  ' 

ao  :•=  a  i  ;  s  i  (/n/)  :-a2;s2(/n/)  :=G2;s3(/n/)  :  =  G4  ;  64 
•  GOTO • A A ;  •  END •  *  ELSE*  • BEG IN • 

RE  S  :  =RE S  ♦-  <  A  1  +  A2  )  /  3  .  o  ;  ^ =  y  +  1  ;  c  :  =  A  +M *  (  E  -  A  )  /D  ; 

:=G2 ; G2 : 

~gi  ; 

o  ■ 

Q- 

0  0  22  f 

BB  :  *  IF  *M  =  2* ( M* / *  2)  *  THEN’ 

*  OEG I N • M : =  M • /  •  2 ; N : =  N - l ; • GO  TO  *  E  e ;  • E  ND ■ ; 

•  IF  *  M->=  1  *  OR  •N-.-4  'THEN* 

00  227 
0022  H 
0023  1 

•  REG  I N • A  p : =S 1  { /  N  /  )  ; 

G0:=G4  ;<  2  :=S2(  /n/  )  ;ci4:  =S3(  /n/)  ; 

•  GOTO  •  AA  ,  ’END4  ;  'END*  ; 

00  2  34 
00236 

00236 

mj : integral  :=r es ; • end • ; 

•REAL*  *  PliOCF.DUP  F  *  D  X  C  V  )  I  •  VALUE  1  V?  1  R  F  AL  *  V  ; 

0  0  2  3  r> 

0024  1 

002  42 

•  BEGIN  •  •  REAL  *P0  ,  A  i  L.A2BJ  •  INTEGER  •  <,  ; 

•IF* V<=  (E(/1/))*TFFN*  • GOTC • MA ; 

pol:po  :  =  poly (  v.ccli ,o ,6); 

O'  43 

00244 

0  0  24  5 

•GOTO* EG; 

ma: pc:  =  ( poly ce (7 1 /) ,coli *  o, 6 ) *v  >/e l/i/ ) ; 

EG:  •  IF • NUCLE AP<0 . 5 • THEN  •  •  G C TO • J P 5 ; A  1 B :  =  o ; A 2H : 

=  o; 

0024  3 

0  0^48 
00246 

•FOR'C-l  ‘STEP*  1  MJNT  IL*  NSS'OC* 

•BEGIN* 

Alfi:=A18+Q(/K/)*Af/</)4 

• 

0.024  8 
OO  249 
CO  2  5  0 

NDEDX  (  ZI  ,  Z  (/K/  )  .  ( M/G  .0664  )  ,  A (/K/ ) ,  (5.172'-3.*(V**2>))  ; 
A20:=A20FG</K/)=FA(/K/); 

•eno  • ; 

C-C-2'  j  | 

0  0  252 
00  2  53 

PO:=PQF( 0. 32G4*F UDGNU *RH  + A  IE I/A2B ; 

JR5: dx:=m/po ; 

*  END • OF  DX; 

- 

0  02  5  4 

00  254 

0  0-2.55 

•REAL* •PROCECURE* READ; 

•  BEGIN •  *  «F AL • V ; I NREAL ( 0  « V )  ; PE AD: =V ; 

0  02  5  3. 
00259 

00250 

•  ENO  «  ; 

.  •  REAL  •  •  PROCEDURE*  MKDX  (  V  }  ;  •  VALUE  »V  *,  *  RE  AL  *  V  ; 

« 

• 

00  252 
QG-26  5 
OO  ?b  5 

*  BEG  IN* •  I NTEGER*  El  'REAL*  A1  C  .  A2CIA1  C:=A2Ct=0.0  ; 

•FUR*  E := 1  *  STEP •  1  • UNT I L  *NSS  *  DC  * 

*  BEGIN*  AI C:  =  A1 CPA ( /E/ J*Cl/E/ ) *MK(7 I .Z(  /£/),( M/0 .066  4 ) 

iA( /E/) tV) ; 

0  0  266 
00267 

0  C-2  6  ft 

A2C:=A2C*G)</E/)*A<  /E/  ); 

*  END  * ; 

MKOX  ;=FLOGNU*  C  V’P  OX  (V  >  *  A1  CT  /  A2C  ; 

0026  9 
002  7  C 

CO  2  7  0 

•end*; 

S  Y  SAC  T {  1 S  ,1 2  » 1  ) ;  SYSAC  T ( 1  .e  .1 32  )  ; SYS  ACT ( 1  . 1 2. 

Q0  273 

0  0274 
CO  27  5 

starta:writetext(  i.*<  •  •  ( *cc*  )•  •)•)  ; 
copy; 

ININTEGER (0»NSS);  •  IF»NSS<0.5*TPEN*  *GCTO*FlN: 

CCnT*A\l  sfjm'cs  nV-SdS 
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sc 


SOURCE  PROGRAM 


SOURCE  STATEMENT 


0C277  I NRE AL ( O  .  RH  )  ; I NREAL ( 0,M)  ;  ININTEGERIO  ,ZI  )  ; 

0028  0  ‘FOR  *  I  := 1  'STEP  *  1  ’ UNT  !L* NSS'DO* C(/l  /)  :=REAO; 

0  0  281  *  FOP  *  I  : -  1  '  STEP*  1  'UNTIL  '  NS5  * 03 '  A  (  /  I  /  )  : -READ; 


0028?  'FORM  :=  l  '  STF.f-*  •  1  *UNT  IL»KSS«DC'  2(/I  /)  :=PEA0; 

00283  I  NlNTEGERtO  .NCM  ; 

0  0  28  4  •  1  F  •  NCNCO  .  1  1  Tl-  Ef  J  ’  »  BEGIN  •  NNEN  :  =-  1  0  JNEN  :  =  t  ;  I  NRE  AL  (  O.FUDGFU  ; 

00287  • END* • ELSE'NNENI-IO ; 

00  286  •  IF*  NK'EN  <0*THEN*E(  /  l/)  1=0.1’  ELSE  '  'FCR'  11=1  'STEP*  1  'UNT  IL'NFN'OO' 

C029-S  F  ( /] /)  :=reaii; _ _ _ _  _ 


00  28  0  ‘FOR  •  I  :  =  1  • STEF*  1  • UN T I L  * N EN  • DC  *  • E EG  I N • 

00280  ai a: =0 ;a?a :=o ; 

00  29  1 _  •  FOR  «  I  1  :  =1  •  STEP  •  I  »  UNT  I  L  *N SS  *  CO  1  «  HE G  I  N  1 _ _ _ 

00  201  aia:=aiaf(q(/ii/)*a(/ii/)  +  (»  I F •  N N F N < C*  ' T H E N •  L I M D EL ( Z I ♦  Z ( / I  1  / ) » 

0  0291  A(/II/),0.1  )  *fudgf_l  *flse  »reao>  >  ; 

C  C  29?  A?  a  =A?A»C(/I  17  )*A  (/I  I/); 


0  C293 
C  0  29  4 
00295 
00296 
O  02  97 
00297 
0  0  299 
0  0  20  1 
00  302 
00  302 
C  0302 
0  0  AQ3 
00  3C  5 
SO  30  6 
00  307 
00308 
0  O  3  C  S 
OC  31  1 


'  E  NO  * ; 

DE< / 1/ ) : =A1 A/A2AJ 

• fno • ; 


•  IF  *NEN= 1 • THEN*  *  GEG  IN  »NEN:=7 ; 

’FOR'  I  :  =  2 ' STEP  •  1  'UNTIL • 7 'DC*  'BEGIN' 
t ( / 1 / ) : =E ( / 1 / ) *  I ;PE(/ i/) :=pe (/i  / )*scrt ( i > : 


•  END '  ;  •  END • ; 

F  1  :=FOFM  AT  (  •  {  •  SNCC  .DDDSS*  )*  )  ; 

«  IF  «  NNENCO  »  THE  N«  'BEGIN'  WFITETEXT  (  1  ,  *  (  *  »(  »CC»  )  »L  INCH  ARP AST CPP I NG» »  I  THv-EL  E 
CTRCNI C^FACTOR*=*AA**' ) • ) ; 

WR1  TE  (  1  iF  1  ,  FUDGED  ;  WRIT  ETC  XT  (  1  .♦(''(  1  CC  *  )'»)')  i _ 

’  End  *  ; 


WF  I TFTEXT  c 1  ,'{"(•  c*  ) '  F ( NEV/ANU  )  +ACF/BX (MEV  .CM/MG • (• 2C  »)••)•) J 

F  1  0  :  =•  F  n P M  AM  *  (  *  —SEC.  DOC  '  -NC  *  )  *  ) ; _ _ _ _ 

•FORM  Z  =  l  'STEP*  1  '  UNT  I  L  'NF>  •  DC”  'BEGIN* 

WRI-TE  (  1  *F1  0  ,E(  /I/>  >  J  V»fi  ITE  U  .  F 1  0  *  DE  C  /  I  /  )  )  ;  W  R  I  TETEXT  {  1  ,  •  (  ••(  *C  *  }  •  •  >♦); 
E ( / T /  )  : -300 . P*50RT( F (  /!/  )/4fc5.5)  ; _ 


00  31  ? 
0031  3 
00  31  5 
00  216’ 
0031R 
00  322 


•  END ’  ; 

NP :  =  NE  N;  I NINTEGER (O.NVMAX) ; 

•FOR*  I :  =  1  • STEP  '  1  *  UNT  T  L • N VMA  X • D  E  * VMAXXt /I  /)  :=READ| 


I  NRE  AL  (  O  *  VM  IN  1  ;  I  N  I  NT  EG  F.R  (  0  ,  N  )  ;  I  M  N  TEGE  R  (  0  .  T  K  )  ; 

•  IF  *  TK  <0 .5 ' THEN '  *  BEG  IN • T  AU { /I  / ) : =0 . 002  ;  T AU < / 2/ >  :=0. 004; TAUl 73/ ) :=0 .01  ; 

♦  FOR •  I : =4  »  STEP'  1  ' UNT IL»1F*OG*TAU(/I/):  =  !0*TAU{/I-3/);TK:  =  16; 


C  0  32  4 
0  0  325 
0R32  7 
00329 
00329 
0032  9 


•END*  'ELSE'  ’ FUR •  I  :  =  1 *  STEP*  I  »  UNT IL •  TK • DO • T  AN ( / 1/  )  : =2  *  12  * ( RE  AD  )  J 
INREAL (0 . ACCC)  ;  IN  INTEGEP ( 0. SUP )  ; 

I NRE  AL  (O . F  LDCNO  )  ;  '  IF • FUOCNUC-O  .1  'THEN  'NUCLEAR : = O  •ELSE* N UC  LEAR: = 1  ; 


V  f  TE  TEX  T  (  1  ,  •  (  *  '  (  •  CC*  )  •  M  ASS*  OF*  T  CN*  (  *  SSS*  )•  CHARGE**3?  **DENS  I T  v  *47  iAND*  A** 
HF**6ACK I NG'  (' CC •)•')*) ; 


00330 
CFO  333 
00333 


F  1  t=F  ORM  AT  (  '  (  *  SNDC.UOCSS  '  )  '  )  ;  F2  :  =FCRMAT  C  '  (  *  SSSSNDDSSS  *)');VfRITECl.Fl*M); 
WRITE! 1 .F2.ZI  );WRTTE(1,F!,  FF )  ! _ 


l  ist«lku  sr0Ris’n> -5  3  4 


. 
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sc 

SOURCE  PROGRAM 

SOURCE:  STATEMENT 

0 

94*  3  35 
©.0  3.!  5 
..  C  0  3  i  7 

•FORM  :=  1  •  STEP*  1  •  UNTIL*  NS  S  *  DC*  • BEG  IN • »  IF  •  l  >  1 . 5*  THEN* 
•BEGIN  *WR  I  Tt  TEXT  (  1  ,  •  {  •  •  (  *C  *  I  •  *  >*  )  ;  WKI  TLTKXTI  1  ,  •  (  •  •  (  •  235* 

»  END  *  I  r.RI  TE  ( 1  ,F2 , 7  (  /  [  /  )  )  ;  v.  T  1  TR  (  1  ,  r  1  ,  A  r  /  1  7  )  )  ; 

)••)*); 

00  34  0 
0034  1 

0  0343 

•end*  ; 

WRI  Tf.TEXT(  1  ,*(*•(•  CC*  )  *M  IN*  VCLCC  IT  Y*=  **+•)•);  W«I  T£(  1  .  F10 
•  l  F  *  N  UCL  F  AR  >0  .  )  *  T  F  EN  •  •  9E6  I  N  •  , 

,  (  VM  IN  *  I  o*  *£>  )  ; 

00  34  3 
0034  4 

0.3  34  5 

WUITETEX  T ( 1  ,  •  f  •  • (  *  CC  *  ) • W ITH *NUCLEAP*T AC  TOR*=****  •  )  •  )  ; 
WRITE!  1  ,F1  *  FUOGNU )  ; 

•  E  NO  •  ; 

00  34  6 

0  0  34  7 

C  034  5 

WR  I  TETEXT  (  1  ,*(*•(•  CC  •  )•*.)•  > ; 
m:=m*o.0664; 

•FOR*  I  :  =  1  •  STEP*  1  •  UNT  IL  *NP*  Dll  *OE  (/  1/  )  :  =  C.  32C4  wrh*DE(  /  1/  ) ; 

• 

CO 34  9 
C035  0 

00351 

CURVEFIT(NP,E  *  DE  «  5  »  CO 1  • RES  )  J 

PRI N  T (CO  1  »RHS  ,5)  ; 

•FOP  •  l  :=0  *  STEP •  1  *  UNT I L *  5 • O O • COL  I  (  /I / )  :=CO! ( 75 

.I/;; 

00  353 

00353 
00  3  54 

•FOR*  I JKL: =1  'STEP*  1  ' UNTIL* NVVAX* DC*  'BEGIN*  VM  AX:  =  VM A  XXI / I JKL/) ; 

DELTA I-LN( VMAX7 VM I N ) /NJ 

f  xpoi /o/ ) :  =  3 ; i nt ? :  =  o.o ; ti me< /c  / ) :=o . o ; vi nt (/o/i :=vmax; 

0  0  353 

C  035  9 

0  0  350 

DF-(  /I  /  )  :  =  POLY  ( E  (  /  1  /  )  ,CCL1  ,0,6); 

*  BEGIN • 

path  :  -range :  -~o  .o  ;  vfi.eff  ( /  o/  ) :  =vma  x  ; 

0  C361 

0  0  351 

■  C  0  362 

'FOP  *  I  I  ;  =  1  ‘STEP  *  1  •  UNIT  I  L*  N*  DC*  •  BEGIN* 

vint(/ii/):-vint(/ii-i/)*exp (-delta >  ; 
ac:  =  dx ( v  i n t ( / » 1/  )  ) ; 

00  35  3 

00354 

00  355 

•IF*  SUP< 1  .9  *  THEN  *  ’ COTC*  GS1  ; 

WRITE  TEXT ( l  ,*(**{  *C*  )•»)•)  ; 

write ( i . fio.ac) ; 

0  0  36  H 

0  036  7 

00.36  7 

c  s  l :  •  F  nd  *  ; 

****-* VELEFF***** 

C  0  3 1  7 

0  0363 

0  C36P 

**PATPLEMC>TH  *  *M£  AN  P  AN  C-  F  77  *  (  TEE  LATTER  IN  MG/CM2)  *(  '  C  *  )  *  * 
•FOR*  I  1  :  -  1  • STEP •  I  *  UNT IL*N*DO*  ‘BEGIN* 

VINTf/I!  /',  ;=VINT(/  I  1-17) *EXF(- DELTA)  ; 

)  *  ) ; 

0  0  5)4 

CO 3 6  7 

0  0  3  7  0 

AC:  =  A3S(  (  7 INT ( 7 I  I— 1/)— VINT (71  17)  )  *DX(  (  VI  NT ( /I  I- 1  7) + V INT (  / I  1/ )  ) 72. 0)7 
(4  00  ,C *ACCC ) ) ; 

•  IF • SIP<1 . 9 • THEN*  *  GOTO*  QS2  ; 

0  0  371 
00372 

00  373 

WR I TETEXT ( 1  ,*(•*(  *C*  )•*)•)  ; 

WRI TE (  1  ,F 10  ,AC  )  J 

QS2  :  S  I  MRULE  (  VINT  (  71  17)  ,VINT  (71  T— 1/  )  ,DX,  INT2,Ol)M,  AC, FA  IL1 

i  : 

00374 

00375 

0  0375 

T  IMF  (  /  I  1  7  )  :=T  I  ME(  /  I  1-1/  )  +  I  NT 21 

AC  :  — ABS(  (  VI  NT {/I  1—17)— V  INT  (71  I/)  ) *  N  KDX  (  (VI NT  (/ I  17 ) +V INT (/  l  I— l/))/2)/ 

(AootAccc ) ) : 

00  3  76 

00377 

0  03  78 

S I MRULE (VINT (711/), VINT(/I1-1/;,NKCX,  I  NT 2 ,DUM ,  AC, FAILS); 
E  XPO (/II  /)  :  =  FX  PQ { 7 I I  —  1 7 )  +  I N  T2 ; 

VELEFF (71  17)  :  =  V  INT ( 7 I  I / ) *£ XP ( — F X PC ( 7 I I 7 )  )  ; 

00379 

00300 

00  3  00 

PATH  :=PATH*-O.S*<  V  I  NT  (  7  1  !-l/)+VINT(/I  17)  )*(TIHE(7I  I  /  )  —  T  I  ME  (71  1-17)  >; 

RANGE: =R ANGE  +  O .5*(VELEFF(/I  I- 1/  )+VELFFF (7  I  1/ )  ) 

*(TIME(7II/)-TIME(/II-l/>); 

V<  •  T  h  *L  SI  t  K[  J  Hf-SaS 


. 


SC  SOURCE  STATEMENT 


SOURCE  PROGRAM 


0  0  381  W«I  TFTCXTd  (  'C*  )••  )'  ) ; 

00362  WP  l  Tlf  tUF10.(V  I  NT  (  /  I  I  -  1  /  )  *|  0  *  *8  )  )  ; 

O  a  <83  W  P  I  TF;  <  1  ,  F  1  O  ,  (  T  I  V  E  (  /  I  I  -  1  /  )  »  5  .  o  /  (  1  «  1  3  )  )  )  ; _ 

0<r  384  Vi R 1  TE  {  3  ,F  1  0  ,  E  XP  (  -E  XPO  (  /  I  I  -1  /  )  )  )  ; 

C  0  3  3  5  WR  I  TF  (  1  *  F  1  O  »  V  L’  LuFF  (  /  I  1-1/)  *1  0**3)  ;  VrR  I  TF  (  1  ,FIO  ,PATHvPH/20)  ; 

093  8  7  VT-U  lr(  1  .FI  O  tFANoE*PH/20)  ? _ _  _ 

00368  'END*; 

OC  >80  V.’R  I  TF  TEX  T  (  1  »•(  ••(  •  CC  *  )  'L  IFF.T  I  yE*******FTAU******v*CUTEr?R*****I  NTF^S'  (  •  C 
■TO  360  ; *  i  *  * )  * )  ; _ _ _ : _ 

GO 30  0  • FOP •  I  I  := 1  » STEP*  1  ’UNTIL*  TK  ‘DC*  ' tFG  IN  •  AC : =HC: =C  .  ? ; 

00  391  •  FOP*  I  :  =  0  *  STEP  •  1  •  UNT  IL'N’DC'  DF.C  (/I/)  :  =  EXP  ( -T  I  ME  (  /  1  /  ) /T  AU  (  /  I  I  /  )  )  ; 

O  0  2  _ *F  OP  *  I  :  -  1  *  ST  EF  '  I  'LNTIL'N  *  Ct* _ ■ 

0  0  3<  2  'BEGIN*'  IF' V£LEFF(  /  I-  1/HOEC.  (  /  I- I  /  )<  I  '-6*  VMA  X  'ANO*  I  >2'  THE  N  *  ’GOTO*  TERM; 

O0  3P3  AAA;  =  (VELEFF(/I-1/  )*TIME(/I/)-VELEFF(/I/)*TIME(/I-l/)  ) / c  T  I  ME( / I / ) 


C  O  30  3 
0  0304 
OO  .30  S 
0  0  33  5 
0  0.  *6 
0  03 '  »  7 
OS  3Q  3 

CC399 
C  0400 
004  0  ! 
o  o4o;~ 
0>?  40  J 
-e.64 ?  .• 
0t;4n7 
0  0'4  O  H 
C0409 
C-C4  1  0 
00412 
0.0  4  1  3 


-T  i. ME  (  /  !-!/)  ) 


EB6:  =  (VElEFF(/I/)  -VELEFF  (  ✓  l-  1  /))/.(  TIME  (/  1/ )-T  I  ME  C  /l-  1  /  )  )  ; 
AC  :  =  AC  +  (  A  AA-t-BEB  *(TIME(/I-l/)+TAU(/Il/)  ))  *QEC  <  /I-l  /  ) 

{  A  A  A  +  FEE*  (  T  IHEC  /  \/  )  TT  AU  (/II/)  )  )  *PEC  (/I  /  )  I _ _ _ 

*  END  '  ; 

AC  :=ac/vmax ; 

V.RI  TE  TEX  T  (  1  ,»("(•  C*  )  "  )  '  )  ; 

>*R  1  TE  (  1  ,  F  10.T*>U(/I  I  /  )  *5  '-131; _ 

UPITEC  1  *  F  I  O  i  AC.  )  ! 

V.PITEI  l.Fin.EFFl  )  ;  V.RITEU  ,  FI 0.E9R2  ); 

*  EN D  *  ;  * E N D ’  ; 


WP1  TETFXT  (  1  ,*<••(• 5C ’)*•)•)  ; 
'END' ; 

*  GC  TO  * RSD2 ; 


FA  IL  1  :  ViF  l  TETEX  rd.'CM'C'  J'FAlLl  *  (  *  C  ’)«•)•)  i  *  GOTO  •  PS  0  2  i 
FA  ILS: WPI TE  TEXT (  l  ,  »(  *  FA  IL5 •(  *C  •)•*)»  )  ; 

FSC.2  :  *  GOTO*  START  A  I 


0  0  414 


fin: • eho • ; 


. 
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